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C. 1 HEAT PIPE ANALYSIS AND DESIGN CODE USER’S MANUAL 


C.1.1 Introduction 

Thi3 section describes the utilization of a digital computer code for Heat 
Pipe Analysis and Design (HPAD). Basically, HPAD calculates the steady-state, 
hydrodynamic heat transport capability of a heat pipe with a particular wick configu- 
ration and working fluid as a function of wick cress- sectional area. Heat load, orien- 
tation, operating temperature, and heat pipe geometry are specified. Both one "g” 
and zero "g° environments are considered. At the User’s option, the code will also 
perform a weight analysis and calculate heat pipe temperature drops. Each of the 
following wick configurations whose cross-sections are shown in Figure C. 1-1 can be 
analyzed using HPAD. 

• Central Porous Slab 

• Circumfere*..l*l Porous Wick 

• Arterial Wick 

e Annular Wick 

• Axial Rectuangular Grooves 

Both the composite and homogeneous modes of operation can be evaluated for the first 
three wick types. The analysis and formulation of the equations used in the program 
are presented in Section C. 1.2. The basic closed-form solution for heat transport 
capability is presented, and the wick properties and self-priming requirements are 
established for each of the wick configurations. A weight analysis and a heat transfer 
model are also developed. This section concludes with the Method of Solution and 
User’s options. 

The program input requirements are described in Section C. 1. 3. 2, and 
the output formate are described in Section C. 1. 3. 3. Nomenclature is listed at the 
end of Section C.l. 3. The flow diagrams, program listings, and sample problems 
are presented in the Appendices. A listing of FORTRAN names with engineering 
quantities is also presented as an Appendix. 
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C. 1.2 Analysis 

C. 1. 2. 1 He jt Transport Capability 


The analysis to determine the heat transport capability of a heat pipe con- 
siders both capillary pumping and sonic vapor limits. Closed form solutions are used 
to predict the steady-state heat transport as a function of wick area in the case of the 
capillary pumping limit or the minimum vapor area for the sonic vapor limit. The 
analysis is performed for the conventional one-dimensional heat pipe shown in Figure 
C.i-2. The following assumptions apply: 

• Uniform heat addition and removal at a single evaporator ana 
condenser section. 

• Uniform wick properties and circular cross-section over the 
entire length. 

• Momentum effects are negligible. 

C. 1. 2. 1. 1 Capillary Pumping Limit 

The closed-form solution for the hydrodynamic heat transport capability 
as determined by the capillary pumping limit can be derived as: 


2 K A (1 + V ) cos 6 F. 

Qj = ; L N, 


r L rr 
p eff 
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Although the individual terms have been discussed previously they will be repeated 
here for easy reference. 


1 . 


The parameter rj is defined as the ratio of the sum of all pressure 
differences due to body forces to the available capillary pressure. 


i. e. , 


= 


r D cos 0 

_E + 

2 H cos 0 


r h 

E 

2 H cos B 
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where H, the wicking height factor, is a property of the working 





fluid and is defined as: 


H - 

fl* 
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The parameter represents the ratio of the flow pressure drop in 
the liquid to the sun of the flow pressure drops in liquid and vapor. 


F = i 

1 A Pj + A P y 


32 K A 
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1 + 


&u 2 

1 h,v 


where the factor <t> depends on whether the vapor flow is laminar or 
turbulent. 


, Re < 2200 
» v 


0 75 

0.0031 (Re ) • Re > 2200 

' v 7 v 
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Re = 
v 
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The turbulence coefficient (0. 0031) has been adjusted so that the 
total vapor pressure drop neglecting momentum effect is equal for 
the laminar and turbulence at Re^ = 2200. This assumption gives 
a smooth transition from laminar flow to turbulence. The maximum 
deviation resulting from the above assumption is about 10%. 


The Liquid Transport Factor is defined as: 
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4. The parameter is the effective transport length defined as: 


C. 1-5 



+ L 

a 
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L + L 
e c 


J e£f 


5. The parameters K, r^, v » A w> and A v are the permeability, 
effective pumping radius, vapor hydraulic diameter, and the wick 
and vapor cross-sectional areas, respectively. These parameters 
are defined by the type cf wick material and wick geometry employed. 
Homogeneous and composite type wicks are treated by determining 
the appropriate values of permeability and pumping radius to be used 
in Equations C. 1-1 and C.l-2. Equivalent wick properties for dif- 
ferent wick geometries are discussed in a later section. 


Once the properties of the wick and working fluid have been determined. 
Equation C. 1-1 is used to calculate the maximum heat (Q^) that can be transported 
without exceeding the capillary limit. When the vapor flow is laminar. Equation 
C. 1-1 can be solved explicitly for either Q t or When the flow is turbulent. 
Equation C. 1-1 becomes an implicit relation for Q^. In this case, the Newton-Raphson 
Method is used to calculate the maximum heat transport. 

The heat transport capability is calculated for both zero "g" and one ”g M 

environments by setting equal to zero for the former and calculating ** (Equation 

C. 1-2) for a specified elevation h (equal to -h) in the latter. The variation of heat 

e 

transport with evaporator elevation (dQ^/dh^) is also calculated for the one "g" lami- 
nar case from: 


K A w F 1 N 1 
db e L eff H 
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In order that a particular wick performs to its full capacity in a gravity 
:teld, it must be capable of complete self-priming. At a minimum, this requires 
that in the static condition and at the specified orientation the capillary pumping avail- 
able during priming must be sufficient to overcome any adverse body forces. Thus, 
for a fully saturated wick, the gravity head required for self-priming becomes: 

h = h + D 01-10 

req e 


/ 
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In the analysis, the self-priming capability of a particular wick is determined from: 



2 H cos 6 
r p, sp 
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where r p = effective pumping radius during filling. Thus, in order to have a fully 
saturated wick: 


h > h 
sp req 


C.l-12 


This program will calculate the heat transport capability for only those wick geome- 
tries which satisfy the self-priming requirement (i. e. , Equation C.l-12). As an 
example, consider the case of an arterial wick (Figure C. 1-1 C) adjacent to the tube 
wall (hp = 0). The self-priming requirement for a horizontal heat pipe with this wick 
is: 

h = D C. 1-13 

req a 


During priming the effective pumping radius is: 

D 

a 

r p,sp ~ 2 
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therefore: 


h = 
sp 


4 H cos $ 
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The code will therefore perform the heat transport analysis for increasing values of 
artery diameter until either: 


D a > 2-^H cos d 
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or the sonic vapor limit or the maximum allowable vapor temperature drop is exceeded. 


C. 1-7 



I 



C. 1. 2. 1. 2 Sonic Vapor Limit 


Generally, the heat transport capability of a heat pipe will be determined 
by the capillary pumping limit; however, in those cases where the working fluid is at 
a low vapor pressure, the sonic vapor limit could become dominant. The minimum 
allowable vapor area that can exist without incurring the sonic vapor limit is calcu- 
lated from: 


v, mm 



2 (7+ 1) 

7 R T 
g 
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When performing the heat transport analysis for a specified wick configuration, the 
transport capability calculated from Equation C. 1-1 is a function of increasing wick 
(i. e. , liquid) area until the wick becomes so thicx that the vapor flow area is reduced 
to the point where the sonic vapor limit is reached. The heat transport portion of the 
analysis is then terminated. The analysis is alsc terminated if the wick is no longer 
self-priming or the vapor temperature drop exceeds a specified value (see Section 
C.l.2.2). 


C. 1 . 2. 2 Wick Properties 

The heat transport capability of any of the five wick geometries shown in 
Figure C.l-1 can be determined using HPAD. Both homogeneous and composite modes 
of operation can be evaluated for the first three geometries. The circumferential dis- 
tribution of the liquid is neglected for the central slab and artery designs. The heat 
transport capability is determined using the equations developed in the preceding 
section. Since these equations are general, equivalent wick properties must be de- 
rived for the specific wick designs and operational modes. In particular, equivalent 
properties must be determined for the following: 

K = Permeability - This is either input as a materia) property; 

or, in the case of the arterial, annular, and groove geom- 
etries, It is calculated consistent with the appropriate fric- 
tion factor and the hydraulic diameter of the liquid (D^ ^). 




i 

f 


y 


n 
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K = 


1AI 

2f Re 


j = Liquid hydraulic diameter - This parameter is defined as 
4A 1 


D h,l P 
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w 


h, v 


H 


req 


= Vapor hydraulic diameter - Tnis is defined analogous to D. 

h,l 

- Effective pumping radius - This iB input as a property of the 
wick material or calculated from the wick geometry. 

- Static elevation head that must be supported in one "g' envi- 
ronment if the wick is to be self-priming. This is calculated 
based on wick geometry. 


p,sp 


= Pumping Radius for self-priming. 


The equivalent wick properties are listed in Table C.l-1 for the various w*ck geom- 
etries. 


C. 1, 2. 3 Weight Analysis 

A weight analysis subroutine which can be employed at the user's option 
is included in this program. The weight analysis is based on containment of the in- 
ternal pressure of the working fluid at a specified maximum temperature. The in- 
ternal pressure is input as the saturated vapor pressure when the maximum temper- 
ature is less than the critical temperature. For temperatures exceedii g the critical 
temperature of the working fluid, the internal pressure is calculated fr im either the 
Ideal Gas Law or the Beattie-Bridgman Equation depending on which is indicated in 
the input. 


For an Ideal Gas: 


P = 


m RT 
g max 
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I r 


where "m" is the total fluid inventory required at the operating temperature, and "V *• 
is the total internal void volume. 


with 


m = m. + m 
1 v 


-i- A v i-/ > i*' A w I - 
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m = P V = P A L 
v * v v • v v 
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When the Beattie-Bridgman Equation is used: 


p = — < v + b >--7 
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where n v” is the specific volume 


v = 


m 
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and 


A = A (1 --) 
o v 
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B = B (1 

o v 
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= 


v T 3 
max 
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A q , a, B q> b, and e are constants for the particular working fluid, which must be in- 
put when the Beattie-Bridgman Equation is required. These constants must be input 
consistent with the following dimensions: pressure in atmospheres, volume in liters/ 
gm mole, temperature in degrees K, and R = 0. 08206 atm-liters/gm-mole-°K. These 
are the units generally found in the literature. 



The weight analysis is performed for a single tube wall thickness which is 
the larger of the specified wall thickness or the minimum wall thickness required to 
contain the pressure with a specified safety factor (S). The weight analysis can also 
be performed parametrically as a function of the radius of a spherical storage volume. 
This volume would be attached to the heat pipe to reduce containment pressures and 
subsequently the system weight. The spherical volume shown in Figure C.l-3 is used 
strictly for containmsnt purposes and should not be confused with the storage reser- 
voir used in gas-controlled heat pipes. A spherical volume is used to simplify the 
analysis; and, although it is impractical because of its fabricability, this model does 
give a good indication of the size of a cylindrical reservoir required to minimize the 
system weight. The wall thicknesses required to contain the pressure are determined 
from: 


t _ LR_s 
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and 


6 = 


P R st S 
2 <S 
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The total internal volume of the system is: 


V = V. + V + V + V 
t 1 v st w 
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v t ■ <™ 2 - 0 - •') v L * i ” V 


a i-32 


The total weight of the system is calculated as: 


m = m + m + m, + n* 
t w np st 
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where 


m 


w 


A >w L 
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m hp = TT(t 2 + 2Rt)f hD L 


hp 
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C. 1-12 


. t 




" [a ? iU « ‘ 4 > 3 - a " R « 3 ' 2 n ll st h f ] /*! 
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. ■> •> 1/12 
h '« R , - (R “ - U-) 
st v st 7 
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The weight analysis is performed for the case of no storage volume by set- 
ting R equal to zero (0). Otherwise, the analysis is performed for increasing values 
si 

of storage radius. The analysis is terminated when either the thickness of the tube 
wall or the spherical shell becomes less than a specified minimum value, or the radius 
of the storage volume exceeds a specified maximum. The minimum thicknesses should 
be specified consistent with fabrication constraints; whereas, the maximum storage 
radius relates to system integration considerations. When no storage volume is em- 
ployed, if die lube thickness required for containment is less than the specified mini- 
mum, the tube weight is calculated for the specified value. 

Also calculated in this subroutine are the temperature drops associated 
with conduction across the heat pipe wall at the evaporator and condenser sections. 
These temperature drops are a function of wall thickness and are therefore affected 
by the size of the storage volume. The equations for these temperature drops are 
presented in the next section. 

C. 1.2, \ Heat Transfer Analysis 

The heat transfer analysis Is based on the thermal model shown in Figure 
C. 1— I. The total thermal impedance R^ is composed of a series of individual resis- 
tances. 


R •» R t R + R + R 
w,o e v \v,e c 
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where the individual resistances are calculated from: 


D In ( -~\ 

0X1 V » lnt ) 


2 K A . 
\v ext 
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Vapor Space 


Wick 


Wall 



Figure C. 1-4. Thermal Impedance Model for Heat Transfer Analysis 
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R 

e 


h A 
e 


int,e 


P T 

R = ir^TT 
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R 


c 


h 

c 


JL _ 

A. ; 
mt,c 
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The internal areas used in the above equations are the actual heat transfer areas as- 
sociated with the particular wick geometry. Keat transfer coefficients for the slab, 
arterial, and grooved wicks are input for the evaporator and condenser sections. For 
the case of the annular or circumferential wicks, the film coefficients are calculated 
from: 


h 

an 




(R. - R ) 

1 V 
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h . 
cir 


eff 


(R - R ) 
' i v' 
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where the conductivities are specified input values. Individual temperature drops are 
calculated from: 


AT, = R. Q C. 1-45 

i i 

and the overall temperature drop is calculated as: 

A T = R u Q C. 1-46 

to 


As mentioned previously, the heat transport analysis will be terminated if the vapor 
temperature drop A T y exceeds a specified maximum value. 
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C. 1,2.5 . Analysis Summary 


The features of this analysis ana the basic assumptions employed can be 


summarized as follows: 


C. 1. 2. 5, 1 Features 


1. The heat transport capability (Q L ef f) or maximum heat transport 
<Q max > can be calculated for five basic wick geometries. 

2. The vapor flow may be laminar or turbulent. 


3. The heat transport capability is determined for both one "g" and 


zero "g' environments. 


4. The effect of elevation on heat transport is determined. 

5. Roth composite and homogeneous modes of operation can be anal- 
yzed for the slab, circumferential, and arterial wick geometries. 

6. The maximum heat transport is calculated as a function of wick 
thickness (or liquid flow area). For the axial groove geometry, 
tLe number of grooves, their width, and the corresponding land 
thickness is determined so that the maximum heat transport is 
realized for a particular aspect ratio. In this case, the maxi- 
mum heat transport is calculated as a function of groove aspect 
ratio. 

7. The minimum vapor area allowed without incurring the sonic 
vapor limit is calculated. 

8. The effective pumping radius required for self-priming is cal- 
culated. 

*9. The fluid inventory, evaporator and coudenser film temperature 

♦When the axial groove geometry is analyzed, the parameters of (5) and the weight of 
0) are calculated for those aspect ratios whose calculated heat transport exceeds the 
specified requirement. 

C. 1-17 




drops and the vapor temperature drop are calculated for the mini- 
mum wick thickness necessary to meet a specified heat transport 
requirement. 

*10. An optional weight analysis can be performed for the minimum wick 
thickness which satisfies the specified heat transport requirement. 
This analysis calculates the total system weight consistent with 
containment of the working fluid at a specified maximum tempera- 
ture. A spherical storage volume may be employed for contain- 
ment. a: the user f s option. 

11. The working fluid can be treated as a saturated vapor, a real gas, 
or an ideal gas (whichever is appropriate) at the maximum temper- 
ature. 

12. The heat transport analysis is terminated at either the sonic vapor 
limit or when the vapor temperature drop exceeds a specified value. 

C. 1. 2. 5. 2 Assumptions 

1. One-dimensional axial fluid flow, i. e. , radial and/or circumfer- 
ential pressure losses, are negligible. 

2. Heat pipe cross-section is circular and uniform over the length 
of the heat pipe. 

3. Wick cross-section is uniform over the length of the heat pipe. 

4. Fluid and material properties are constant over the length of the 
heat pipe. 

5. Uniform heat addition and removal. 

G. Axial heat conduction in the tube wall is negligible. 


*Wheu the axial groove geometry is analyzed, the parameters of (i) and the weight of 
(j) are calculated for those aspect ratios whose calculated heat transport exceeds the 
specified requirement. 
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( 7. Vapor pressure drop due to gravity is negligible. 

S. Momentum losses in the vapor are negligible. 

9. Liquid flow is always laminar. 

10. Nucleate boiling and entrainment limits are not considered. 

These assumptions are .consistent with most heat pipe analyses and are generally not 
prohibitive. The last three assumptions become invalid when very high heat loads 
or operation with a fluid at a low vapor pressure are required. These conditions are 
generally encountered with liquid metals at operating temperatures substantially below 
their boiling point, and care should be exercised when using this analysis under those 
circumstances. 

C. 1. 3. 0 Program Description 

C. 1.3. 1 General 

A listing of the program is presented in Appendix A. The program was 
written in FORTRAN V and was designed to operate on the UNIVAC 1108 system. The 
FORTRAN names and the physical quantities they represent are listed in Appendix B. 
Storage requirements are on the order of 50, 000 words (octal). 

The program logic is illustrated in the flow diagram contained in Appendix 
C. Basically, the program reads the input data, circulates equivalent wick properties, 
performs a heat transport and thermal analysis as a function of wick thickness, per- 
forms a weight analysis if required, and outputs the data. 

The deck setup as shown in Figure C. 1-5 consists of job control cards, 
the program source deck (w'hich may include a fluid property data acquisition code), 
additional control cards followed by the input data and program termination cards. 

As indicated above, the program has two major cptions. First, a fluid 
property data acquisition code can be utilized in conjunction with the IIPAD source 
program. This eliminates the need for inputting the various fluid properties at the 
specified vapor temperature as part of the data deck. The second option is the weight 
analysis. Either option is exercised by using the appropriate integer in the fourth 




input card of the data deck as described in the next section. 


C. 1. 3. 2 Input Description 

Table C. 1-2 describes the entries to be made on the various input cards 
and indicates when each of the optional cards are to be excluded. The FORTRAN 
name, format, and units to be used are indicated for each entry. A listing of sample 
input data for a homogeneous circumferential wick and for the axial groove wick is 
presented in Table C.l-3 and C.l-4. 

C. 1. 3. 3 Output Description 

The program outputs essentially all input data. This is followed by the 
heat transport analysis. The self-priming requirement and the minimum allowable 
vapor area (based on not exceeding the sonic limit) are printed. A table in which the 
film and vapor temperature drops and the calculated heat transport performance 
parameters are listed as a function of wick thickness or groove aspect ratio is printed 
next. This table is followed by a summary of the parameters associated with the min- 
imum wick thickness or aspect ratio’s which satisfy a specified heat transport require- 
ment. If the requirement cannot be satisfied by the particular design, the following 
statement will appear at the end of the preceding table: "No Area Exists to Satisfy 
QMAX Requirement”. If a weight analysis is requested, the tube thickness required 
for containment and a weight breakdown are listed in tabular form as a function of the 
wick storage volume radius for the minimum thickness or those aspect ratios which 
satisfy the heat transport requirement. This is then followed by a table which lists 
the conduction temperature drops and the syrtem’s temperature drop as a function 
of the storage volume radius. A listing of typical output data is presented in Appendix 
D with the sample problem, 

C. 1 . 3. 4 Nomenclature 

The nomenclature used in the Heat Pipe Analysis and Design Code section 
is given in the following tabulation. 


Table C.l-2 Input Data Description 


Input 

Card Fortran 


Format 

Name 

Descr ’>tion 

2A6 

HDl 

Headings (working fluid) 


HD2 


5A6 

HD3 

HD4 

Headings (heat pipe iterial) 


HD5 
HDl 2 


10A6 

HDG 

HD7 

HD8 

Headings (type of wick) 


HD9 
HD10 
HD11 
HDl 3 



HD14 
HDl 5 
HD16 


613 

MORE 

Control Point, integer <1 for last set of 
data, otherwise integer >1 


OPWT 

Control Point, integer < 1 without weight 
analysis . otherwise integer > 1 


DATA 

Control Point, integer > 1 no data acquisi 
tion code, otherwise integer < 1 


OTOVV 

Control Point, type of wick geometry 


1 slab 

2 arterial 

3 annular 

4 circumferential 

5 axial grooves 

STATE Control Point, integer >1 use Beattie- 

Bridgma;! Equation, integer < 1 use Ideal 
Gas Law 

FLUID Control Point, type of working fluid (DAC) 
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Unit 


o 




Card Fortran 


No. 

Format 

Name 

Description 

Unit 

5 

3F10.5 

TEMP 

Operating tempe rature 

°K 

- 


TMAX 

Maximum system temperature 

°K 



DT 


Maximum allowable vapor temperature 
drop 

°K 

6 

6E10.4 

RHOL 

Liquid density 

kg/m 3 



RHOV 

Vapor density 

kg/m 3 



XLAMD 

Latent heat of vaporization 

w-s/kg 



SIGMA 

Surface tension 

N/m 



XMUL 

Dynamic liquid viscosity 

kg/m-s 



XMUV 

Dynamic vapor viscosity 

kg/m-s 

7 

2F10.5 

XMW 


Molecular weight of working fluid 

kg/ mole 



GAMMA 

Ratio of the specific heats 


8* 

5E10.4 

ASO i 

1 

Constants for Beattie-Bridgman Equation 




SA | 

I 

(pressure in atmosphere, volume in liter/ 




BSO 


gm-mole, temperature in °K, R = 0.28206 




SB 


atm 1 ite rs/gm-mole K) 




SC J 




9 

E10.4 

PHI 


Wetting angle 

degrees 

10 

6F10.5 

QMAX 

Maximum heat transport 

W 



PERF 

Performance factor 

- 



FSAFE 

Safety factor 

- 



HIGH 

Elevation between the cor denser end and 

m 


evaporator end 


* Card 8 is needed only when the value of card 4 (5) is an integer smaller than 1 
and the value of card 4 (2) is an integer greater than 1. 



Input 


Card 

No. 

Format 

Fortran 

Name 

Description 

Unit 



SIGN 

Sign convention , +] positive elevation 
-I negative elevation 

- 



VOID 

Void fraction of the wick material 

- 

11 

3F10.5 

XLEV 

Length of the evaporator section 

m 



XL AD 

Length of the adiabatic section 

m 



XLCO 

Length of the condenser section 


12 

4EI0.4 

XOD 

Outside diameter of the heat pipe 

m 



TMIN 

Minimum wall thickness of the heat pipe 

m 



DELM 

Minimum wall thickness of the storage 
reservoir 

m 



RMAX 

Maximum allowable radius of the storage 
reservoir 

m 

13 

3E10.4 

RHOM 

Density of the heat pipe material 

. ' 3 

kg/m 



STRES 

Yield stress for the heat p*pe material 

N/m 2 



WALLK 

Thermal conductivity of the wall material 

W/m °K 

14* 

213 

I 1 

Mesh size of the coarse wick material 

- 



I 2 

Mesh size of the fine wick material 

- 

15 

5E10.4 

XKP 

Permeability 

2 

m 



RPE 

Effective pumping radius for heat transport 

m 



CRPE 

Effective pore radius for self-priming 

m 



HESL 

Evaporator film coefficient of slab wick 
heat pipe 

W/m 2 -°K 



HCSL 

Condenser film coefficient of slab wick 
heat pipe 

W/ m 2 -°K 


* Card 14 to Card 15 are needed only if the value of Card ^(4) is equal to one. 

\ 
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Input 

Card 


Fortran 



No. 

Format 

Name 

Description 

Unit 

16 

4E10.4 

D 1 

Diameter of the coarse wick material 

m 



D 2 

Diameter of the fine wick material 

m 



RHOW1 

Density of the coarse wick material 

kg/m 3 



RHOW2 

Density of the fine wick material 

kg/m 3 

17* 

13 

OTAR 

Control Card 





1 for closed artery 

2 for open artery 


18** 

5E10.4 

XMCWD 

Minimum wick thickness (i*e. , diameter 
of the artery) 

m 



STTNO 

Step increment for wick thickness 

m 



H 

Pedestal height of arterial wick 

m 



HEAR 

Evaporator film coefficient of arterial wick 

W/m 2 -°K 



HCAR 

•» 

Condenser film coefficient of arterial wick 

W/m 2 -°K 

19*** 

6E10.4 

XMCWD 

Minimum v/ick thickness (i.e. , diameter 
of the artery) 

m 



STINC 

Step increment for wick thickness 

m 



RPE 

Effective pumping radius for heat transport 

m 



H 

Pedestal height of arterial wick 

m 



HEAR 

Evaporation film coefficient of arterial wick 

W/ m 2 -°K 



HCAR 

Condenser film coefficient of arterial wick 

W/m 2 -°K 


* Cards 17 to Card 19 are needed only if the value of Card 4 (4) is equal to two (2). 
** Card 18 is not needed if the value of Card 17 is equal to 1. 

*** Card 19 Is not needed if the value of Card 17 is equal to two (2). 
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Input 

Card 

No. 

Format 

Fortran 

Name 

% 

Description 

Unit 

20 

2E10.4 

D 1 

Diameter of the wire of the screen mesh 

m 



RHOW1 

Density of the wick material 

kg/m 3 

21+ 

3E10.4 

XMCWD 

Minimum wick thickness (i.c., distance 
between the wick and the tube wail) 

m 



STINC 

Step increment for wick thickness 

m 



XKLIQ 

Thermal conductivity of the working fluid 

W/m-°K 

22* 

2E10.4 

D 1 

Diameter of the wire of the screen mesh 

m 



RHOW1 

Density of the wick material 

kg/m 3 

23** 

3E10.4 

RPE 

Effective pumping radius for heat transport 

m 



XKP 

Permeability 

o 

imi 

m 



XKEFF 

Effective thermal conductivity of the liquid 
and wick in circumferential wick heat pipe 

\V/m-°K 

24** 

2E10.4 

D 1 

Diaiieter of the wire of the screen mesh 

m 



RHOW1 

Density of the wick material 

kg/m 3 

25*** 

8E10.4 

ARMAX 

Maximum value of aspect ratio of the groove 

m 



ARMIN 

Minimum value of aspect ratio of the groove 

m 



WMAX 

Maximum value of the groove width 

m 



WMIN 

Minimum value of the groove width 

m 



TLIvlAX 

Maximum value of the land thickness of the 
groove 

m 



TLMIN 

Minimum value of the land thickness of the 
groove 

m 


* Card 21 and 22 are needed only if the value of Card 4 (4) is equal to three (3). 
** Cards 23 and 24 are needed only if the value of Card 4 (4) is equal to four (4) 
*** Cards 25 and 26 are not needed if the value of Card 4 (4) is equal to five (5). 
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I 


Input 

Card 

No. Format 

Fortran 

Name 

Description 

Unit 


HEGR 

Evaporator film coefficient of the grooved 
heat pipe 

W/m 2 -°K 


HCGR 

Condenser film coefficient of the grooved 
heat pipe 

W/m 2 -°K 

26*** 213 

NDAR 

Number of groove aspect ratio 



NDWD 

Number of groove widths 



*** Cards 25 and 26 are not needed if the value of Card 4 (4) is equal to five (5). 
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nitrogen 

ALUMINUM «LLOY 6’J61“T6 

HOMOGENEOUS WICK AGAINST THE WALL (200 MESH SCREEN) 
2 2 2 A 1 2 


b0. 

300. 

10*0 

U.OOOOE+OO 



^.0 

1.0 

2.0 

0.2 

0.8 

0.2 


I.27u0t-028.9000E-048.9000E-OAS.8000E-02 
2.7000c ♦032. 7600E*081.6000E*02 
7.7000E-116.8500E-051.4000E-01 
5.33o0t-052.850Cc>03 


Table C. 1-3. Sample Input Data for a Homogeneous Circumferential Wick 
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NITROGEN 

ALUMINUM ALLOY 6061-T6 
AXIAL RECTANGULAR GROOVED HEAT PIPE 
0 2 0 5 1 2 

80. 300. 10.0 

8.0000E»025.9000E»001.9500E*058.2000E-031.4480E-045.34oOE-08 
28.0 1.4 

c.ooooe*oo 

0.005 1.0 2.0 .00254 1.0 0.7 

0.2 O.o 0.2 

1.2700E-028.9000E-048.9000E-045.0800E-02 
2.7000E*032.7600E*081.6000E*02 

1.5000E*005.0000E-017.6200E-044.0600E-045.1000E-043.8100E-043.QOOOE+027.6000E*02 
3 3 


Table C. 1-4 Sample Input Data for a Axial Groove Wick 


NOMENCLATURE 


Symbol 


Description 


A 

A t Aq, B, Bq 
D 
Dh 
Fl 
H 
K 
L 
N 
P 

P w 

Q 

R 

Re 


T 

V 


Area 

Constants for Beattie-Bridgman Equation 

Tube diameter 

Hydraulic diameter 

Pressure drop ratio 

Wicking height factor 

Permeability 

Length 

Transport factor 

Pressure 

Wetted perimeter 

Axial heat flow rate 

Radius, thermal resistance 

Reynolds number 

Gas constant 

Safety factor 

Temperature 

Volume 


a, b, c 
f 

g 

h 

k 


m 



w 


Constants for Beattie-Bridgman Equation 

Friction factor 

Acceleration 

Heat transfer coefficient, elevation 

Thermal conductivity 

Mass 

Pumping radius 
Thickness 
Groove width 


a 

0 

y 

6 

c 

v 

6 

X 

u 

v 

f 

d 


Groove half angle 

Heat pipe orientation with respect to gravity 
Ratio of specific heats 

Groove depth, thickness of wall of storage volume 

Porosity 

Gravity factor 

Contact angle 

Heat of vaporization 

Dynamic viscosity 

Kinematic viscosity 

Density 

Surface tension 
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Subscripts 

Symbol 


Description 


n 

an 

c 

cir 

e 

eff 

ext 

hp 

i 

int 

1 

max 

min 

req 

sp 

st 

t 

V 

W 

w 


Adiabatic 
Annular wick 
Condenser 

Circumferential wick 

Evaporator 

Effective 

External 

Heat pipe 

Index, individual 

Internal 

Liquid 

Maximum 

Minimum 

Required 

Self-priming 

Storage volume 

Total 

Vapor 

Wick 

Wall 
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a te 


i 



\ 

Appendix A. Flow Diagram of 
Heat Pipe Analysis and Design Code 


ppFC.FnTNJC P h ^ 


VIOT pnAt^ 













calculation 



parameters 


THK Of HEAT PIPES 


PERFORM The THERMAL 
analysis FOR each 


WALL THICKNESS 
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Appendix B, FORTRAN Names 




Fortran 


Name 

Description 

Units 

HD1, HD2 

Headings (working flu id) 


HD3, HD4, 
HD5, HD12 

Headings (heat pipe material) 


HD6, HD7, 
HD8, HD9, 

Headings (type o.C wick) 


HD10, FD11, 
HD13, HD14, 

■ ' 


HD15, HD16 



MORE 

Control Point, integer = 0 for last set of data, 
otherwise inter ger = 2 


OPWT 

Control Point, integer < 1 without weight analysis, 
otherwise integer > 1 


DATA 

Control Point, integer < 1 no data acquisition code, 
otherwise integer >1 


OTOW 

Control Point, type of wick geometry 

1 slab 

2 arterial 

3 annular 

4 circumferential 

5 axial grooves 


STATE 

Control Point, integer >1 use Beattie -Bridgman 
Equation, integer < 1 use Ideal Gas Lav/ 


FLUID 

Control point, type of working fluid 



1 hydrogen 

2 nitrogen 

3 oxygen 

4 water 

5 ammonia 

6 methanol 

7 acetone 

8 freon- 21 

9 sodium 

10 potassium 



11 lithium 



12 mercury 

o 

TEMP 

Operating temperature 

K 

TMAX 

Maximum system temperature 

°K 

DT 

Maximum allowable vapor temperature drop 

°K 

RHOL 

Liquid density 

kg/m* 

RHOV 

Vapor density 

kg/ m s * . 

XLAMP 

Latent heat of vaporization 

W s/kg 

SIGMA 

Surface tension 

N/m 

XMUL 

Dynamic liquid viscosity 

kg/m s 

xmuv 

Dynamic vapor viscosity 

kg/m 3 
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Fortran 

Name 


XMW 

GAMMA 

ASO 

SA 

BSO 

SB 

SC 

PHI 

PIH 

QMAX 

PERF 

FSAFE 

HIGH 

SIGN 

VOID 

XLEV 

XL AD 

XL CO 

XOD 

TMIN 

DELM 

R MAX 

RHOM 

STRES 

WALLK 

II 

IND 

QTRAN 

XLHP 

XLEFF 

QPRED 

XI D 

AE 

AC 

RHEAD 

XNL 

11 

12 

XKP 

RPE 


Description Units 

Molecular weight of the working fluid kg/mole 

Ratio of the specific heats 

Constant for Beatfie-Bridgman Equation 


Wetting angle 
Wetting angle 
Maximum heat transport 
Performance factor 
Safety factor 

Elevation between the condenser end and evaporator 
end 

Sign convention, +1 positive elevation, -1 negative 
elevation 

Void fraction of the wick material 

Length of the evaporator section 

Length of the adiabatic section 

Length of the condenser section 

Outside diameter of the heat pipe 

Minimum wall thickness of the heat pipe 

Minimum wall thickness of the storage reservoir 

Maximum allowable radius of the storage reservoir 

Density of the heat pipe material 

Yield stress for the heat pipe material 

Thermal conductivity of the wall material 

Indicator 

Counter for performance parameters 
Maximum heat transport requirement 
Length of the heat pipe 
Effective heat pipe length 

Maximum heat transport requirement (including 
the performance factor) 

Inside diameter of the heat pipe 
Internal area of the evaporator 
Internal area of the condenser 
Required head of self-priming for slab wick, 
circumferential wick, and annular wick 
Liquid transport factor 
Mesh size of the coarse wick material 
Mesh size of the fine wick material 
Permeability 

Effective pumping radius for heat transport 


degrees 

radians 

W 


m 


m 

m 

m 

m 

m 

m 

m 

kg/m 3 
N/m 2 
W/m °K 


W 

m 

m 

W 

m 

m 2 

m 2 

m 

W/m 2 


m 

m 


t 
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Fortran 



Name 

Description 

Units 

CRPE 

Effective pore radius for self-priming* 

m 

HESL 

Evaporator film coefficient of slab wick heat pipe 

W/m*°K 

HCSL 

Condenser film coefficient of slab wick heat pipe 

W/nr *K 

D1 & D2 

Diameter of the wires of the screen mesh 

m 

RHOW1 

Density of the coarse wick material 

kg/hi 3 

RHOW2 

Density of the fine wick material 

kg/m 3 

XMCWD 

Minimum wick thickness 

m 

CWL' 

Wick thickness 

m 

CPCOR 

Capillary pumping pressure for self-priming 

N/m 2 

CCORE 

Self-priming head 

m 

CPGHD 

Capillary pumping pressure for heat transport 

N/m 2 

BODY 

Body force 

N/m 2 

OTAR 

Control card 

(3) for closed artery 
(2) for open artery 


STINC 

Step increment for wick thickness 

m 

H 

Pedestal height of arterial wick 

m 

HEAR 

Evaporator film coefficient of arterial wick 

W/rn? °K 

HOAR 

Condenser film coefficient of arterial wick 

W/m 2 °K 

DMAX 

Maximum allowable diameter of artery for se)f- 
priming 

m 

XKLIQ 

Thermal conductivity of the working fluid 

W/m °K 

XKEFF 

Effective thermal conductivity of the liquid and 
wick in circumferential wick heat pipe 

W/m °K 

SONIC 

Sonic velocity at operating temperature 

m/s 

A SONIC 

Minimum allowable vapor area for sonic limit 

m* 

NL 

Indicator 

- 

EHEAD 

Equivalent pore radius for self-priming 

i: 

B 

Half of the wick thickness 

U\ 

THETA 

Angle 

rad *ms 

C 

Chord length 

■ >.i 

AV 

Vapor area 

n* 

CHEAD 

Net capillary' pressure for heat transport 

N/m 2 

AW 

Wick area 

m 2 

WP 

Wetted perimeter of liquid area 

m 

DHV 

Hydraulic diameter of the vapor 

m 

ARC 

Arc length 

m 

HAE 

Evaporator film conductance 

W/K 

HAC 

Condenser film conductance 

W/K 

DHL 

Hydraulic diameter of the liquid 

m 

QLMAX 

Maximum heat transport factor in "l-g" 
environment 

W m 

QLMOG 

Maximum heat transport factor in ’’O-g*’ 
environment 

W m 
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Fortran 


Name 

Description 

Units 

REYND 

Reynolds number 

- 

DIFT 

Vapor temperature drop 

# K 

QMAXC 

Maximum heat transport in rt l~g” environment 

W 

QMAOG 

Maximum heat transport in environment 

W 

DQDH 

Slope of heat transport versus evaporator elevation 

W/m 

AETD 

Evaporator temperature drop 

*K 

ACTD 

Condenser temperature drop 


XQ(IND) 

Temporary storage space for QMAXC 

W 

XAftND) 

Temporary storage space for AV 

m* 

XB(IND> 

Temporary storage space for CWD 

m 

XL(IND) 

Temporary storage space for AW 

m* 

XQLflND) 

Temporary storage space for QLMAX 

W/m 

XARC{IND) 

Temporary storage space for ARC 

m 

XDIFT(IND) 

Temporary storage space for DIFT 

K 

AREA 

Vapor area 

m* 

AREAL 

Liquid area 

m* 

VTD 

Vapor temperature drop 

•k 

LAYER 

Number of layers ut screen mesh 

— 

XMFLD 

Mass of the working fluid 


VOW 

Actual vjlume of tl* wick 


WTW 

Weight of the wick 

kg 

ETD 

Evaporate" temperature drop 


CTD 

Condenser temperature drop 

°K 

ei j 

Constants for turbulent vapor flow 

— 

C2 > 
C3 ) 
CHECK 

Counter for iterations when vapor is turbunnt 


QT(I) 

Temporary storage space forQ max 

w 

A(N) 

Constant for turbulence calculation 

N/m* 

DA(N) 

Derivative of turbuient transport equation 

- 

QTEMP(II) 

Temporary storage space for Qirru: 

w 

ARMAX 

Maximum value of aspect rail ) 

. - 

ARMIN 

Minimum value of aspec ; * to 

- 

WMAX 

Maximum value of the groove, w^th 

m 

WMIN 

Minimum value of the groove /i-it h 

m 

TLMAX 

Maximum value of the land thickner^ 

m 

TLM1N 

Minimum value of the land thickiie&s 

m 

HEGR 

Evaporator film coefficient of the grooved heat pipe 

W/m*°K 

HCGR 

Condenser film coefficient of the grooved heat pipe 

W/m* *K 

NDAR 

Number of groove aspect ratios 

- 

NDWD 

Number of groove widths 

— 

DNAR 

Number of groove aspect ratios 

— 

A RING 

Increment for the aspect ratio 

— 

DAR(I) 

Aspect ratios 

— 

DAWSP(I) 

Maximum allowable width of the groove at an aspect 

m 


ratio for self-priming 
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Fort ran 






Name 

Description 

Units 

! 



DREQ 

Required pumping height 

m 




DNWD 

Number of groove widths 

- 




\V1NC 

Increment for the groove width 

m 




PA 

Groove width 

m 




XYZ(J) 

Groove width 

m 




DEP(J) 

Groove depth 

m 




SUV 

Internal radius of the grooves 

in 




NMAX 

Maximum number of grooves 

- 




NM3N 

Minimum number of grooves 

- 


i 


XK 

Number of grooves 

- 




DXL (K) 

Land thickness of the groove 

m 




BASE 

Base thickness of the groove 

m 




SAV 

Vapor area 

m* 




SAI 

Vapor area 

m* 




AGO 

Area of one groove 

m* 




AMO 

Meniscus area of one groove 

m* 




ALO 

Liquid area of one groove 

m a 




GVOWl(K) 

Total liquid area 

m e 




ALAND 

Total land area 

m* 




DD1(K) 

Maximum heat transport factor for groove 

W m 





geometry in "1-g" environment 


| 



EEl(K) 

Maximum heat transport for groove geometry 

w 





in "1-g" environment 





FF1(K) 

Maximum heat transport factor for groove 

W 111 





geometry in "0-g" environment 





GG1(K) 

Maximum heat transport in "0-g" environment 

w 




HH1(K) 

Slope of heat tr:ms|>ort versus evaporator elevation 

W 111 





in "l-g” environment 





DGM(K) 

Mass of the working duid 

kg 




D\VICK(K) 

Weight of the land of the groove 

kg 




NG1(K) 

Number of grooves 

- 




NG2(J) 

Number of grooves 

- 




BB2(J) 

Groove width 

m 




CC 2(J) 

Land thickness 

m 




DD2(.T) 

Maximum heat transport factor for groove geometry 

W m 





in "l-g" environment 





EE2(J) 

Maximum heat transport for groove geometry in 

w 





"l-g" environment 





FF2(J) 

Maximum heat transport factor for groove geometry 

W 111 





in "0-g" environment 




- 

GG2(J) 

Maximum heat transport for groove geometry in 

w 





"O-g" environment 





HH2(J) 

Slope of heat transport versus evaporator elevation 

W ill 

I 




in "l-g" environment 


1 . 



P P2(J) 

Mass of the working fluid 

kg 
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Fortran 

Name 


Description 


Weight of the land of the groove 
Vapor temperature drop 
Total liquid area 
Land thickness 
Number of grooves 

Maximum heat transport factor for groove geometry 
in ”l-g M environment 

Maximum heat transport for groove geometry in 
"l-g M environment 

Maximum heat transport factor for groove 
geometry in "O-g" environment 
Maximum heat transport for groove geometry in 
•*0-g n environment 

Slope of heat transport versus evaporator 
elevation in "1-g" environment 
Vapor temperature drop 
Total liquid area 
Temporary storage space 
Temporary storage space 
Case number 

Wall thickness of the storage reservoir 
Wall thickness of the heat pipe 
Gas constant 

Wick area 

Internal radius of storage reservoir 

Weight of storage reservoir 

Required wall thickness of storage reservoir 

Minimum wall thickness of storage reservoir 

Density of the heat pipe material 

Density of the working fluid 

Internal radius of the heat pipe 

Internal diameter of the heat pipe 

Length of the heat pipe 

Maximum allowable radius of storage reservoir 

Outside diameter of the heat pipe 

Maximum system temperature 

Minimum wall thickness 

Mass of the working fluid 

Specific volume 

Coefficients of Beattie, Bridgman Equation 

Internal pressure of the system 
Internal pressure of the system 


atm 1 

gmole°K 

ft* 


m 

lb/ft 5 

lb/ft 5 
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Fortran 

Name Description Units 


BTRED 

TRED 

BVHP 

VVTL 

WTH 

WT. 

BINCR 

BVST 

BVTT 

DEXT 

DINT 

AEXTE 

AEXTC 

EWTD 

CWTD 

TTD 


Calculated wall thickness of the heat pipe 

Calculated wall thickness of the storage reservoir 

Volume of the heat pipe 

Weight of the working fluid 

Weight of the heat pipe 

Total weight 

Increment for the radius of the storage reservoir 

Volume of the storage reservoir 

Total volume of the system 

Outside diameter of the heat pipe 

Inside diameter of the heat pipe 

External area of evaporator 

External area of condenser 

Temperature drop across the evaporator wall 

Temperature drop across the condenser wail 

Total temperature drop of the system 


ft 

m 

ft* 

kg 

kg 

kg 

ft 

ft* 

ft* 

m 

m 

ni* 

in* 

*K 

*K 

*K 
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dimension xci7ooi,xai20oi,x8(20oi,xli2ooi * xqluooi .Qtm 

DIMENSION 41 *>,DA( M*UTFMP!5I. ADIFT! 200>,AARe!200) 

DIMENSION OAR! 2001 .DAWSP1200)*0Ri-C<20O),xy2t 2 001.UEPC7001,L>XL(200I 
01 HE NS I UN 001(2 00 t*EEl(20U)*FF 1 1 2001 , GG U 200 I ,HHI ( 200) « DGM ( ?(JU I 
01 MENS ION OGW( 2001 «DW I CM 200) *NG i ( 200 1 *NG2 1 2 001 *082! 2UQ I *CC 2! 700) 
01 MENS ION DD7|?00)»EE2t2Q01*FF2! 2G01* GG 212001* HH 71? 00 )*PP 2(200) 

01 MENS ION 0Q/< 200 ),NG( 2001*001 I 200 1 • 0G2( 200 I ,0.01(2 001 ,l)G4(?001 
01 MENS ION DG4(2GO),NCASEt?0O),DV!O(2CCl,RK2(2OQl 
01 MENS I CN GV'jW!l200)*GVO*?l2n0),GVUW)(200) 

INTEGER OPWT,OTGW*GTAR,DMA,SIArE,fLUlU 
5005 FORMAT U0A61 
498 FORMAT (bill 
5010 FORMAT (8F10.5) 

5009 FORMAT (8E10.4) 

61 FORMAT (//// 24H INPUT 0 A T A (* 10A6,LH)I 

63 FORMAT I//19H TUBE MATERIAL *5A6) 

65 FORMAT ( /2E 20. 7 * 10X* F4 .2 1 
TO FORMAT I/4E20.7) 

176 FORMAT (//2X,10A6) 

607 FORMAT ( /E20. 7, 15X*E20.7* 24X.E20. 7 1 

609 FORMAT I//23H MIRE DIAMETER OF THE ,I3,39H MfcSH SCREEN WIRE 

10IAMETER OF THE ,I3*I5H MESH SCREEN(M)1 

610 FORMAT l/E20.7 t 21X*E20. 71 

611 FORMAT 1//17H DENSITY OF THE , !3*4QH MESH SCREEN! KG/M-K-M 1 DENS 
l ITY OF THE * 1-3* 22H MESH SCREEN! KG/M-M-HI I 

76 FORMAT (//I8H WORKING FLUID - ,?A6 *4H AT *F9.2*3H(KM 

622 FORMAT </4X*F6.1) 

623 FORMAT i/4X,FB.3) 

891 FORMAT !/4X»F10.2) 

86 FORMAT (Z///76H SELF-PRIMING REQUIREMENT 
l UN 1G WJ T H ELEVATION 0F.E20. 7.4HIM) I 1 
634 FORMAT ( /620.7 * 20X , E20.7, I 3X *E 20. 7 1 
637 FORMAT (/£20.7*69X*£20.7> 

202 FORMAT (/E20-7.10X.E20.7: 

50 FORMAT (////51H HEAT TRANSPORT ANALYSIS! 

1 • lOAb* iHl I 
10 FORMAT ( * l • ) 

60 FORMAT (/////• HEAT PIPE ANALYSIS AND 
ID E S I G N IHPAD)*) 

62 FORMAT (////• HEAT PIPE PROPERTIES*) 

720 FORMAT (//* THERMAL CONDUCT IV I T Y( W/M-K )• 1 

64 FORMAT .(//* YIELD STRE$S(N/M-MI DENSITV(KG/M-N-Hl CONTAIN 
ItR SAFETY FACTOR* I 

66 FORMAT I//* OUTSIDE 01 AMETEM Ml KIN- WALL THICKMESSIMI • 1 
66 FORMAT !//' LENGTH OF THE HEAT PIPE SECT IONS* 1 

69 FORMAT (//• EVAP(N) AOJM) COHO 

I IHl TOTAL (HIM 

400 FORMAT !////• GEOMETRY OF THE STORAGE 

IVOIUU'I 

71 FORMAT (//* RAOIUSI M) OF MAX. STORAGE VOLUME 1 1 
73 FORMAT (////• WICK GEOMETRY*) 

605 FORMAT (//* Y ; JiO FRACTION OF THE WICK*) 

606 FORMAT (//• p:RME4BILITY(M-M| FFF. PUMPING RAO. FOR HEA 

IT TRANSPORT!*! EFf. PUMPING RAO. FOk SELF-PR IMZNG(M) • 1 

608 FQkMAT (//• MIN. THICKNESS OF THE WiCKlMI*) 

614 FORMAT (//• MIN. DIAMETER OF THE AKFFRYIH1M 

615 FORMAT I//* WIRE DIAMETER OF THE SCktEN(M)') 

617 FORMAT 1/7* DENSITY OF THE SCREEN W IRE I KG/M-M-M) • I 


618 FORMAT I//* FFF. PUMPING RAO. FOR HEAT TRANSPORT (M) • I 

619 format (//• MIN. OISTANCE BETWEtN THE WALL AND THE WICK(M)*) 

620 FORMAT !//• PfcRMEARlLITYIM-Ml EFF. PUMPING RAD. FOR HEAT TRA 

JNSPORT ANO SELF-PRIM1NG(MI* | 

rr format (//• fluid properties* i 

621 FURMAT (//• MOLECULAR WEIGHT! KG/MOL E ) * I 

4622 FORMAT (//• THE RATIO OF SPECIFIC HEATS! CP/CV I • I 
72 B FORMAT (//* THERMAL CONDUCT I V I TY( W/M-K I • I 

7)0 FORMAT (//• EFFECTIVE THFRMAL CONDUCTIVITY OF THE LIQUID AND THE 
1WICK(W/M-K)*> 

624 FURMAT (//• CONTACT ANGLE (DEG) * I 

79 FORMAT (//• SURFACE TENS ION(N/MI LATENT HEAT ( J/KG) * I 

81 FORMAT (//• L IQ. DENSITYIKG/M-H-MI VAPOR DENSITY I KG/M-M-HI • | 

6 i FORMAT (//* L 10. VISCOSITYIKG/M-SI VAPOR VI SCOSITY (KG/M-S I • I 

104 FQKMAT (//• LIQUID TRANSPORT FACTOR (W/M-M ( •) 

888 FORMAT (////• PERFORMANCE REQUIREMENT 5* I 

889 FORMAT I//* ELEVAT IONIH t * ) 

4202 FORMAT (//• PEDESTAL HElGHTfNI 1 I 

890 FORMAT (//• MAX. TEMPERATURE! K) * I 

143 FORMAT (//• MAX. ALLOWABLE VAPOR TEMPERATURE OROPIKIM 
85 FORMAT (////• OUTPUT D A* T A * 1 
633 FORMAT !//* MAX. CAPILLARY HEAD FOR SELF-PRIHINGI Ml KE0U1RED 

I HEAD(H) EQUIVALENT EFF. RAD. FOR SELF-PRIHINGINIM 

635 FORMAT (//• MAX. ALLOWABLE DIAMETER OF THE ARTERY FOR SELF-PRIMIN 
IGIH1M 

636 FORMAT (//• MAX. ALLOWABLE DISTANCE BETWEEN THE WALL AND THE WICK 

1 FOR REQUIRED HEAD! M 1 REQUIREO HEADCN) • I 

200 FORMAT !//• SONIC LIMITATION 1 ) 

201 FORMAT !//• SONIC VELOCITY(M/SI MINIMUM VAPOR AREAIM-M 

1 I FOK SONIC LIMITATION* 1 

4630 FORMAT !//• THICKNESS OF VAPOR TEMP EVAPORATOR CON 

10ENSER <G«It IG-ll !C-0) !G>0) 

2 (C*1IM 

4631 FORMAT (//• DIAMETER OF VAPOR TEMP EVAPORATOR CON 

lOENSER (G-il (G-ll CG-OI 1G-0I 

2 (G»l )• ) 

4632 FORMAT (//• DISTANCE BETWEEN THE VAPOR TEMP EVAPORATOR CON 

1DENSER (G-l) !G*lt (G-0) !G-0) 

2 (G-l)*) 

4633 FORMAT !* WALL ANO THE WICK DROP TEMP DROP TEM 

IP DROP (QLEFF)MAX QMAK IQLFFFIMAX UMAX 

2 (OQ/DH) * ) 

U FORMAT (/* IN) IK| (K) 

1!X) (N-M) IWI IW-MI | W) 

2 I W?M| ■ 1 

3630 FjRMAT I* THE WICK DROP TERM DROP TEN 

IP DROP (QLEFF)MAX QMAX (QLEFF 1MAX QMAX 

2 (DQ/DH)*) 

23 FORMAT (//• NO AREA EXITS TO SATISFY XMAXIW) REQUIRMENT 1 ) 

31 FORMAT !//• OMAXIW) REQUIREMENT I INCLUDING THE PERFORMANCE FACTOR 
1) IS*) 

126 FORMAT !//* THE PERFORMANCE FACTOR IS*) 

26 FORMAT (//» THE MAX. VAPOR AREA! M~N) TO SATISFY QMAXIW) REOUIREMC 
1NT ISM 

424 FURMAT !//• THE MIN. THICKNESS(H) OF THE WICK TO SATISFY QKAXlWl 
IRE QUIRE ME NT IS*) 

425 FORMAT !//• THE MIN. 01 AMETER(H) OF THE WICK TO SATISFY QMAXIW) R 
1EUUI REMENT IS*) 




42s format (//• r me if* x. distance (mi bei-een The wall an o inf wic* rrj 
l SATISFY QMAXlM J REOUIMENENT ism 
HOI FUAMA! I/EJ 6. J f >X,«E14.J) 

132 FORMAf (/7X,F4.2I 

555 FORMAT t/E20.7«lH(*(3*l?H L A Yf ft ( S I OF , 13, 1H- , I ) , 1 4H MESH SCRUNI I 
rro FCRMaT (/E 20.7,9X,E20.7,9;!,E?0.7) 

#04 FORMAT I/H6.3,3X. JF l4.3,8X«*MN.A. f 2X.M4. 3, *X *4HN. A. ,2X , E 14. J*9X, 

1 4MN a A • | 

722 FORMAT l/5X,E20.7,16X,t?0.7) 

V* F OAF* AT l/F?0.2,24X,E20.T) 

99 FORMAT C/4X, l J, 5X, MyJ* ,?6.2, 7X, UO. J, U.EIQ. J,2X,F 10. l,4X,fc 10.1,4 
lX,UG.3,4X*E10.3 v 5X,fVt.3l 

91 FORMAT I//* THE MIN. LIQUID ARfc At M-MJ TO SA T T SF V QMAXfWf PEOLlPFM 
IENT ISM 

93 FORMAT (//• MASS OF THE NORA I MG FLUID RE OUl AM) I *GI M 
145 FORMAT !//• VAPOR TEMPERATURE DROPlY.) EVAPORATOR TfcMP. DRDPIM 
1 CONDENSE A TEMP. OAORiriM 

151 FORMA T I // * FLOW IS NOW TUR&Ul ENT (REYNOLD NO. GREATER THAN 2 

1200 I a ' I 

203 FORMAT |//« MINIMUM MIC* AREA VIOLATES ThE SONIC LIMITATION* NO 
1H6AT IRANSROAT AND WEIGHT ANALYSES ARE PRESENTED • » 

13 FORMAT 1 / / • 7 ME MAI . AND *lN. VALUES OF ASHCT RATIO OF THE GROQV 

L t ( CROON E DERTh/OROOYE wIDTHPI 

14 FORMAT t//* THE MAX. ANO M IN. VALUES OF THE GROW NlCTH(M»M 

15 FORMAT (//* THE MAX. AND MIN. VALUES OF Tn£ L*M0 TMICFNESSIM) ' I 

721 format t//* evaporator film coefficient! w/m-m-xi condense* f 

IRM COFFFlCifNflw/M-M-M* I 

16 f t/RMA T I / / • THE ASPECT RATIOS OF THE GROOVE MAX. ALLOWABLE G 

1 R 00 V E WiOTH FUR S 1 1 F -PR 1 M INOI M 1 *» 

IT FORMAT (* (GROOVE OEPTH/GRCQvC WIOThIM 

35 FORMAT (//» 

1 10*11 10*11 (6*01 (0-01 
2I0-1IM 

36 FORMAT !/» CASE GROOVE ASPECT RATIO GROOVE LAND 

I I0LEFF1MAX UMAX tOlf'F)MAX OMAX 0 

2Q/W* I 

|7 FORMAT I* NO- NC. I OERTH/wIOTHl WlOTHtMl TMKVNES 

ISfMl 1W-M| ( «V | (M-MJ <«) ( 

?W/M1»I 

54 FORMAT (//• THE CASES WHICH SATISFY THE 0*AX<W» REQUIRE ME NT ARC L 

I I ST AS FOLLOWING* 1 

55 Forma I (/• 

l it* ii* » 

56 FORMAT I/* CASE GROOVE ASRICf RATIO G/CWVfc LANC 

i QMAX MASS OF TMF VAROR TEMP, EVAPORATOR CONOfcN 

2SER * 1 

57 FOwMAf (• 1*0- NO. tOtRTH/wlUTH) MlLtHlM) fH|CRNfSS 

!■) (W1 FIUILIXOJ D*OF»i.l TEMP DROP!*! I‘.M« 0 

<ROP(riM 

771 FORMAT (//• EVAPORATOR T£MR. ORfJPtX ) CJN.vFNSM Tfwp. UROPHO*) 
59 FORMAT (/4X.I3.5X.I3, 7X.F6.2.3X, 7U4.3I 
500 wt AD (5,59051 H0UHD2 
♦CUM*U 

■Rile 16*101 

wf AO 1 5* 5005 1 HO 3 , HD4 ,h 05 * h U 7 

w( AG (5*50951 hC 4 »hG 7 »moh,mC 9 ,hOIC*m 1 .U ,nOl J,moi a ,hD 1 5, HO 1 6 
-fAO (5,4991 MORE »OPwf , DA T A.GTu*. 5 T A f t , F LUIO 
h f AD (5*50101 Tf"P,TwAX,Of 


IF (OATA.Gr.il GO TO 5011 

wC AO (5.50091 RHOL , RHOV, XL AMO, SIGMA *XMUt#XMUV 
READ (5,50101 XMw, GAMMA 

GO TO 5023 
50 U LGUPP* 2 

CALL OAC I FLUID, L SO PP.LOtJT *RMOL ,RHUV, XL A"D, SI GMA , XMUL , XNUV, TfcMP.P 

1 V,XMw, GA hMA) 

IF U00I.9T.1I GO TO l 
5923 IF (STATE. EG. U GO ID 5022 . 

MEAD (5,5009 ) ASO,SA,BSO*Sfc#St 
5022 MfAO (5,50091 PHI 

PMI*PM1*3. 14159/180.0 
P I K* Ph (*180.0/3.14159 

READ (5*50101 OMAX.PEKF.FSAFE.Hl&f.STGN.VOID 
MEAD (5.50101 XLEV.XLAD.XLCQ 
READ (5,5009) XOD, .MlN.DELM.RMAX 
READ (5,50091 RHOM, STRE S , wM.U 
11*1 
1N0*0 

CTRAN-OMAX 

XI.HP-UEV*XlAD*XLCO 

XltfF-XLAD4XLtV/2.0*XLCO/'2.0 

QPKED«GMAX*PF*F 

O.GRREG 

xrD*XOD~7.0RTMIN 

AF-XID*J.14I59*XLEV 

AC-X10*3.14159*XtCO 

»HEAD-XID*SI9N»HIGH 

XNL- S I 9MARXL AMD4RMOL / XMUL 

R-XIG/2.9 

HS*MIGH*S1GN 

GO TO ( 330*331* 332,333* 12I*0TDW 

330 REAC (5,49# ) (1,12 

MEAD 15,50091 XXP.RPE * CRRE.HESL .HCSL 
READ 15,50091 01 .0? ,RM0W1 , RH0»2 
IF 111. £0. 121 GO TO 9000 
XMCwD* 2.0* ( 01*2.0*02 1 
GO TO 9902 
9005 XMtWD-2.0*01 
9002 twn«XNCWD 

CCC0R-2.0*SIGMA*C0SIRH11/CRRE 
CCORE'CPC0P./(RK)L*9.fcl 
CRGHO*2.0*SIGMA*COS(PHI t/APE 
8 OOY • RHOL *9 . S • ( X I D4 5 1 GN*H I GH I 
CHEAO*CRGHO-8GOY 
GO TO 334 

331 k«AD ( 5*493) OTAR 

GO TO (1332, 1331) *‘jT A* 
mi AC AO (5*50091 xmCWO.SUNC.H.HEAR.hCAR 
GO T,; 1331 

133? RF AD (5,50991 XMC WO, S T INC , RPE.H.HEAR.HCAR 
1333 ACAD 15,50991 DI.PhOMI 
CWO-XMCWD 

Ah»h*SIGN*HIGh 

AH*AHM6.0*SIGMA«CnS(RHl t/(RHGL«9.B))«*0.V-AH)/2.0 
CCuPC^DMAX*AH 
U'J lo 3 34 

31/ Pf ADI-MV0C9I XMCwOtSUNC.XrUO 


J 



Hf A0(5,5U09) DlttfHQMl 
OD-mCwD 

UMAx- 2.0 *SIGhA*COS(PH| )/|RHr;L*9.8MxlO»HlGH*:,lGNI ) 

CtoRt»GMM 

RHtAO“XIO*HIGH*5IGW 

GO TO 336 

»3i RfcAO<5,suo9i *kp,r?e,xm ff 
REA0I5, 60091 Dl.RHOWl 
xmGwd*?.o*di 
Cwu-XPCwO 

CPvH0»2,0*SIGma*CGSIPh|I/RPE 

bGOT*RHUl.*9.?*<XlO»SlGM*HtGH» 

Ch c *U-CPGHD-«ODv 

CCOKt-CPGHO/IRHOl*9.0» 

P.Ht AQ«X I D»Hl GH*S I GH 
D2* 0,0 

336 !.OHIOUGAMMA*I£wPwl, 8*1565. 3J*32.2/XMw)**0.5)*0. 3068 

AS0MIC»P£RF*GMAX*<2.G*(GANMA*1.0n**O.5/(XLAMD*RHGV*SOMlCI 

M.-1 

•*tT£ (6,101 
WRITE ( 6# 60 I 

want It. 61) HD6,H07,HD3,HDV,HOi0.HDll,HD13,HD16,HC15,HDlt 
WRITE (6.621 

WRITE (6.631 HD3.MD6.HD5.HC12 
WMU (6.720) 

WPirt (6,70) WALL* . 

WRITE 16,66) 

WRITE (6.65) SntES.RHOH, ESAFE 
WRIlfc (6.66) 
iMU (6,70) XOO.TMH 
WRITE 16.60) 

WRITE (6.69) 

WRITE (6.70) IL£V,XLAD,XIC0,XLHP 
WRITE 16,600) 

WRITE 16,71) 

WRITE (6.70) R WAX 
WRITE (6.10) 

WRITE (6,73) 

WRITE (6,176) HD6. H07.hOB.h09.HJ10.HC1 1 .HO 1 3, MO 16. HO I 5, HD 16 
WRITE 16.605) 

WRITE (6,623) VOID 

GO TO (601,602.603*606 | T OTOW 

601 WRIT* 16,606) 

WRITE 16,607) XXP.PRE.CRPE 
WRIIc 16,600) 

WRITE (6,70) XHCwD 
«Pil£ 16,60V) U.I2 
WRITE (6,6(01 01,02 
WRITE (6,611) 11,12 
WRITE (6,610) P.HCil »RHOW? 

GO 70 705 

602 GO TU (613, 612), OTAR 
612 W* 1 Tt (6,616) 

WRITE (6,70) XNC.tr!.’ 

WRITE ( 6,61 6l 
WRITE (6,70) (;t 
wRITt (6.617) 

WRITE (6, TO) wHOwl 




VO To 7V6 
613 WUfTtr (6,6181 

WRITE 16,70) MPE 
WHITE (6,616) 

WRITE (6,70) XMCWO 
W«|Tt 16,615) 

WRITE (6,70) 01 
WHITE (6,617) 

WHITE 16,70) RHOWl 
GO Tu 706 
603 WHITl (6,619) 

WHITE (6,701 XMCwD 
-RITE (6,616) 

WRITE (6,70) 01 
WRITE (6,617) 

WHITE It, 701 KhOWl 
GO TO 705 
606 WuITt (6.62Q) 

WRITE (6,607) XRP.RPE 
WRITE (6,615) 

WRITE (6,70) 01 
WRITE It , 61 M 
WRITE (6.70) RHOWI 
705 WRITE 16, 101 
WRITE (6,77) 

WRITE (6,78) HD1 ,HD2, TEMP 
WRITE 16,6211 
WRITE (6,622) XMW 
WRITE (6,662?) 

WRITE (6,623) GAMMA 
GO TO (723,327.7?5»72M,OTQW 
723 WRITE (6,721) 

WRITE (6,722) HESl.HCSL 
GO TO 727 
327 WRITE (6,721) 

WR1IL (6,72?) HEAk.HCAR 
GO TU 727 

725 WRITE (6,728) 

WRITE (6,70) XKLIO 
GO TO 72? 

726 WRITE 16,730) 

WRITE (6,70) XrfcFF 

727 WRITE U,6?6) 

WHITE (6,673) PIH 
WRITE (6,79) 

WHITE (6,70) SIGMA, XlftHr 
WRITE (6,01) 
wfilTl (6,70) RHOt.RHOV 
WRITE (6,83) 

WRITE (6,70) XHUUXMUV 
WRITE (6,(86) 

WHITE (6,70) XM4. 

WRITE < 6,10 
WRITE (6,8081 
WRIT! (6,311 
WRIT, (6,70) GPREO 
■RITE (6,126) 

WRITE (6,11?) PERF 



C.I-49 


WRITE (6,6891 
WHITE (6,01 US 

GO TJ 14200, 4201, 4200, 4200), gTOW 
4201 »»ITfc 16,42021 
WHITE (6,701 H 
4200 Wft|U 16,6101 

WRITE (6,091) mx 
j(K|Tc (6,143) 

write (6,691) or 

701 WRITE 16,101 
■KITE (6,8)1 
wRITfc (6,06) MS 
GO TO (630»63l»632»630),QTQW 

630 WRITfc (6,633) 

EMEAtW.0*SlGWA/(RH0L*9.8*RHEAD> 

WRITfc (6,634) CCQRfc,RM£AD,E.HEAO 
GO TO 613 

631 WRITE (6,636) 

WRITE (6,70) DMA* 

GO TO M3 

632 WRITE (6,636) 

WRITE 16,637) CCQRE.RHEAO 
313 GO TO 1305, 306, 307, 308), OTOW 
103 &«CW0/2.0 

THETA-ACOS(B/Rt 
C“R*Slh( THE T A| 

AV»2.0*A*R*THfTA-2.0*8*C 
GO TO 74 

306 AV"3.14159*(XID**?-X*CWD**2)/4.0 
GO TO 74 

30 7 AVO. 14139* (X|D-2.0*XHCWD 1 ••2/4.0 
GO TO 74 

306 AV-1. 14159* I XI0~2.0*XWCwD) **2/4.0 
74 WRITE (6,200) 

WRITE (6,201) 

WRITE (6,2021 SONIC, ASONtC 
IT (ASOMC.GT.AVI GO TO 3 

WRITE 16,50) H06,H07,HD8*H09,H010,H011 ,HD1 1, HD 14, HO 15, HOI 6 
GO TO (423, 626,62 T, 626 1, OTOW 
623 WRIT! 16,4430) 

WRITE (6,3630) 

GO TO 629 

626 WRITE (6,4631) 

WRITE 16,1630) 

GO TU 029 

627 WRITE (6,4632) 

WRITE (6,4633) 

629 WRITfc (6,11) 

5 GO TO (320, 321, 322, 323), OTOW 
320 0-CWU/2.0 

CH’f AO-CRGHD-^OGr 


wm«iNO«i 

IHE ( A* ACOS ( 8/ H ) 

C»R*SIN(THETA| 

AV*2.0*R***TMfcrA-2.0*B*C 

iw»3.l4l69*R*R-AV 

-P-4,0«(C**»*Th) (A) 

GHtf>4. J*A*/WI» 


KV»0HV/2-0 

ARC* 2,0* THETA/ 3 ,14 159 
HAt»HESl*AKC*AP 

hac«hcsc*arc*ac 

GO TG 324 

321 AV3. 14159*1 XlD**2-CW0**2l/4.0 
WP"3.14169*(XU*CW0I 
0HV»4 ■ 0*6 V/WP 
RY-DHY/2.0 

AW»3.L4159*C*D**2/4.0 

QHL-CWD 

XKP«DHl*0Hl/32.0 
ARC* 1 *0 
HAE*HEAR#AE 
hAC*HCAk*AC 

IF (CwO.GT.CmAX) GO TO 762 
GO TO 1351,332), OTAR 
351 CPGH0*2.0*SIG«A*C0S(PHI)/«PE 
GO TO 353 

362 CPGHD-2.0*SIGWA*COS(PH1)/ICWO/2,OI 

363 dOOr-RHOL*9.0*(AH*CWDI 
tHeAo-cPGHo-eoor 

I NO* 1N0* 1 
GO TU 124 

122 AV-3.14159*(XIO-2,0*CWO|**2/4.O 
WP»3.\4159*(X10-7,0*CWD> 

DHV-4,0*AV/WP 

RV-0HV/2.0 

0HL-2,0*CW0 

XXP-UHL*OHL/48.0 

AW»3.14159*XIO*X10/4.0-AV 

CP0H0«2.0*SIGMA*C0S(PHM/CW0 

B0DY-RH0L*9.a*(X!0*S!GN*HlGH| 

CHEAO-CPGHO-bOOr 
(NO* |NCt*l 
ARC" l • 0 

HAC«XKLIO*AE/CWO 
HAC"XKUQ*«C/CWD 
IF (CwQ.GT.ONAK) GO TO T62 
GO TO 324 

32 3 Av»l. 14169* (X ID-2. 0*CW0l **2/4.0 
*P-3.14159*<XID-2.0*CW0I 
0HV*4.0* AW/WF 
KV-DHV/2.0 

&W-3.l4169*Xt0*Xl0/4.0-AV 

CHEAC-CPGHO-BODY 

IhDMhD+l 

ARC-l.G 

HAt"XKffF*AE/CWO 
HAC* X* tF fc * AC/CWG 

324 wLHAX>0.6*CHt A0/(XHUL/<2.0*RHOL*AW*XKP*XLAHO)*4-0*XHUV/(AV*RHOV*XL 
l AHO*R V**2 • 0 1 1 

GLNUGsO.6*CPGHC/(XMUl/(2.04RHUL*AW*XRP*XLAWD)*4.0*XmjV/( AV*RHOF*XL * 
IAWG*SW**2.0)) 

R» YNO*?.0*(0LNAX/Xl EFF|*RV/|AV*XWU9*XLAR0) 

1 F UEYhG.tr. 2200. 01 GO TO 40 

GIFT - t8.0*7W , JY*GLfAZ/tAV*RHOw*XlAMD*RV**2n*TEMP/(XLA*'0*«H0V) 

- W AXC- CL MAX/Xl IFF 




OMAOg-OIMUG/XUFF 

0Q0M*-O.5*l9*b*«MCJl*SIGN/XUFF I / ( XNUl/ ( 2. U**HOL •AW*X*W*XL A HD) 
1XHUV/IAV*RH0V*XIAMD*WV**2.01I 
OUOH»COOH* (Ql MAX/ ABS(GLHAX) ) 

aeto-gmaxc/mae 

ACTO-QMAXC/HAC 

WRITE (6*0011 CtoD*D!F| v AETO»MCrD t ULMAX*GMAxC*gLMUG*QMAOG*OODH 
62 XQ ( I HD ) «OHAXC 
XA(I NOl-AV 

xeuscMCwo 

Xl(lNDl-AW 

xuluno»*olhax 
XAR cl 1ND)* AAC 
XOIFT ( !NDI*OIFT 
(F (OIFT.GT.OTI GO TO 1000 
GU TO (360*366* 366 * 362 ) *OTOW 
360 IF III. tO. 121 GO TO 9001 
STlNC»2.0*IOt«D2) 

GO TO 366 
9001 $T|NC-2.0*D1 
GO TO 366 
362 $T|MO2.0*Dl 
366 CWD«C wD+ST I NC 

GO TO (3TI,572.573*3T6)*0T0W 
?71 if (CWD.GT.XI01 CO TO 1000 
CO TO 316 
572 TCWD-CWD»H 

IF UCMD.GT.XJD> GO TO 1000 

IF ICUO.CT.OMAX) GO TO 1000 

GO TO 310 
373 rCMD-2.0*CMD 

IF irCVD.GT.XfO) GO TO 1000 

IF (CWD.GT.D»*X| GO TO 1000 

GO TO 310 
576 ICMD-2.0*CW0 

IF ITCmD.GT.XIDI GO TO 1000 
316 IF (AV.GE.ASONIC) GO TO 5 
1000 I MAX* 1 

IF UhO.fO.l) GO TO 762 
l’r ( X0( i .* .GE .CPAE 0 1 GO TO 6 
DO 21 J* 1*1 NO 
IF (XGUJ-QMEO) 21*22*26 
22 I MAX* J 

AA£A*XA(|MAX) 

CMD-X04 IMAX > 

AAEAl-XLt IMAX) 

ARC-XAACUMAXI 
V?D-XO|FT(|MAXl 
GO 10 25 
26 | WAX* J 
|M«|MAX~t 

AAEA»XA(|M>MXAI JMAXJ-XAI1MI 1*1 QPRED-XOI !*))/( XO(tMAX>-XQ<|M I) 

CMO *XB( |M|«(XBI |HAX>-XB(tH)!*(OPXEO-XOl IM)|/(XQ( I MAX l-XQI I Ml I 
AREAl>Xl.HM*m(inAXl-XM !MIIM0RR€0-XQ( IK )) / ( XOl IMAX 1- *0(1*11 
AMC-XAHCI IM)«(XAKC( | MAX t-X ARC U K )) * I OPR E 0- X(J ( J M H/ ( X(jt I MAX J~XO( I W| 
II 

Vlb-XC]FmMI«IXDIFUIMAX|-XO|FmMn*(OPRED-XC( IMM/IXUIMAXI-XGI 
l|M| I 


GU Tu 25 
2t CONTINUE 
762 MKJTF 16*231 
WRITE (6*311 
WRITE (6*701 QFAEQ 
WRITE (6*1631 
WRITE 16*6911 07 
GO TO | 

25 w<UTE (6*31 1 

WRITE (6*70) OPKEO 
«KITE 16*126) 

WRITt (6*132) PEAF 
WRITE (6*26) 

WRITE (6*70) AREA 
AW^AKEAL 

GO TO (620* 62 1*627*620 1 *OTOW 

620 WRITE (6*6261 
lAY£H*(CWD-02l/f2.0t|Ol«02)J* 

XlA-tAYER 

WRITc (6,5551 CMD*LAYER,I!«I2 
XMFLO-RHOL *VC ( D*AW*XLHP+RK3Y*AREA*XLHP 
YOW-(I.O-VOIOI*AW*XLHP 
GO TO 666 

621 WRITE (6,625) 

XMFL0*r<H0L*(AW*XLHP*2.0*Dl*V0ID«H*XlHP|*AAEA*ftWV6XDW 
VOW-Q. 0 
GO TO 623 
627 WRITE (6*626) 

XMFLD-RHOl*AW*XLHP*RHOV6AAEA#XLHP 

YOW-O.O 

623 WRITE (6*70) CWD 
666 WRITE (6,911 

WRITE (6*70) AW 
WRITE (6*93 » 

WRITE (6*70) KKFLD 
WRITE (6,163) 

WRITE (6,691) 0T 

GO TO ( 310,31) ,312,313), OTOw 

310 WTW-AW*XtH?w(Di*RHOWl*D26RHOW2>/(Oi602> 

HAE-HC5l*ARC*A£ 

MAC-HCSL*ARC*AC 

VOW« ( I , 0-VOI D) *AW*XIHP 

GO TO 316 

311 WTW-Rr2)Wl*XLriP*(3.16l596CW062.0*H|*2.0*Ol 
MAE-HEARWAF 

HAC-HCAR«AC 

GO TO M*, 

312 WTM-Rt jW1*XLhP*3.16159*(X 1 0-2* 0*CW0 1 *2* 0+01 
HAE-XRU06AE/CWD 

haC-X*UQ*AC/CWO 
GU TO 316 

313 wrw-3.16139*IXIO-CWDI«UMHOWI6XU(P#2.0*01 
mae-xkeffrae/cwo 

HAC-XREFF6AC/CWD 
316 ETO-OTRAN/HAC 
CTOwTRAN/HAC 
WRITE 16,165) 

WRITE (6,770) VTO, ET0,C1 0 




CALI WT 'NA I DELM.FSAFE , XOD, TH INfOPwT ,KHOL ,RHGM r RMAX, X 1 0, XIHP, KM 
lflD.XMW,lMAX,$TRES,WTlrf,NOUW,OTOW,VTD,t TU,CTO,WALLK.QrRAN,Xlttf,XLCU 
2»HAC ,HAE , $lATk*AS0«SA t BS0»SBi SC, VOW, AW ) 

CO TU t 

AO IF (Nl .G7.2) GO TU AB 
NL-HL*5 
WRITE (6,151) 

GO TO (A|tl»A(12,A!13,AllAI .OTUW 
All! HftITt 16, A63Q) 

WAITE (6,3630) 

GO TO Alls 
A 1 12 WRITE (6, A63 1 ) 

WRITE (6,3630) 

GO TO AU5 
All) WRITE (6, A632 ) 

WRITE (6,A6)3| 

GO • TU Alls 
AUA WHITE (6.A630) 

WRITE (6,3630) 

AllS WRITE (6, Ml 
At Cl "A8SI CHEAO I 

C2-XMULMXIHPHLA0I/I2.0*RH0LAAW*XKP*XLAHD> 

C)»0.06SS*XMUV**0*?5*( XLEV/3*l'**XLA9/L«56+XLCU/$* IA) /IRHQV*RV**l*2 
15*(XLAhD*AV)**l.TS) 

CHECX-0.0 

urui- 2 .o*ci/(c?*cj) 

0 A 1 CHECK»CMECK ♦ 1.0 

*_ IF (CHECK. Gf. IS. 0) GO TO AS 

i H-l 

C/1 A(H)-C)*QT(N)**1.75*C2«0MN)-C1 

— DA(N)»l.7S*C3*QT(N)**0.7S*C2 

IF (OA(N)I A3,AA,A) 

AA 0T|H)-QT(N)4O.O2 
CO TO A1 

A) OMHM)-QT(N) A(NI/DA(N) 

QT ERF t « I MQT IN) 

IF (A BS (t QT (2I~QT (IMS OT fill «LT* 0*011 GO TO AS 

UTIKi-OTIK+ll 

GO 10 Al 

AS IF CI.GT.l) GO TO BIB 
11-11*1 
O-EAO-CPGHO 
GO T J AB 
BIB IMII-l 

CMAXC»QTEMP|l) 

XOIFTI lNDnC)»(,TEMP<lt*M,7S*TEMP/(XLAMD*RH0VI 

OIFfsXOlFTI INDI 

AETO-CTEPMI 11/HAE 

ACTO-UTEXPlD/hAC 

IF UDlFMIMDl.GT.Of ) GO TO 1000 

WHITE (6, BOA) CWO, XOIFTI INOf*AEfO,ACTO,OTEHPI l),0TCHP(2) 

GO TU A2 
A A-<EA*XA| I ) 

CWO'Xrtl 1 I 
AR5-Al«XlUI 
VTIWOIFTIU 
AHC»XAHCI 1) 

UTRA'4>XQ( 1) 


GU TO ?S 
3 WHITE (6,20)1 

GCJ TO 1 

12 READ IS, 6009 ) ARHAX»ARMINiWMAX.W*tN v TlMAX,TLHIH c H£Gft,HCGR 
HE AO ( St A9S ) NOAKtNDWO 
WRITE (6,10) 

WRITE (6,60) 

WRITE (6,61) HD6, HD7,HOB,HD9.HDIQ,HOll ,HD1 3,HDIA,H01S,HD!6 
WRITE (6,62) 

WRITE (6,67) H03,H0A,HDS,HD12 
WRITE 16, 

WRITE (6,66) 5TRES,RH0K,FSAFE 
WRITE 16,66) 

WRITE (6,70) X0D.TM1H 
WRITE (6,60 
WRITE (6,69) 

WRITE (6, TO) ALEV, XL AO, XL CO, XL HP 
WRITE (6.A00) 

WRITE 16.71) 

WRITE (6,701 RMAX 
WHITE 16,10) 

WRITE (6,73) 

WRITE (6.17A) HO6,HD7,HO8,.HO9,HDlOtHDll,HOl),MDlA,HOl5«H0l6 
WHITt (6.13) 

WRITE (6,70 ARMAXvARMIN 
WRITE (6, 1AI 
WRITE (6,70) WMAX,WMIN 
WHITE (6,16) 

WRITE (6.70) TLMAX,TLMIH 
WRITE 16,10) 

WRITE (6,771 

WRITE 16,70) HOI, HD2, TEMP 
WRITE (6,621) 

WRITE (6,622) XMW 
WRITE ( 6, 62A I 
WRITE (6,62)1 P1H 
WRITE (6, A622 ) 

WRITE (6,623) GAMMA 
WRITE (6,721) 

WKITt (6,722) HEGR.HCGR 
WRITE (6,79) 

WRITE (6,.0) SIGMA, XLAMO 
WRITE (6,81) 

WRITE (6,70) PHGL,PHOV 
WRITE (6,03) 

WRITE (6,70) XMUL.XHUV 
WRITE 1 6, 19A ) 

WRITE 1 6, 70) XNL 
WRITE (6,101 
WRITE (6.868) 

WRIT; (6,31) * 

WHITE (6,70) QPAEO 
WRITE 16,126) 

WRITE (6.1)2) PERF 
WRITE (6,889) 

WRITE (6,70) HS 
wRITfc |6,R90I 
WRITE (6,891) (MAX 



« > 


ftrtlTt ( 6,1431 
WRITE ( 6 * 09 M OT 
WRITE ( 6,101 
WRITE ( 6 , 05 ) 

WRITE ( 6 , 06 ) MS 
WRITE 16,161 
WRITE ( 6 , 17 ) 

ONAR-flDAR -1 

ARlNflARUAX-ARMNl/ONAR 
PAR 1 "ARM IN- Aft INC 
DO Si I-l,NDAR 
OAK 1 »OARUAR 1 NC 
0 RRUI- 0 RR 1 

OAwSPf I )■ I -HIGH*S IGN* 1 (H JGH*SIGN )** 2 * 0 .O*$IGNA*CO! 1 (PHI 1 *OAR 1 1 1 / ( RH 
101 * 9 . 0 ) )** 0 . 5 I/( 2 . 0 *DAR( ()l 

oaequi-high*sigh*oawspid*oaai n 

MITE ( 6 , 34 ) OARUI,OAWSMU 
13 CONTINUE 

WRITE ( 6 , 200 ) 

WRITE ( 6,2011 

SONIC - 1 1 CAH»A*TEHMl. 0 * 1 545 . 33 * 32 . 2 /XNWI** 0 . 5 ) * 0.3040 
ASONIC-PEAF*QRAX*( 2 . 0 * 1 GAMNA* 1 . 0 ))**Q. 5 / (XLAHD*RHOV* SONIC) 

WRITE ( 6,2021 SONIC, ASGNIC 

MITE ( 6 , 50 ) H 06 ,H 07 , HD 9 ,H 09 , MOlOfHOll, HOI 3 , HO 14 , HD 15 , HOI 6 
MITE 16 , 35 ) 

WRITE 16 , 36 ) 

WRITE 16 , 37 ) 

DHWD-NOkO-l 

W 1 NC- IWHAX-WNINI/DNWO 

00 in I- l ,NOAR 

PA-WHIN-WiNC 

DQ 1 1 4 J-l.NDWD 

PA-PAtWIHC 

XYZIJI-PA 

IF (XVZ(J).CT.OAWSP(III CO TO 121 
OEPI J)-DAR|l)*XYZt J 1 
3 ftV-( XlD- 2 . 0 * 0 EPI JH/ 2.0 
NNAX- 2 . 0 *SRV* 3 . 14139 /(XY£IJ)*HHIN) 

NHIN- 2 . 0 * 5 R Y* 3 . 14159 / 1 XY 21 J )+TlRAX )+l.O 

00 115 K-NH 1 N,KNAX 

IND-R 

XK-K 

DXUK)» 2 . 0 *SRV* 3 .L 41 S 9 /XK-XYZW> 

6 ASE- 3 . 14159 *XI 0 /XK- 0 XLIK) 

SAV« 3 . 141596 SRV*SRV 
SA 1 * 3 . 14159 *X 1 0 *X 1 D/ 4.0 

A«-(SA!-SAV-XK*DXUK)* 0 EKJI-XK*(O. 3 * 8 ASei** 2 l/XX 

ANO» IO.S- 3 . 14159 / 0 . 0 ) *XYZf J )** 2 

AtO-ACO-AhO 

AW«ALO*XK 

GVC’HiU)-AN 

AL ANO-S A l-SAV»AGO*XX 

WP- 2 . 0 *( 0 EPWI* 0 . 325 * 4 A$€)*RR 

DHL * 4 . 0 *AW/WP 

AV‘SA|-AW'AUNO 

IF UV.U.ASOMC) GO TO 118 

K Y»SK Y 

XX P- ( Oht*DHL 1 / 32.0 


4 PF-XYZIJ) 

t.'GHIW. 0 *SIGNA*CU 5 |PH| 1 /RPE 
<JUOY«ftHOL* 9 » 4 * (HI GM*S IGN* 0 EP(J 1 ) 

CHfcAu-CPGHD-dUDY 

uo ro nr 

3040 Cl»A*S«.CMtAD> 

C 2 -XMUL*(XLHP+XtAO)/l 2 . 6 *«HOL*AW*XKP*XLANOI 

C 3 » 0 . 0655 *XNIJY** 0 . 25 * I XL EV/ 3 . 14 +XLAD/ 1 . S 6 +XLCU/ 3 . 14 ) / I AH 0 V*AY** 1.2 
16*1 XL ANO*AV >** 1 . 75 ) 

CH£Ck* 0.0 

0 UU- 2 . 0 *C 1 /IC 2 *C 3 ) 

3041 CMtCK-CHECK* 1.0 

IT ( CHECK. GT . 20 . 0 ) CO TO 3045 
N-l 

A ( N| >C 3 *QT ( N) ** 1 • 75 *C 2 * 0 T ( N )-Cl 
DA 1 N>-l. 75 *C 3 *QT(Nl**O.T 5 *C 2 
IF IDA(N)) 3043 , 3044,3043 

3044 QMNI-QT(N>* 0.02 
GO TU 3041 

3043 OT(N*l)*QMN)-A|NI/OA(N| 
o:enpmii«qmni 

IF (ABSHQTt 21 -QTC 1 ) l/QT( 1 ) I.LT. 0 . 001 ) CO TO 3045 

QT(N)-QT(N*ll 

GU TO 3041 

3045 IF m.OT.ll CO TO 3010 
11 - 11*1 
CHEAD-CPGHO 

GO TO 3046 
>010 11 - 11-1 

X0(FT(tN0)-C)*QTCHF(l|6*U75*TENR/(XUN0*RHQY) 

IF IXOIFTIINOI.CT.OT) GO TO 110 

DDlfXl-v.O 

EE 1 IRI- 0 TENM 11 

FF 1 (K 1 - 0.0 

CC1(X)-QTEMP(2) 

HM11R)=0.0 

OGN< X ) -RHOL *XLHF*AW*RM 0 V 6 UHP*AV 
OwICX ( X )-RHOM*XLHR*AL AND 
NCUK)-K 
GQ TO 115 

117 0 Dl(X)-O. 5 *CHEA 0 /(XNUL/( 2 . 09 RHOi*AW 6 XKR«XLANDI* 4 . 0 *tNUV/UV*RHOV 6 J! 
1 LAND*RV«* 2 . 0 )I 
tt lift ) -001 (Kt/XLEFF 

F F 1 ( K ) - 0 . 5 *CPGH 0/ 1 xr,)L/ 1 2 . 0 *AHCL*AW*XKP*XIAHD) * 4 . 0 *XNUV/< AVMMOV** 
1 LANQ' RV** 2 . 0 ) ) 

GG 1 tKt«FFHK)/XLEFF 
«EYNn- 2 . 0 *EEMX)*RV/IAV*XNUV*XLAR 0 ) 

IF (REYNQ.GT. 2200 . 0 ) GO TO 3040 

XDIFT(K)-|S 1 0 *XMUV*OOl!K)/(AV*RHOV*XiAND*RV** 2 l)*rEMP/(XLAM*RHOV) 
IF UOIFMKt.GT.OTI CO TO 116 

HHUKI— O.5*(9.0*RHOL*SIGN/XLEFF)/(XHUL/12.O*RHOt*AW*XKP*XLANOI*4. 
10*XNUV/(AV*RHOV*XLANO**V**?.0| ) 

HHlUl-HHlUl*1OOl(X)/A05(OOllRn) 

NG 1 IKI-K 

DGHlK)-RHOt*XLHP*AW»RHOV*XLHP*AY 
UWICMR)»RH 0 W*XIHK*AI.AND 
GO TO 115 
114 NGUKI-R 


C.I-53 




IDIFMM-O.O 
DXUKPO.O 
DD1U1-0.0 
tt If X)*0.0 
FF 1UI-U.0 

CG1UP0.0 
HH1U»*0.0 
UGh(K»*0.0 
DtalCK (K | •0.0 

I IS CUNT J hu£ 
ll*N*IN 

00 119 t’-NHIN'hlU* 

IF M8-EE1UJI 120,119,119 

120 Ll'l 
AB^EUt) 

119 CONTINUE 
NG2I Jl-U 
PB2( ji-xtm jj 
CC2I Jl-OXUUI 
0D2IJI-D01IUI 
t£2 WI»EE MU | 

FF*(J>»Ff MLL ) 
GG2(J)-GGLfU J 
HH2I JI-HHMU I 
FP21 J|*DG»IU> 

002 C J I • Ow I C k ( 1 1 1 
*«21 JI-XD1FTUL) 
CVOW2U)-CVOMMLL1 
GO TO m 

121 NG2IJI-0 
BR21JI-0.0 
CC2 1 11*0.0 
D02IJI-0.0 
EE2I JT-G.O 
FF2I JPO.O 
GG21JI-0.0 
HH2UI«0.0 
PP21 Jl-0.0 
0021 JI-O.G 
*K2( JI-0.0 
0V0W2IJ1-0.0 

119 CONTINUE 
N«* l 

&C-EE2INNI 

00 131 H* 1 , NDtf D 

IF I BC-EE2I N| ) 1)9,133,133 

1)9 MK>M 

BC*CE2(H| 

m CONTINUE 

*721 I >«Bft?im«l 
OGMU ) ■ XYH 1 I 
0*1 (||-Cr.2tNN) 
ill I l'NG2(NN) 

□Gil M-002TNN) 

0021 I l-U 2<NN| 

DOST I >»ff ?(*»! 

00911 MG02IN-I 



DG9I 

OONI n-PP?INHI 
UWIC*<ll-0y2tHH) 

NCAScI J J*t 
ovrom-Afc2ihwi 
GVGW3 I I } >CVQm2(HN) 

11) CONTINUE 

00 139 IN- 1 ,NGmR 
IH»IN 

C0*DG2( IN| 

00 137 IKNMN.NOA* 

IF (C0-DG2I INNI) 138,1)7,137 
138 IM-Jmh 

CD-0G21INN) 

1)7 CONTINUE 

NST&a-NGI INI 
ST0N1-0GW1 IM 
STOP 5*001 (INI 
STOft9>DG2l JM| 

ST085-0G3UNI 
STOR6-OG9I IN} 

STOR7.DOM1HI 
STrjftft.DGNI INI 
STQR9«0 UICk| |M) 

HGIO-OAM INI 
STOU-DXlflH) 

STQ12*DVT0( IN) 

S TUI 3»GVUN 3 ( IN | 

0VTD1 IN1-0VT0IIN1 
NO! 1H|*NG( IN) 

DGwt INI-DGNl INI 
0G1 (INI-OGHIN) 

DG2I INI-0G2I INI 
0031 1*1*0031 INI 
0G9I INI-DG9UN1 
0G5I INI-0G51 INI 
OGN( IN)»DG“( IN I 
ON IC* 1 INI "QUICK I IN) 

OARI IH|»DAR{ INI 
OXMINPOXK IN) 

GVUNJ(|N)*GV0N3IIN» 

NG( INPN5T0R 
OONI INI-ST0R1 
0011 INI-STUt) 

0G21 lN|«S r 0P9 
UG3I IM*STURO 
0G9( INl^TOP.9 
DG*< IN)«SfClR7 
LGH( I'.|»STURh 
D wICM IHl-STOkV 
OAJII INI-GTOlO 
UXLI I N I • G IU 1 1 

OWTOUNPiTOl? 

OVUWJl IN) » r ,TfJl 3 
139 Cunt j mji 

UfJ 9 K PI, SOAR 

1 7 1 I [oOM i n l***G(MtCAB|| I «OGUt 1 1 •OKI IIS«f)GI|| I t DG2 1 1 1 *0G3 


t>Q * ro 


i ■ 9* CONTINUE 

'} *KITE 16, ill 

ii WRITE (6,70) QPHl D 

i WRITE 16.1261 

i’i WRITE <6,1)21 PE*I 

LCA5E-0 

f 00 16\ 1*1, NONA 

IF (002111. IT. CPRE01 GO TO 162 
LCASt-I 

161 CONTINUE 

162 IF UC6SE.EQ.G) GO TO T62 

1 MRITt <6,661 

•; NR lit 16,66) 

j ’ NR 4 T.6 <6,661 

jj'i WRITE 16.6M 

It MAE-HEGR*AE 

,L ' HAC-HCGR*AC 

00 60 1-l.LCASE 
ETD-0G2UI/HAE 
CTD-0C2( 1 )/HAC 

WRITE 16,69) NCA$EI I ) ,NGI 1 ) ,OAR| 1 1 ,OGtf 1 1 ),OXL f I ) , 0G2 < 1 1 »DGHl I ) ,OVT 
lO<n,fcTO,CTD 
6« CONTINUE 
STQRt-STRES 
00 96 I-l.LCASE 
STRES-STORE 
VT0»0/T0<U 

wTN.uwicmn 
NCAAE-t 
XMELD-DCMII) 

V0W«-GVGW3(ll 
ETD-DG2 < I l/HAE 
CT0-DG2(I)/HAC 

CALL NT ANA IDElN,FSAFE,XaO»TMIN.OPWT # AHOL,AHOM t RMAX,XID,XLHP,XM 

lFLD.XHW,TMAX.STRES,WTW,NCAAE.QTO‘SVTQ.ETO,CiD,WALLK*QTRAH,XLEV,KLC 
20,HAC,HAE,STATE,ASU, SA, ft SO, SB, SC, VOW, AM ) 

96 CONTINUE 

1 IF <M0RE.GT.1I CO TO 600 
STOP 
END 


SUBROUTINE WTANA < DEL«, FSAFE, XOO, TNI N,OPNT,RHOl ,KM0M,RMAX,X10.XLHP 
l,XNFLD,XHW,TEMPtSTRtS»WTU,NCAAE,OTOW,VTD,ITD,CTD,WALLX,OTRAN,XLEV, 
2XLC0,HAC,HAE»Sf ATE,ASO,S BSQ, SB, SC, VOW, AW) 

DIMENSION THICX<20),STH!K<20) 

INTEGER 0PWT.0T3W, STATE 
B7 FORMAT 1*1*) 

01 FORMAT I////* WEIGHT ANALYSIS FOR THE 
1 S Y S T E MM 

)6 FORMAT <//• (EQUATION OF STATE FOR REAL GAS < I .E. ,BEATT IE-BR IDGmA 
IN EQUATION) IS USED TP CALCULATE THE INTERNAL PRESSURE)*) 

37 FORMAT I//* IEOUATION OF STATE FOR IDEAL GAS IS USED IQ CALCULATE 
I THE INTERNAL PRESSURE ) * 1 
6 FORMAT I//25H CASE NUMB6 R,AX,|1) 

66 FORMAT iff* SAFETY FACTOR MAXIMUM Tf MPf RATURt < K ) • ) 




66 FORMAT l/FI1.2,6X,F16.2l 

66 FORMAT <//• LIQUID NT. (KGI WICK WT.IKGI • ) 

66 FORMAT < /E 1 6 . 7, E ) 8. 7 ) 

2) FORMA I iff* LIQUID WT.UGIM 

TH FORMAT iff* RAO. I STORAGE ) TUBE THICKNESS Sf. THICKNESS PRESSU 

UF HEAT PIPE wT • STORAGE WT. TOTAL WT.M 

79 FORMAT <• IN) IM) <H) <N/M 

l-H) <KG) I XG| <KG)M 

303 FURMAT (/7EI5.7I 

92 FURMAT iff * NO STORAGE VOLUME IS NEEDED* I 

11 FURMAT <////• THERMAL ANALYSIS FOR THE 
1 S Y S T E M *| 

13 FORMAT iff * INSIDE DIAMETER uF THE TUBEINI • I 

16 FORMAT iff * EVAPORATOR LENGTHIM) CONOENSER LENCTHIM) • I 

15 FURMAT </El5.7,l2X,El5.7! 

17 FORMAT iff * THERMAL CONDUCTIVITY OF THE TUBE WALK W/M-K) • I 
19 FORMAT iff* RAOIAL HEAT TRANSFER ACROSS THE TUBE WALL IM) * ) 

25 FORMAT iff* VAPOR TEMP. OROPtK) EVAPORATOR FILM TEMP. OROPtR 

1) CONDENSER FILM TEMP. OftQPIKPt 

26 FORMAT < /E 1 5. 7 , 1 6X ,E 1 5.7,23X,E 1 5.7) 

27 FORMAT iff* RAO. < STORAGE) TUBE THICKNESS TEMP. DROP ACROSS 

1 TfcMP. DROP ACROSS TOTAL TEMP.* I 

29 FORMAT <• (HI (M) THE EVAPORATOR WALL IK) 

1 THE CONDENSER WALL I K 1 DROPIK)*) 

62 FORMAT (/2EIS.7t2X,E15.7,!6XiEL5.7«6X«E!S.T) 

IF I0PWT.LI.1I GO TO 111 
WRITE 16,8?) 

WRITE (6,08) 

STHlKCll-0.0 

THICK11I-TMIN 

NNN*0 

CONST-0. 00206 

II-l 

IND-0 

BAW-AW*1D«76 

RS-0.0 

WTS-0.0 

DELR-0.0 

0DELH*DELH*3.2B1 
0RHOH-RHOM*O. 06246 
0RHOL-RHOL*O. 06243 
BR-X 10*3-281/2.0 
BX (0-X 10*6.281 

STRES-<STRES/( 706.1*9.0) )• 166.0 

BXLHP-XLHP*3.201 

BRMAX-RMAX*).261 

0X00* XOO* 3. 20 l 

BTEMP-TEMP*l.e 

BTMIN-TMIN*3.701 

BXMFL»XMFLD*2.205 

IF (STATE. EQ. LI GO TO 32 

WRITE (6,36) 

SPVGL- 1 1 3. 1 6 L 59*XIO**2/6.0)*XLHP»1000.0-VOW)/< XMFLD*1000.0/KMW) 

A- A SO* ( 1.0-SA/SPVUL I 
B*BSQ* | l.O-SB/SPVQL I 
EPSI«SC/ISPV0L*TEMPA*3) 

PRf S-CONST*TEMP*< 1.0- EPS I )•< S*VOL*0)/< SFV0L**2)-A/<SPVOL**2) 
0PRES-PRES*16. 7*166.0 



CO TO 33 

3 2 6PRES-&XHFl.*134S.li*PTewP/(3.14139*ttH*BK*exLHP*XMW-V0W*XHw| 

WRITE 16,37} 

33 BtREO*BPRES*OR*FSAFE/STRFS 
IF 16TRED.GT.BTM1N) GO TO 16 
INO* 1 

BTRED--BTN1N 
16 IR£O*6TREO*O,3O60 
THICK ( l 1* TREO 

dVHP* 3. 1 4 139*BXLHP*BR* BR-VOW 
P*E5»6PRES*7Q3. 1*9.6/ 14-.0 
WTL-XHFLD 

WTH-HHOM*3. 14159*1 X I 0*7 AED*TKE D**2.0 )*XLHP 
IF (QT0W.EQ.5) GO TO 100 
WT*WTH*WTL*NTh*WTS 
GO TO 101 

100 WTM«WTH*WTW 
WT-WTH*WTL*WTS 

101 IF (0T0W.NE.5I GO TO 3 
WRITE (6,41 NCAAE 

3 WRITE (6,461 

WRITE 16,43) FSAFE.TEWP 
IF (OTOW.EO.31 GO TO 21 
WRITE 16,64) 

WRITE (6,63) XKFLOtWTW 
GU TU 22 

21 WRITE (6,23) 

WRITE (6,63) XMFLD 

22 WRITE (6,78) 

WRITE (6,79) 

WRITE (6,303) RS, TREO, OELR, PRES, WTH, WTS,WT 

IF (QPWT.LT.l) GO TO 2 

IF ( INO.EO.ll GO TO 91 

BINCR-ICRHAX-BR 1/9.0 

IF (BINCR.LT. 0.0) GO TO 2 

8YHP-BR*BR*3.14159*BXIHP-V0W 

BAHAX*8R-BINCR 

(JO 30 1-1,10 

0RMAX«&RHAX*R1NCR 

AS*BRWAX*0. 3048 

QHl»0RNAX-(BRHAX**2-BM*2l**0.3 

&VST-4.0*1.14139*6AHAX*M/3.0-3.l4139*bHl*11.0*BH**2*BHl**2)/6.0 

bVTT*bVHP*bVST 

IF (STATE. EO. 1) CO TO 34 

SPV0L«BVTT*0. 02632*1000. 0/(XHFLO*lOOO.O/XHW) 
a«ASl*11.C-SA/SPVDL) 

0-BSU*( l.O-SH/SPVOL ) 
f PSl-SC/(SPVfJL*IEWP**3) 

PRES= CONST *TEHP*( I .O-EPSl )•( SPVOL*fl ) / ( SP V0L**2 )~A/ t SP VUC**2 ) 
HPrtES‘PRES*14. 7*144.0 
Ml IU 33 

14 tJPRES-HXNFL • 1 343 .T1*RU HP/ <BVTT*XNN I 
33 BTRFC-t;PR£S*RK*FSAFE/STRES 
IF (BTrfED.GT.ftTXlN) GO TO 31 
HTP.EO-PTMIN 
NHN= 3 

31 «0tlk*nPKfS*HRMAX*FSAF6/l2.0*STRfcSI 
IF IriCFLR.GT.hDELH) CO TO 32 


ODELR-BOELN 
32 RHAX=6RHAX*0.304B 
TREO*flTRE 0*0. 3048 
0tLR*«DELR*0.3O48 
11*1*1 

STHIKI II |*RMAX 
THICK ( 1 1 1* TREO 
PRtS-UPRES*703. 1*9.8/ 144.0 
HI *BM1*0# 304B 

WTH-RHOH* 3. 14159*1 XI0*TRtD*TRED**2l*XLHP 

WTS- (4.0*3. 14139*(RNAX*OELR1**3/3.0-4.0*3.I4139*R)(AX**3/3.0)*RHO«- 
12. 0*1. 141 39*RWAX*HI*DELR*RH0N 
IF ICT0W.EQ.3) GO TO 30 
WT*LTH*WTL*WTW*WTS 
GO TU 31 

30 WTH*WTH*WTW 

WT *WIH*WTl *WTS 

31 WRITE (6,3031 R$,TREO,DELR,PAES, WTH, WTS.WT 
IF (NNN.GT.U GO TO 2 

30 CONTINUE 
GO TU 2 

91 WRITE (6,92) 

2 WRITE (6,87) 

WRITE (6,11) 

IF ( OTOW.NE .3 1 GO TO 12 
WRITE (6,4) NCAAE 
12 WRITE 16,13) 

WRITE (6,303) XIO 
WRITE (6,14) 

WRITE (6,131 XLEV.XLCO 
WRITE (6,17) 

WRFTE (6,303) WAUK 
WRITE (6,191 
WRITE (6,303) QTRAN 
WRITE (6,23) 

WRITE (6,26) VTD,ETD«CTD 
WRITE (6,271 
WRITE (6,29) 

DO 28 J-1,11 
DEXT-XID*2.Q*THICK(J) 

OlNT'XlD 

AEXTE ”DEXT*3. 14) 39*XLEV 
AfXTC«OEXT*3.14l39*XLCO 

E Wl D a QTRAN*OEXT *( ALOGI DEXT/OINT I )/( 2.0*NALLK*AEXTE I 
CNT0-QTHAN*OEXT*(AL0G(DEXT/0|NT))/(2.O*WAUK*AEXTC) 
TTD-VTD*ETO*tTD*EWTO*CWTO 

WRITE (6,82) STHIK( J ) .THICK (J) ,E WTO, CWTO, TTD 
28 CONTINUE 
111 RETURN 
END 
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The first sample problem consists of determining the heat transport capa- 

-2 -4 

bility of an aluminum (6061-T6) heat pipe (1. 27 x 10 m O.D, x8.9xl0 m wall). 

A homogeneous circumferential porous wick (200-mesh stainless steel screen) is used 
in the heat pipe. The heat pipe working fluid is nitrogen at an operating temperature 
of 80°K, and the tube wall material is aluminum alloy 6061 -T6. Heat load, elevation, 
maximum system temperature, and acceptable range of vapor temperature drop are 
specified as: « 


Q 

^req’d 

h 

T 

max 

AT V 


5W 

0.00254 m 
300. 0°K 
10. 0°K 


i 


The length of the evaporator, adiabatic section and condenser of the heat pipe are 
specified as: 

L = 0.2 m 

e 

L = 0.8 m 

a 

L = 0. 2 m 

c 

The spherical storage reservoir, if needed for containment purposes, is specified to 
-2 -4 

be 5. 8 x 10 m in diameter and 8.9 x 10 m wall thickness. Finally, the effective 
thermal conductivity of the wick is assumed to be 0.14 W/m-°K (i.e. , thermal con- 
ductivity of nitrogen). The associated data cards are listed in Table C. 1-3. The 
resulting computer output data follow. As indicated in the printout, the maximum 
heat transport capability in a one "g" environment is approximately 0.98 watts, which 
is less than the required heat load. Therefore, weight and heat transfer analyses are 
not presented and the program terminates with the following statement: "No area exists 
to satisfy QMAX (w) requirement". 

In the second sample problem, the performances of a nitrogen heat pipe 
with an axial rectangular grooved wick are considered. The maximum and minimum 
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allowable aspect ratios, land thicknesses, and widths of the grooves are specified as: 


^ max 

— 

1.5 


* . 
mm 

= 

0.5 


L 

max 

= 

5.10 x 10' 4 

m 

L min 

= 

3. 81 x 10" 4 

m 

W 

max 

= 

7. 62 x 10 -4 

m 

W . 
min 

= 

4.00 x 10~ 4 

m 


The heat load requirement is: 


Q, 


max 


0.005 W 


2 o 

and the evaporator and condenser film coefficient are assumed to be 380 W/m - K 
and 760 W/m 2 -°K, respectively. These values are based on measured data. The 
other pertinent data remain unchanged and are specified in the first sample problem. 
It should be noted, however, that the properties of nitrogen are input here because of 
the option used in the control card. The associated data cards are listed in Table 
C. 1-4. The resulting computer output data follow. Since the maximum heat transport 
capability in a one "g" environment of each specified aspect ratio is larger than the 
heat load requirement, the weight and heat transfer analyses are performed for each 
aspect ratio. The results of these analyses are presented in the computer output. 
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C. 2.1 


DATA ACQUISITION CODE USER'S MANUAL 

ft 


Introduction 


This section describes the utilization of the Data Acquisition Code (DAC). 
Basically, DAC generates the thermophysical properties of various heat pipe working 
fluids at saturation conditions. The fluids selected are those most commonly used and 
span the range from cryogens to liquid metals. Both constant properties and tempera- 
ture dependent properties are stored in the property subroutines of the DAC. Proper- 
ties are stored as tabular data points or, in the case of the liquid metals, in functional 
form. Also, derived properties (e. g. , liquid transport factor) are calculated and out- 
put by the code. 


* This code can be used alone to obtain fluid property data or can be used in 

combination with the Heat Pipe Analysis and Design (HPAD) Code and thereby reduce 
the amount of input data required. 

The program input requirements are described in Section C. 2. 2. 2, and 
the output formats are described in Section C. 2. 2. 3. References are listed at the end 
of Section C. 3. The flow diagram and program listing are presented in the Appendices. 
A list of Fortran names with physical or engineering quantities is also presented as an 
Appendix. 


C. 2. 2 Program Description 

C. 2. 2. 1 General 

The program logic is illustrated in the flow diagram contained in Appendix 
A. The FORTRAN names and the physical quantities they represent are listed in Ap- 
pendix B. Storage requirements are on the order of 50,000 words (octal). A listing 
of the program is presented in Appendix C. The program was written in FORTRAN V 
and was designed to operate on the UNIVAC 1108 system. 

Basically, the program reads the input control data, generates the thermo- 
physical properties of the specified fluid, and outputs the property data. The deck set- 
up as shown in Figure C. 2-1 consists of job control cards, the program source deck, 
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> . .. 



Control 



Figure C. 2-1. Program Deck Setup 


C. 2-2 






additional control cards followed by the input control card and program termination 
cards. The program has two major options. First, the thermophysical properties 
can be determined at a specified temperature or at specified intervals over a specified 
range. If the specified temperature range is out of property data limits, the program 
will automatically adjust the range to be consistent with the available property data. 
The second option consists of having the property data output in scientific units (SI) 
or in both scientific and engineering units. 

C, 2, 2, 2 Input Description 

Table C. 2-2 describes the entries to be made on the various input cards 
and indicates when each of the optional cards are to be excluded. The FORTRAN name 
format, and units to be used are indicated for each entry. A listing of sample input 
data for different cases is presented below in Table C.2-1. 

Table C.2-1. Sample Input Data 

2 

4 

1 2 

3 k ! 0 • 0 

u 

4 

2 1 

2<sO*0 3l?0, 1C, 


C. 2. 2. 3 Output Description 

The program outputs essentially all the thermophysical fluid properties 
as well as derived properties needed for the Heat Pipe Analysis and Design Code. 
These properties are divided into two categories in the output. First, the constant 
properties which consist of the following: 

• Molecular weight 

• Gas constant 

• Ratio of specific heats 
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Table C. 2-2. Input Data Description 


Input 

Card 

No. 

Format 

Fortran 

Name 

Description 

Units 

1 

13 

MORE 

Control Point, integer = 0 for last set of 
data, otherwise integer = 2 

■* » 

2 

13 

FLUID 

Control Point, type of working fluid 

- 


.1 Hydrogen 

2 Nitrogen 

3 Oxygen 

4 Water 

5 Ammonia 

6 Methanol 

7 Acetone 

8 Freon-21 

9 Sodium 

10 Potassium 

11 Lithium 

12 Mercury 


3 

213 

POINT 

Control Point, integer = 1 for Arid proper- 
ties at a temperature, otherwise integer = 2 




UNITS 

Control Point, integer = 1 for M.K.S. units, 
integer = 2 for both M.K.S. and Engineering 
units 

- 

4* 

F10.4 

TEMP 

Temperature (for fluid properties) 

°K 

5** 

3F10.4 

TMIN 

Minimum temperature (for fluid properties) 

°K 



TMAX 

Maximum temperature (for fluid properties) 

°K 



DT 

Temperature increment (for fluid properties) 

°K 

6*** 

3F10.4 

TMIN 

Minimum temperature (for fluid properties) 

°F 



TMAX 

Maximum temperature (for fluid properties) 

°F 



DT 

Temperature increment (for fluid properties) 

o 

F 


Card 4 is needed only when the value of card 3 (1) is an integer equal to unity. 
Card 5 is needed only when the value of card 3 (2) is an integer equal to two. 


Card 6 is needed only when both values of card 3 are equal to two. 
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• Normal melting point 

• Normal boiling point 

• Critical temperature 

• Critical pressure 

The second category listed below consists of the temperature dfjpendent properties and 
includes the derived data. 

• Vapor pressure 

• Liquid density 

• Vapor density 

• Surface tension of liquid 

• Latent heat of vaporization 

- • Thermal conductivity of liquid 

• Dynamic viscosity of liquid 

• Kinematic viscosity of liquid 

• Dynamic viscosity of vapor 

• Kinematic viscosity of vapor 

• Ratio of kinematic viscosity (vapor/liquid) 

• Wicking height factor 

e Liquid transport factor 

• Sonic heat flux 

• Sonic velocity 

A listing of output data associated with the sample Input data of Table C. 2-1 is pre- 
sented in Table C. 2-3. If the specified temperature is outside of the data range it will 
be indicated in the output by "TEMPERATURE OUT OF TABLE RANGE". 
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Table C.2— 3» Sample Output Dato 
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ChiMCU fonn\ti,k • H ao 
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NORMAL HLLT|N% rot 
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yoj« f>o 


it nrCMAlORt 
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MfOM PRESSURE 

irsn 


LU'JlO Ul'iSlTT 

iuh/i t-f wn 


fCMRCRATUMi.pLPLNOCMT PROPERTIES 


VAPOR OENSITT 
itRH/r i-» wtl 


SURFACE TENSION 
iLhf /» » I 


.__a*J-H2 


tATCM *tAT 
teiu/LiJM 

IhiKHAU CONDUCTIVITY or UttlUtt-. 

l» tu/nR-t r-r l 

OTfc*H|C RIStOSITT OF LIQUID 

ILIH/FI-HII 

MH h *Uc W|icOilTT Of t|QU|Q_ 

IFI-FT/hHj 

DYNAMIC VIS COSlft or VAPOR 
Ufat/FT-Hit) 
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lRiu/M«-r t-m l 
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Table C.2-3* Sample Output Data (Continued) 
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«.» A T A ACttOJSITlON C M t lD«(l 
Mt K» St UNITS 

FLUID »WATER CHEMICAL FORMULA »H26 

C 0. N. S . T , A .M T F R 0 F C B T It S 


, O 

-S3 

ix 


w 

■'9 

8 


? 8 



molecular height ikg/hoLC) gas co n stantin-m/k&.k > 

I **000 *IA| tUGROlOl 


RATIO of SFCCIFIC HEaTS 
1.330 
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Table C.2-3 ! Sample Output Dato (Continued) 
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O 
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o 

3 
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o 
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o 
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fO 
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II t'l'LKATUKt 

lifTl ¥ !*.» 

7*n.p 

*♦*•0 

JuO * 0 

310*0 

320*0 

330*0 

3' n«r 

3b||.|> 

VAk(jk 

1 n/m-« ) 

•♦7*u3 

• l *♦ U* 

•3b*0* 

•b**U* 

• iu*ns 

•14*0% 

•24*0S 
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•70-01 

•11*00 

• 1 7*00 

• 23*00 

surface UnsIO* 






1 n/h) 

• 7*-o i 

.71-01 

.71-01 

• 4***1 1 

•44-01 

•44-01 

•4*-0| 

•43-01 

lATChT hI*t 

•2S*qT 

•24*07 

•7**07 

• 7*4«07 

•2**0/ 

• 1**07 

•23*07 

•23*07 

iHtpniL condultivitt or uuwid 
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* 4 | *00 
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OtN|»lC V I SC Ob 1 T ▼ 01 LlOUIO 
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.43-03 
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•b«-03 

•*♦-03 

•*3“01 
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0 Hikh|{ VlbCO^HT or V*PQ* 

•*7-Qb 

•|0-q* 

• | 0*04 

MI-0* .. 
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RlktMTK mCQ&l" Or v*P0* 
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•*4.07 

• | 7*0® 

•7o*ot 
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FLUID SELECTION* 
1 HYDROGEN 

7 

ACETONE 

2 

NITROGEN 

8 

FREON -21 

3 

OXYGEN 

9 

SODIUM 

4 

WATER 

10 

POTASSIUM 

5 

AMMONIA 

If 

LITHIUM 

6 

METHANOL 

12 

MERCURY 


CALL PROP SUBROUTINE 
CONSISTENT WITH THE 
TYPE OF FLUID 


LAST SET 
OF DATA 

"'''''''Tyes* 
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INO*l 

UNITS COi<T CARD' 
ENCRG UNITS 


RE AO CONTROL CARO 
FOR TABLE ft UNITS ) 


CALC PROP er 
FORMULA or 

interpolation 



IND * I 

UNITS CONTCARD = 
ENCRG UNITS 
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TRANTER -'f JTi\ENT 
DATA FROM ThE 


prop subroutine 


transfer esc rrwp- 

I'E PENDENT PROPERTIES 


»HOM PROP SUBROUTINE 


CACC TkC DEBivfO 
TEMP-OCPENOENT 


TRANSFER PERTINENT 
DATA FROM THE 


PROP SUBROUTINE 




UNITS OR ENGRG 



UMTS OR ENGRG 


TABLE RANGES 
EXCEF05 TABLE 


table range a 

EXCEEOS TABLE 



SUBROJT.NE POUT 


SUBROUTINE POUT 2 


SUBROUTINE POUT 3 








Appendix B. FORTRAN Names 
and Associated Physical or Engineering Quantities 



Fortran Name 


Description 


Units 


MORE 

FLUID 


LOUT 

LSUPP 

ARHOL 

ARHOV 

A XL AM 

ASIGM 

AXMUL 

AXMUV 

TEMP 

APBESS 

XMW 

GAMMA 

ORHOL 

ORHOV 


Contro; Point, integer = 0 for last set of data, 
otherwise integer = 2 

Control Point, type of working fluid 


1 

Hydrogen 

2 

Nitrogen 

3 

Oxygen 

4 

Water • 

5 

Ammonia 

6 

Methanol 

7 

Acetone 

8 

Freon-21 

9 

Sodium 

10 

Potassium 

11 

Lithium 

12 

Mercury 


Out of temperature range indicator 
Control Point, suppress DAC external output 
Liquid density 
Vapor density 

Latent heat of evaporization 
Surface tension 
Dynamic liquid viscosity 
Dynamic vapor viscosity 
Temperature 
Vapor pressure 
Molecular weight 
Ratio of specific heat 
Liquid density 
Vapor density 


kg/m 3 
kg/m s 
W s/kg 
N/m 
kg/m s 
kg/m s 
°K 

N/m 2 

kg/mole 

kg/m 3 
kg/m 3 




Fortran Name 

Description 

Units 

OSIGM 

Surface tension 

N/ni 

OXMUL 

Dynamic liquid viscosity 

kg/m s 

OXMUV 

Dynamic vapor viscosity 

kg/m s 

OPRES 

Vapor pressure 

N/m* 

OXLAM 

Latent heat of evaporation 

W s/kg 

POINT 

Control Point, integer ~ 1 for data at a tempera- 
ture, otherwise integer - 2 

- 

UNITS 

Control Point, integer - 1 for M.K.S. units, 
integer - 2 for M. K. S. and British units 

- 

HD 

Headings (working fluid and chemical formula) 

- 

GAS 

Gas constant 

ft lbf/lbm°F 
atm liter/gm K 

TMELT 

Normal melting point 

o o 
K, F 

TBOIL 

Normal boiling point 

o 

o 

TCRIT 

Critical temperature 

o o 
K, F 

PCRIT 

Critical pressure 

N/m 2 , psi 

NK 

Control Point 

- 

OTEMP 

Temperature 

°K 

OTHCL 

Thermal conductivity 

W/m°K 

OXNUL 

Kinematic liquid viscosity 

raVs 

OXNUV 

Kinematic vapor viscosity 

m* /s 

RATIO 

Ratio of kinematic viscosities (vapor/liquid) 

- 

OXNL 

Liquid transport factor 

W/m 2 

OH 

Wicking height factor 

m 2 

NBRIT 

Counter for units used 
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Fortran Name 

Description 

Units 

SONIC 

Sonic velocity 

m/s 

OSHF 

Sonic heat flux 

W/m* 

BTEMP 

Temperature 

•A. 

o 

F 

BPRES 

Vapor pressure 

psi 

BRHOL 

Liquid density " 

- __ lbm/ft 5 



BRHOV 

Vapor density 

lbm/ft 5 

BSIGM 

Surface tension 

Ibf/ft 

BXLAM 

Latent heat of vaporization 

Btu/lbm 

BXMUL 

Dynamic liquid viscosity 

lbm/ft hr 

BXMUV 

Dynamic vapor viscosity 

lbm/ft hr 

BXNUL 

Kinematic liquid viscosity 

ft 2 /hr 

BXNUV 

Kinematic vapor viscosity 

ftVhr 

BXNL 

Liquid transport factor 

Btu/hr ft* 

BH 

Wicking height factor 

ft* 

BSONI 

Sonic velocity 

ft/hr 

BOSHF 

Sonic heat flux 

Btu/hr ft 2 

BTHCL 

Thermal conductivity of liquid 

Btu/hr ft °F 

TMI3 

Storage space for temperature 

o 

X 

o 

TMA4 

Storage space for temperature 

O o 
K, F 

T1 

Storage space for temperature 

O 

o 

*1 

T2 

Storage space for temperature 

Cn 

o 

X 

0 

XTEMP 

Temperature 

o 

K 

RHOL 

Liquid density 

kg/m 3 
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Fortran Name 




Description 

Units 

RHOV 

Vapor density 

kg/m 3 

SIGMA 

Surface tension 

N/m 

XLAMD 

Latent heat of vaporization 


W s/kg 

XMUL 

Dynamic liquid viscosity 

kg/m s 

XMUV 

Dynamic vapor viscosity 

kg/m s 

THECL. 

Thermal conductivity 

W/m °K 

PRESS 

Vapor pressure 

N/m* 

MM 

Number of data points 


IND 

Counter for output format 


NIND 

Counter for interpolation method 


TMIN 

Minimum temperature 

o 


K 

TMAX 

Maximum temperature 

°K 

DT 

Temperature increment 

°K 

A1 

Storage space for temperature 

°K 

A2 

Storage space for temperature 

°K 
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zz-z'o 


INTEGER FLUID 
2U READ (5*101 MURE 
RfcAU (5,101 FLUlf\ 

10 FOR HAT (3131 ■£> 

LUP r v; 4* 

lsupp*q ^ 

CALL DAC ( FLl^O, LSUPP.l OUT, ARHOL 

1 MUV ,TfcMp,APRfcS» XMW, GA^M A I 
IF (HCRE.GT.l) CO TO 
STOP 
tNO 


• GAVMA I 

ru 


■RH0Y,AXLAM,ASIGM,4XHUL,AX 


1 S I GH.UXHUl , OX MUV, T E MP, OPK E S, XMW , GAMMA ) 
303 AKHOl*ri«Hni(U 
ARHOV*OKHf iV I 1 1 
ASIGH-OSIGMt 11 
AXMUL-UXMULC II 
AXMuv*nxmjvm 
APRF5*GPKPSU I 
AXL AH = OXL AM ( 1 | 

105 RETURN 
END 


SUBROUTINE OAC (FLUlD,LSUPP,lUUT,ARHOL,AAHOV,AxLAM,ASIGN,AXKUL,AX 
1MUV.TEMP, APRES.Xh ..GAMMA) 

01 HENS ION ORHOL<ei,ORHOV(«l,OflGM 81 ,0XHUU«1,0XHUV(8I *0PRtS(Bl 

DIMENSION QXLAMm 

INTEGER FLUIO, POINT, UNITS 

IF U5UPP.lt. 21 GO TO 301 

NCR 1 1 -0 

POINT-! 

UNITS* 1 

301 GO TO (U2,3,6,5,6,7,fl»9,10.U,12) t FLUl0 

1 CALL HYDRO I LSUPP, l OUT, HER 1 T , POINT, UN I TS »URHDL ,0RH0V,QXLAH,0 

1 SIGH, 0XHUL,0XMUV, TEMP, OPR ESfXHW, GAMMA | 

UO TO 303 

2 CALL NITRO ILSUPP.LGUT ,NB* I T.PCI NT , UN I TS, ORHOL •ORHOV.OXLAM.O 

ISI CM, OXMUL »OXMUV, TEMP, OPRES i XMW, GAMMA I 

GO TO 303 

3 CALL OXYGE (LSUPP, LOUT ,N8RI T, POINT, UNI TS.ORHOL .QRHOV, OXLAM,U 

lSIGM,OXMUL t OXMUV, TEMP, QPRES, XMW, GAMM) 

GO TO 303 

5 CALL WATER ILSUPPtlGUT ,N0RI T, POINT, UNI TS.ORHOL, URHOV.OXLAN.O 

ISICM.OXHUL.OXHUV, TEMP, OPRES, XMW, GAMMA I 
GO TO 303 

5 CALL A*MON (LSUPP , LOUT ,NBRI T , PU INT, UN I TS,ORHOL ,ORHOV, OXL AM,0 

IS! GM,OXHUL,OXMUY, TEMP ,0PRES,XMVJ, GAMMA 1 

GO TO 3D3 

t CALL HETHA (LSUPP, LOUT, NBA I T, POINT, UN I TS,ORHOL,ORHOV,OXLAM, 0 

1 SIGM.OXMULwOXMUV, TEMP, OPRES,XMW, GAMMA 1 
GO TQ 303 

7 CALL ACETO (LSUPP, LOUT, NBA I T, POINT, UNI TS.ORHOL, QRHOV, OXLAM.O 

ISIGH.QXHUL tOXHUV, T6HP , C’KFS, XMW, GAMMA | 

GO TO 303 

6 CALL FRE21 ILSUPP,t0UT,NBRIT,PUlNT f UNlTS,0RMQL.0RH0V,0XLAM,0 

1SIGM,0XMUL 9 0XHUV,TERP,0PHES 9 XMW,GAHMA| 

GO TO 303 

9 CALL SOOIU USUPP, LOUT, NBRIT, POINT, UNITS, ORHOL,ORHOV,OXLAM,0 

ISIGH.OXHUL.OXHUV, TEMP, OPRES, XMW, GAMMA | 

GO TO 303 

10 CALL PUTAS (LSUPP, LOUT.NBRIT, POINT, UMTS,CRHOL,ORHOV,OXLAH t O 

1 S|GN,OXMUL,OXMUVt TEMP, QPRfcS.XM*, GAMMA 1 
GO TO 303 

U CALL L I T H I (LSUPP, LOUT, NBRIT,PCINl.UN|fS,OKHOL t aRHnv, OXLAM.O 

ISIGH.OXHUL ,OXHUV,TEHP,OPRES,XHWf GAMMA) 

GO TO 303 

12 CALL MIRCU I L SUPP .LOUT ,NBR 1 T, PL I NT ,UN1 TS»l)RMOL , ORHDV, OXL AH,0 


SUBRlUVl: t POUT KHO, XKW ,THELt ,TBOIL • TCR ITfPCRl T »NBR| T .GAMMA.LSUT 
DIMENSION HD(2«2) 


35 FORMAT |*1M 

36 FORMAT (////• 

37 FORMAT I / ///9H 
55 FORMAT | //• 

Irt FORMAT UUi* 

1 S*l 

39 FORMAT U///* 
1 


DATA ACQUISITION CODE IOACI 
FLUID -,2A6,9X. IBHCHEH1CAL FORMULA »,2A6> 

M. K. S. UNITS* I 

CONSTANT MOimi 


GAS CONSTANT IN-M/KG- 


1 BM-F | 
SI F0RMA1 


MOLECULAR WEIGHT (KG/MOLE) 

RATIO OF SPECIFIC HEATS* I 

50 FORMAT I/FI5.3.22X, F15.7, I7X.F15.3) 

51 FURMA*’ I iHf* NORMAL MELTING POINTUI NORMAL BOILING POINT (AC t 
1 CRIMCAL TEMP(R) CRITICAL PRESSURE (N/M-MI * ) 

52 FORMAT (/F13.3,7X,FI9.3,7X,F21.3,6X,EI6.3I 

56 FORMAT (//• ENGINEERING UNITS*) 

50 FORMAT (////• MOLECULAR WEIGHT (LBM/MOLE ) GAS CONSTANT CF7-LBF, 
RATIO OF SPECIFIC HEATS*) 

(////' NORMAL MELTING PUINUF) NORMAL BOILING POINMF) 
1 CRITICAL TEMPIFl CRITICAL PRESSURE! PSI )* I 
IF (LSUPP.GT.il GO TO 55 
IF (NBRIT.EQ.il GO TO 53 
W UE 16,351 
W *Tc » 6 , 36 ) 

WRITE (6,551 

WRITE 16,37) I HDIl, 11,1-1,2), HDM, 2) 

WRITE (6,36) 

GAj«( 0.08 205 /XMW) *10* 0**5/0*9669 
WRITE (6,391 

WRITE (6,50) XMW, GAS, GAMMA 
WRITE (6,511 

WRITE (6,52) TMELT,T0OIL,TCRIT,PCR(T 
GO TO 55 

53 WRITE (6,351 
WRITE (6.361 
WRITE (6,56) 

WRITE (6,371 IHDU, 11,1-1, 21, H0I1, 21 
WRITE (6,38) 

WRITE 16,50) 

GAS-1555. 33/XMW 

WRITE 16,50) XMW, GAS, GAMMA 

WP1TC (6,51) 

IM£LT«TMELT* 9 . 0 / 5 . 0 - 559. 67 
TBUll«TBOIL* 9 . 0 / 5 . 0 - 559.67 
TCR I T-TCR 1 1 - 9 . 0 / 5 . 0 - 559.67 


/ 


Vy * 



pckit.pcrit*i4.t/io.o**s 

WRITt (6,42) IHEU, TBUIL,TCRII,PCRll 
44 RETURN 
END 


SUBROUTINE POUT 2 INK, OTEMP,OPHES,URHUl,aRHQV,U$IGN,aKlAH,UXMUl,tlXN 
IUV, GAMMA, XHN,NBR IT, OlHCL, POINT, LSUPP) 

DI PENSION aTEKP(SI,OPRbS(8l,OKHUUet,ORHuV(ei,OSIGM(m,OXLAH(B» 

01 PENS ION OXMULI8l*OXHUVm,OXNUUBI ,OXNUV( 6! «KAT 10! 8) ,OXNL Ul 
DIMENSION OHiei.SONICtei.O'.HFIBUUTHCLUl 

UI PENSION BTEPPI8I, 8P"e5I6J,BRH0U8),BRHuVI8»,BSIGM(a) ,OXLAM(0» 
DIMENSION DXMUl m,BXMUV(6>,BXNUl(Bl ,8XNuVm,BXNMB) .OH<tW 
DIMENSION BSONMO),BOSHMeUBTHCUU) 

INTEGER POINT 


66 FORMAT « * 1 • ) 

67 FORMAT!* 

I 0 E N T P 
6B FORMAT (//• 

70 FORMAT | /• 


temperature 


PROPERTIES*! 

TEMPERATURE*) 

VAPOR PRESSURE*) 


72 FORMAT I /• LIQUID OENSITY*) 
74 FORMAT I /• VAPOR OENSITY*) 


99 FORMAT I /• 
121 FORMAT « /• 


76 FORMAT I /• SURFACE TENSION*) 

7ft FORMAT I /• LATENT HEAT*) 

B'J FORMAT ( /• THERMAL CONDUCTIVITY OF LIQUIOM 
62 FO. MAT < /• OYNAMIC VISCOSITY OF LIOUIOM 
84 FORMAT C /• KINEMATIC VI SCC5IT V OF LIQUID*) 
66 FORMAT I /• DYNAMIC VISCOSITY OF VAPOR*) 

90 FURMAT ( /« KINEMATIC VISCOSITY OF VAPOR* I 
92 FURMAT I /• RATIO OF KINEMATIC V 1 SCOS I T 1 E S • ) 
9S FORMAT ( /• NICKING HEIGHT FACTOR* I 
97 FORMAT I /• LIQUID TRANSPORT FACTOR*! 

99 FORMAT I /• SONIC HEAT FLUX*) 

21 FORMAT < /• SONIC VELOCITY*) 

69 FORMAT t L7H (KELVIN) ,23X,8FLQ.l) 

71 FORMAT I 14H (N/M-M) ,?6X,8E10.2) 

73 FORMAT I1SH IKG/M-M-M1 ,2SX,8E10.2I 

7* FORMAT (LSH IKG/K-H-M) , 2SX , flG 10. 2 * 


77 FORMAT U1H (N/M| 


,29X,8EL0.2> 


79 FORMAT (1 1H (N-S/KG) ,29X,8EI0.2) 


81 FORMAT ( I9H (W/H-Kl 
63 FORMAT I16H (KC/M-S) 
8b FURMAT I17H l*-M/S) 
89 FORMAT ( 16H (KG/H-S) 
91 FORMAT U7H < H-M/S ) 


,?ix,aeio.7> 

, 24X , 8E 10 .2 ) 

,2 3X, BE 10. 2 1 
,?4X,66 10.2) 
,23X,8E10-2> 


94 FORMAT I20H ( VAPOR/l I QU 1 0 1 ,Z0X,8M0.2I 

96 FORMAT UOH (M-H) , 30X, 8E 1 0 .2 ) 

98 FORMAT I14H IN/M-M) ,26X t 8EI0.2) 

100 FORMAT (KM IW/M-M) ,2BX,BE10.7> 

119 FORMAT I18H ( H/S I ,22X,8E10.2) 

103 FORMAT I l?H (FAHRENHEIT) ,23X,8F10.W 

104 FORMAT I13H IPS!) ,27X,8E10.2) 

1 OS FORMAT II8M UBM/FT-FT-FT) ,22X,8E10.2) 

106 FORMAT I 1 AH 1 1 BM/ F T-F T-F T I ,22X,BfclU.2> 

107 FORMAT U4H (LBF/FT) , 26X, fE ’ 0. 2 1 

108 FURMAT li'iH I B TU/L8M | , 29X , DU 0.2 > 


109 FORMAT I19H ( BTU/MR-FT-F ) ,2lX,8tl0.2l 

110 FORMAT U7M UBM/FT-HR) ,23X,Sfl0.2) 

111 FORMAT I17H ( F l-FT/HR ) ,?3X,6E10.2I 

112 FORMAT I1BH I LBM/f T-HRI ,22X,8EL0.2) 

113 FURMAT I16H 1FT-FT/HRI , 24X, 8). 10.2 1 

114 FORMAT ( 12H IFT-FT) ,?8X,8M0.2) 

ns FORMAT (IBH I BTU/HR-FT-F I ) .22X,BE10.2) 

120 FURMAT I 19H I BTU/HR-FT-FT ) ,?tX,JEI0.2l 
116 FORMAT (LBH (FT/HR) , 22X, BE 10.2 ) 

IF IlSUPP.GT.l) GO TO • 1 17 
IF (POINT. EQ.l) GO TO 11 
IF INK.EQ.II GO TO 117 
11 WRITE (6,66) 

WRITE 16,67) 

DO 6 S 1 * 1 , NK 

OXNUL 1 1 I'OXMUL 1 1 l/ORHOL 1 1 ) 

cxNuvm*oxMuvm/ORHovm 
MAT 101 1 )*OXNUV( 1 ) /OXNUL ( 1 ) 

OXNL I n*OSlGMU)*OXLAH< | )«OXHQLU 1/OXMUL 1 1 1 
UHm*OSIGM(n/l9.«*ORHOLim 

SONIC! II* ! (GAHMA*01 EMPI 1 )*1.6* 1545. 33* 32.2/XMW)**0. 5| *0. 3048 
USHF ( I ) «OXL AM! I UORHGVI n*SONIC(l>/{ 2. OM GAMNAU . OM **0.5 
6S CONTINUE 

IF ( NBR IT -01 101 , 102,101 
102 WRITE ( 6,601 

WRITE (6,69) (OTEMP! Dtl-ltNK) 

WRITE (6,70) 

WRITE (6,71) ( OPRE S( 1) » I* 1 ,NK ) 

WRITE 16,72) 

WRITE (6,73) IORHOLI I 1*1*1, NX! 

WRITE (6,74) 

WRITE 16, 7S) (ORHOVd )»I-1,NK) 

WRITE 16,76) 

WRITE (6,77) (0S1GM(I),I*1,NR) 

WRITE 16,78) 

WRITE (6,79) ( OXL AMI I), 1*1, NX) 

WRITE (6,80) 

WRITE (6,81) ( QTHCL (l),I*l,NK) 

WRITE ( 6 , 02 ) 

WRITE 16,83) IQxmJLI 1),I«1,NK) 

WRITt 16,841 

WRITt ( 6 , t>S ) (OXNUL (I ),|-l,NK> 

WRITE (6,86) 

WRITE 16,89) (OXMUV(II,I-ltNK) 

WRITE (6,90) 

WRITt ( 6 , 91 ) (OXNUVm,l«l,NK) 

WRITt ( 6 , 92 ) 

WRITE (6,94) (RATIO! I ) , I»l?NK) 

WRITE (6,95) 

WRITE (6,96) (0H( I),!«1,NK) 

WRITE (6,97) 

WRITE 16,96) (QXNL(lt,l-l,NK) 

WRITt 16,99) 

WRITE (6,100) (QSHF ( 1 ) , I ■ 1, NK) 

WRITE 16 , 121 ) 

WRMt (6,119) (S0N1C(II,I-1,NK» 

GO TiJ 117 
101 on 118 1 * 1 , NK 


b2-2'0 


* 

1 

1 


drtM»M ii-uterpi 1 1 * 9 . 0 / 5 . 0 - 459., .7 

iJPKESl ) l-OPRFM I )*0. 001)1451 
BWHOLf I ItCHHULI |I*U.U6243 
dRHOV ( I |*()RHl)V t I )*0. 06243 
BTHCL ( I )*OtHCL( I 1*0. '*780 
»JS!GH( I ) sUS I GM{ | ) *0.005 709* 12*0 
dXIAMt | )-OXHM| 11*0.0004 JU2 
liXHULU l*UXMlrl I 11*2. 205* 3600.0/ 3. 281 
3XKUVU)*(IXMUVm*?.;05* 1600. 0/3. 281 
rtXNUL C 1 ) -PXNIJt I ||* 1600. 0*3. 2dl**2 
flXNUVI i )-l)XNUV( 1 1*3600. U*3.?81**2 
rtXNl ( I l-OXNL I 1 )*3.* 129/ 3.281**? 
dH( 1 1 *UH( 11*3. 281**2 
tISONI ( I 1* SON ICC I)*3.2dl* 1600.0 
118 BUSHF ( I l-OSHF | I )*3.4l 29/ 3.28 1* *2 
WRIU (6,601 . 

WRITE 16,103) t BTEMP1 I), I * 1 ,N* ) 

WRITE 16,70) 

WRITE (6,1 06 1 (6PKESI l),l-l,NKI 
WRITE (6,721 

WRITE (6,105) IBRHOK ] ), 1*1, HR) 

WRITE 16# 76 ) 

WRITE (6,106) I BRHOV (I), 1*1, NX | 

WRITE (6,76) 

WRITE (6,107) (BSICM(I) t !-l,NRI 
WRITE (6,76) 

WRITE 16,108) (BXLAH(l),|*t,NK) 

WRITE (6,30) 

WRITE (6, 10W) (BTHCL! I), I-l.HRI 
WRITE 16,0?) 

WRITE (6,1101 <fiXHUl(II,I«l,NKI 
WRITE (6,04) 

WRITE (6,111) (BXNULl 1), I* l# NX) 

WRITE (6,66) 

WRITE (6,1121 (eX*UV<Il v l«t,HKI 
WRITE (6,90) 

WRITE (6,113) IBXNUVl I), 1*1, NX) 

WRITE (6,92) 

WRITE (6,94) (RAT ICtI),l*l,NK) 

WRITE (6,951 

WRITE (6,114) <BH(t),l-l,NK) 

WRITE (6,97) 

WRITE (6,115) (BXNK I ),I-1,NK) 

WRITE (6,99) 

WRITE (6,120) ( BOSHF (11,1-1, NX) 

WRITE (6,1211 

WRITE (6,1161 ( BSON 1(1), 1*1, NX I 
117 RETURN 
ENO 




SUBROUTINE ROUTJ ( TM13, TMA4,NBR l T, UNITS, POINT, TEMP , T 1 , T?,LOUM 
INTEGER UM1TS,»»0|NT 
10 FORMAT 1*1*1 

43 FORMAT I//55M TEMPERATURE ( S I EXCEEDS TABLE RANGE, THF RANGE IS FR 
10M,F» .1 ,6H(XI TO, Fit.) ,4H(K). ) 

12 FORMAT I//20H THE TEMPERATURE I S, Ffl . 1 ,4H( K ) . ) 



I 


LOUT-*? 

n FORMAT W/MII THF TfcMPtKMUKl KA\GE IS FRUK,F8.U6H(K| TO, Ffl. 1, 4H 

iui.i 

4 format (//55H TEMPERATURE! SI EXCEEDS TABLF RANGE , THE' RANGE IS FK 
lUM, FtJ.l ,6H( F ) TO, F 8. 1 ,4H( E ) . ) 

?4 EUKMAT (//?0H THE TEMPERATURE I S* F 8. 1 ,4M( F 1 • ) 

15 FORMAT I//31H THE TEMPERATURE RANGE IS FRUM,FB.l,6H(F) T0,f8.1,4H 
1(F).) 

WHlTl (6,(01 
IF (MjRir.FQ.l) GO TO 3 
WRIT* (6,43) I MI 3, TMA4 
IF | POINT. fit. 1 ) GO TJ 11 
WKIT( 16,1?) TtMP 
60 TO 22 

11 WRITE (6,13) T 1, T 2 
2? IF (UN1TS.EQ.1I GU TO 62 
RE TURN 

3 IMI3-TM13*9. 0/5. 0-459. 67 
IMA4-TMA*.*9. 0/5. 0-459. 67 
1 EMP*. TEMP *9. 0/5.0-459 .67 
(1-11*9.0/5.0-459.67 
' T?-T/*9.0/5. 0-459. 67 

WHITE (6,4) TMl 3, IMA4 
IF (POlNT.NE.il GO To 14 
WRITI (6,24) TEMP 
60 TO 62 

14 WRITE (6,15) T1,T2 
62 RETURN 
ENO 


SUBROUTINE HYDRO ( LSUPP, LOUT, NBRI T, POINT, UNI TS,ORHOL,ORHCV,aXLAM,0 
1SIGH.0XMU1 tOXMUV, TEMP, OPRES»XMW» GAMMA t 
01 MENS ION XMUK 3CI,XMUV(30),HD(2.2l,XlAMOOOI,OXNUL(8),OXNUVt8> 
DIMENSION XTEMP ( 30), PRESS! 301,RH0V( 30) ,RHOL I 30) , SIGMA ( SO) 

01 MENS ION OTHCL (B I ,OXLAMI 8 ) ,QTEMP( Bl *OPRESf 8) ,ORHUL I il,ORHOV(B I 
DIMENSION OS (GM(6 ),THECL( 30) 

INTEGER POINT, UNITS 

15 FORMAT (//• NOTE INPUT TEMPERATURE RANGE HAS BEEN AOJUSTEQM 

16 FORMAT • • 1 • ) 

MM-9 

IND-0 
NBRIT-0 
GAMMA-1.41 
XMW-2.0? . 

TCRI T-«32.9 
TMEL 1-13.81 
TBOI L-20.27 

PCRIT- 107. *6. 895*10. **3 

DATA (HOI I,l),l"l«2l,HCl It 2 ! /CHHYDROGEN, 2HH2/ 

DATA (XTEMPU), I- 1,9, /14. 0,15. 0,17. 0,19. 0,21. 0,21. 0,25. 0,27. 0,29.0 
1 / 

UATA (PRESSIII, 1-1,51/ T.40EO, l. 30E+4, 3.46E+4, 7.55E*4, |.16E*5,2.2 
13E*5,3.27r«5,4.27E«5,6.l0E*5/ 

DATA (KHOLI I ), 1-1,91/ 7.44E*lt7.42E*l,T.36E*l v 7.22E*UT.02E*i,6.7 

I7F*1,6.41E*1,5.99L*1,5.53E«1/ 

UAIA (RHiJVI | I, 1-1,9)/ l.5?E-I«2.26E-L«4.90t-l, 9.451-1, 2.6TE+0,?. TOE 



/■■ ■■ 


I 



MU,4.08b*C‘,6.0 jE*0. 9.02640/ 

DATA ISIGHAIII, 1-1,91/ 2.94fc-3,2.81t-3,2.49H-3,2.14t-3,l. 77e-3,(.4 
1 JE-3, l.06F-J,7.10E-4, J.20E-4/ 

DATA (XIAHOS II, 1*1,91/ 4. 4 3E *5 , 4. 52E *5 , * . 4 St *b . 4.4 36 *5 , A- J4fc ♦ 4,4. 2 
13F»5,4.0TE«6,3.62Et5,3. JSE»S/ 

DAM (XNUI (tl, 1*1,91/ 2*)4E-5,2.l8€-5,l.7 8e-4,l.46E-5,l.2 , »F-5,l.l 

l4E-5,l.0TE-5,UQIE-5,9. 306-6/ 

OATA (XHUVt 11,1-1,91/ 7*006 •7,7.506-7,8.606-7,9.706-7, l.O^E-6, U2 

11C-6, l. 34C-6, 1.466-6, 1.606-6/ 

DATA CHECK 1 1, 1-1,91/ 1.03E-1, I.G6E- l , l . ME - U 1. 156-1 , I.20E- t,l.? 
. 166-1, 1.29E-1, 1.346-1, U38E-1/ 

NiNO-l 

16 aSUPP.GT.il GO TO 12 
RfcAO 15,30) POINT, UNITS 
JO FORMAT <(0131 
12 IF ( POINT- l ) 31,32,31 

32 IF (N6R1T.6Q.1I GO TU 64 
IF (LSUPP.GT.il GO TO 64 
READ (5,331 TEHP 

33 FORMAT (6F10.4I 

64 IF ITENP.LT. XTEMPtUI GO TO 34 
IF I TEMP.GT .XTEHPIMMI t GO TO 34 
GO TO 11 

31 READ (5,331 TMJN, TMAX.DT 
IF I N6R 1 1 .HE .11 GO TO 5 
THIN*(TM*M*459.67 1*5. 0/9.0 
TMAXMTHAXM59. 61 1*5. 0/9.0 
0T-DT65. 0/9.0 

5 IF MMIK.GT.XTEMPIMH)) GO TO 34 
IF ITHAX.LT. XIENPIUI GO TO 34 
IF (THIN. GE.XTEHP 111! GO TO 43 
INO-1 

TH1N-XTEHPIU 

43 IF (THAX.LE.XTEHPIHnn GO TO 44 
INO-1 

THAX-XTEHPIMH) 

44 TEHP-TH1N 

11 CALL P0UTltM0,Xnw,THUI,T601L,TCRlT,URIT.NBRtT,GAHHA,LSUPP| 

IF IP01NT.NE.il GO TO 20 

IF I TEHP.IT.XTFHPI 1 > ) GO TO 63 

IF lUhP.GT.XTEHPIHM)) GO TO 63 

20 IF (POINT. 60.1) GO TO 21 

IF ITHIN.G1.XMMPIHH)) GO TO 10 
IF ITHAX.LT. XltHPtlM GO TO 10 

21 IF IlND.tO.OI GO TO 14 
WRITE 16,161 

WRM6 (6,19) 

14 IF t POINT- 1 ) 45,46,45 
46 OT-O.O 
THAX-TEHP 

45 TIMP-UHP-OT 
99 00 5/ K-l,8 

khp-iihp^ot 

IF (K.NE.I) GU TO 68 
IF (TFHP.GT.IHAX) GO 10 63 
66 IF ( T fcHP.GT* FMAX I GU TO 61 
H“»K 

00 68 I-NINO.HH 


IF (xTtHPm-UHPl 66,59,60 

59 Il-I 
OIFKPIKI-TEMP 
UPME SIR l-PRFSSIII) 

ORHOLlKMRHUim I 
UHHUVIlO-BHUVfll I 
USIGMIKI-SIGMAI II) 

UXLAHIRMXLAHOIIII 
OAHUKRMXHULdll 
jXMUtf(*:l»XHUtfl III 
OTHCLCRMTHfcCLUII 

IF (POINT. SQ.U GO 10 61 
GO TU 57 

60 IF U.tO.ll GO TO 59 
1 1* l 

ICJM-l 

OUHP(M-TEHP 

UPR6SIM -EXP( (( ALOG(PRESSC 111 I- ALOG( PR6SSI (Ql 1 I ♦TEMPMXTEMP I 1 1 )*AL 
lUG(PRFSS((Olt-XTEHP<IOl«ALOG(PRCSS((lllll/(XTEHP(IIt-XT€HP(tOII| 
QRMOKR l-RHQLI 1 1 MIRHOL ( 1 1 l*RHQL( 10 1 !•( TEHP-XTEHPI HI|/tXTENP(lll- 
1XTEKPM0II 

ORHUVIK l-EXPI I (ALOGI RHOVl 1 1 1 l-ALOGI RHOVl (01 IMICNPd KTEHPI I IMAL 
10G( RHOVl 101 l-XUHPt lOMALOGI RHOVl III 1 1 1/1 X TEMPI K IT- X TEMPI 10 » 1 1 
OTkCL IKI-THECLI 1 1 MITHECLi 1 1 l-THECLt 101 Ml TEHP-XTtHPI I1II/(XTFMP(| 
1 1 1 -XT t HP ( 10 1 ) 

U01GH(M«SIGPA( MMiMCHAtHl-SlGMAI (OllMUMP-XUNPIU ll/IXURPI I 
lIl'XUHPUOtl 

OXLAHIK.I-XLAHOI 1 1 |*(XL AHOt II 1-XLAMOI IQ 11*1 TEHP-XT6MPI It)l/(XTEHP(t 
1 1 1-X 1 6 HP I Kll I 

OXHULIR MXHUL 1 1 1 Ml XMOLI 1 1 l-XMUL (101 !•( TEHP-XTEHPI 1 1 1 l/IXTEHPI III- 
iXTEHrilOll 

OXHUVIKMXHUVI ] I Ml XHUVt 1 1 l-XMUVt 10 1 }*( 1 EHP-XTEMP ( 1 1 II/IXTEMPI1H- 
IXTEHPIIOM 

IF (POIHT.EO.il GO TO 61 
GO To 5T 
68 CONTINUE 
51 CONTINUE 

61 CALL POUT 2 INK,QTENP,OPRFS,OAHOLfORMOV,OSIQH,OXLAM,OXHUUOXHUV.tAH 
IHA,XHU,NftRI T fOTHCL, POINT ,LSUPPI 

IF (PU1HT.EQ.1I GO TO 10 
IF (TEHP.LE.THAX) GO 10 99 
10 IF INHRIT.NE.OI GO TO 62 
NttKll-1 

IF (UNHS.EQ.2I GO TO i2 
GO TU 62 
34 Al-XUHPUI 
A2-XTEHPIHH) 

IF (PU1NT.EQ.1I GO TO 55 
T6HP-0.0 
GU 10 66 
55 ININ-0.0 
I MAX-0.0 

66 CALL POUT) (A1,A2,N8RIT, UN ITS, POINT, TEMP, THIN, TNAX, LOUT I 
6) IF (UNITS. NE. 21 CU TO 62 
IF (NHKIT.NE.O) GO TO 62 
vrtHlT-l 

IF (POlNT.NE.il GO TU II 
UU TU 64 


/»? «r: nns 
tun 



SUBROUTINE NITRU USU?P* 10 UT *N»R ! T*PblNT*'JNm,OH»Dl *GRHUV, JXLAM*U 
ISlGN.LXMUL.OXPUVdEHPtOPHES.XMW.CAMMAI 
01 «CNS (UN XNUK 3 U|*XMUV( JO I ,HD( 2,7 ) . XL A Ml>( JO I .f.X«UL(«l ,UXMUV(il 
DIMENSION XT IMP! JO > * PR£ SSI 30 1 ,RHOV( Ju I , RHUL * 3 U ) , S I GHA ( 30 ) 

DIMENSION OfMCU( 8 ),liXLXH( 8 )#OT€KPl 8 ),UPRt S I 8 ) ,OKtH)L ( 8 ) «ORUUV ( d ) 
OlHEHSlOn 0 $ FCM ( 8 I • THE CL ( JO) 

INTEGER POINT, UNITS 
16 FORMAT |* 1 *> 

15 format (//• wort input tempera turf range has been adjusted*) 

MM* | i 

|N 0*0 

NBRIT -0 

GAMMA* 1.606 

AMW- 28.02 

TCR | T* 125.96 

TMELT- 63.16 

raOIL*T 7.35 

PCRI T* 3 . 395 * 10,**6 

OATA |HDI 1 * II, 1 - 1*2 ),M 0 ( 1 , 2 ) / 8 HNnROGtN, 2 HN 2 / 

DATA (XTEMP(l),!*l f 13 )/ 65 . 0 , 70 . 0 , 75 . 0 , 60 . 0 , 85 . 0 , 90 . 0 , 95 . 0 , 10 u. 0,10 
15 . 0 * 110 . 0 , 115 . 0 , 120 . 0 , 125 . 0 / 

P OATA (PRESS! I )» IM* 13 )/ 1 . 70 E+ 6 * 3 .B 2 < ♦ T. 50 E * 6 , l. 32 E* 5 , 2 . 20 Ea 5 , 3.5 

r \3 l 6 E» 5 * 5 . 40 E» 5 * 7 . 63 E» 5 ,I.nEA 6 *l. 52 EA 6 k l . 99 Ea 6 , 7 . 58 E* 6 *J. 25 L* 6 / 

ro OATA (RH 0 Lm*l*l* 13 l/ 6 . 66 EA 2 * 8 . 6 SE* 2 * 5 . 26 t« 2 *B. 02 E« 7 *T. 78 E* 2 t ?.52 

Q) If J 2 , T. 26 C*Z* 6 . 96 fc* 2 * 6 . 91 E + 7 , 6 . 26 £* 7 , 5 . 82 E ♦/, 3 . 261 ♦ 7 * 6 . 3 IE * 2 / 

D'W A (RHUVI 11 * 1 * 1 , 131 / 9 . 206 - 1 , I.B 3 E+ 0 * 3 . 80 E * 0 * 5 . 8 ?F + 0 , 9 . 656 * 0 * 1.52 
lE + l* 2 . 21 E*L* 3 .rDE + l* 6 . 33 F+l* 5 « 90 E*l* 7 « 95 t:+l* 1. 1 36 + 2 * 2 . 106 * 2 / 

DATA «S 2 GMR<l|* 1 *l,m/l. 27 E- 2 * 2 .O 6 E- 2 * 9 . 3 OE- 3 ,e. 206 -)* 7 .l 66 - 3 9 6 .l 
15 E- 3 t 5 .l 3 E- 3 * 6 . 156 - 3 * 3 . 20 E- 3 * 2 . 28 E- 3 * 1 . 60 E- 3 , 6 . 80 t- 6 , 1 . 70 E- 6 / 

OATA (XL AMO ( 1 3 * 1 * 1 , 1 31 / 2 . 1 1 E+ 6 , 2 . 06 E*S, 2 . 016 + 5 , 1 . 96 E+ 5 , 1 , 806 + 3 * 1.6 
10 E+ 5 * 1 . 72 F+ 5 *i. 6 ?E+ 5 # 1 . 69 E+ 5 # 1 . 35 E+ 5 , 1 . 186 + 5 , 9.306 ♦ 6 , 6 . 606 * 6 / 

0 A 1 A (XMUl I 1 1 * l H I* k 3 l/ 2 . 79 E- 6 , 2 . 19 E- 6 , 1 . 76 f.- 6 * 1 . 6 IE- 4 * 1 . 206 - 6 * 1.06 
lE- 6 , 9 . 20 E- 5 , 8 . 60 E- 5 * 7 . 9 Ufc- 5 * 7 . 50 E- 5 * 7 . 10 fc- 5 . 6 . 80 f- 5 , 6 . 50 f- 5 / 

. DATA (XMUVf ll,l*l* 13 l/ 4 .* 0 E- 6 , 4 . 90 E- 6 , 5 . 386 - 6 , 5 .tJlf- 6 *S. 2 eE- 6 , 6.70 
16 - 6 . 7 . 1 OE- 6 , 7 . 55 E-j, 8 . 1 51 - 6 * 9 . OCE- 6 * 1 . 05 E- 5 , 1 . 216 - 3 * 1 . 616 - 5 / 

DATA (TH£CL< I).I*l, 131 / 1 . 6 GE- 1 * l . 526 - t • I . 4 +f-l * l. J 6 E- 1 , 1. 2 76 - 1 * 1.1 
196 - 1 , 1 . 10 E- 1 , 1 . 0 ' E-l* 9 . 20 E- 2 , 6 . 201 - 2 * T. 10 E- 2 , 6 . 00 E -2 » 4 . BGE- 2 / 

NINO- l 

IF ILSUPP-CT. 2 ) GO TO 12 
RE AO ( 5,301 POINT, UNITS 

30 FORMAT (10131 

12 * IPUJN 7 - 1 ) 31 , 32,1 

32 IF (NCR IT .£ 0 * 1 ) GO TO 66 
IF aSUPP.GT.il GO TO 66 
READ ( 5*331 UMP 
53 FORMAT ( 8 F 10 . 6 I 
66 IF (TEMP.LT. XTEMPMII GO To 36 
IF (TfMP.GT.XTEMP(MM)l GO TO 36 
GO TU U 

31 READ 15 * 33 ) IMJN, 7 MAX* 0 T 
IF (NCRlT.NE.lt GO TO 5 
THIN* (rMINA 459 .tr I < 5 . 0 / 9.0 



fMAX»(TMAX+469. 671*3. 0/9.0 
Ul *01 *5. 0/9.0 

5 IF (TNIN.GT.XTtHPfHM)) GO TO 36 
IF ( IPAX.LT.XTf HP ( ( I ) GO TO 36 
IF lfMIN.GE.XTL HP (1)) GO TO 63 
I IU* l 

TMl9>XIEMPUI 

63 IF I IMAX.LE.XTFMPIMMI » GO TD 66 
1 HU* L 

TMAX-XIEMPIHM) 

66 TEMP*TM(.i 

II GALL POUT 1 (HO* XMW* TMEL T.TBOIL* TCK(T*PCR(T *NSR t T* GAMMA *LSUPPI 
IF (PUINT.NE.II GO TO 20 
IF I TEMP.LT .XTLMPI ID GO TO 63 
IF 1 1 f HP. GT. X TEMP I MM) ) GO TU 61 
70 IF (POINT. EQ.tl GO TO 21 

IF MHiN.GT.XTEMPIMMH CO TU 10 
IF ITMAX.LT.XTEMPfl)) GO TO 10 
21 IF (IND.EO.O) GO TO 16 
HKITt (6,16| 

WRITE (6*15| 

16 IF IPUINI-ll 65*66*65 
66 OT-O.O 
TMAX-IEMP 
65 TEMP-TEMP-OT 
99 DO 57 R-l*8 
lEMP-ltMPfDT 
IF iK.NE.ll GO TO 88 
IF I TIMP.GT.TMAX) GO TO 63 
68 IF (UMP.GT.TMAX) GO TO 61 
NK-K 

DO 68 I *NI NO*MM 
IF (XTLMPI | l-TEMPI 68*59*60 
49 (I«| 

UItMP|R|*TEMP 
l)PKES(K)* PRESS! Ill 
ORHOL IK ) -RHOL Jill 
ORHOVIRl-RHOVim 
USIGMIKUSIGMA1ID 
UXLAHIIO-xLAMDIII) 

UXMUL I K 1«XMUL 1 1 1 1 
OXMUVIfC l*XMUV< III 
UThCL IK9-THECL (III 
IF (POINT. EO.ll GO TO 61 
GO TU 57 

60 IF II.EQ.I) GO TU 59 
Il-I 
10 * 1-1 

O.E/PIKI-TEMP 

O p RFS(K|»FXP(( (ALOG( PRESS! If ) l-ALOCl PRESS! 1011 I* TEMP* (XTEMP! I 1 9* At 
10GIPRESSn<J))*XTEMP(l01tALOG(PRESSt It II) )/( XTEMP (I! t-XTENPt 101 1) 
UKHOLdU-RUOL t 11 l+IRHOLi II l-RHOlt 10) IMTtMP-XTEMPlIt 1 t/IXTEMPt III- 
IXTFM»»I10)) 

URMOy | K ) a xp( ( | ALOGI RHOV (III l-ALOCt RHOVI 10)1 1 tTEMP* ( XTEMP 1 1 1 )*AL 
|UG( RHUVI IOM-XTEMPI I0)*AL0G( RMQVI 1 1 Ml I/I XTEMPI 1 1 1-XffMP( |0) ) I 
UTMCL(K)«lHtCL(IIH|THECL(ll)'IHECL 110)1*1 TEHP-X TEMPI 1 1 1 1/IXTEMPI I 
1U-XHMPUUU 

USIGM(K|*S(GMAim«(SIGMA(KII-StGHA( (0 1 1*( T£NP-XIEMP( II) l/CXTEMPU 


1 


Ut-XHHP cum 

0XLAMUI-X1 AMOt II)«(XLAMC( lll-XLAHUl I(J» |*( TE^P-X I (rPPI I f • |/|Xf| PPf I 

it i-xicMptnm 

DXMULIM-XMULI m*(XKUL(ll 1-XMUU IQI lMTtMP-XTEMM t IlimTTtfPI 1 1 1« 
IXUMHJIO; I 

UXNUV ( K ) *XMUV( 1 ll*(XHUV(ll l-XMUV I 101 I • ( TEMP-XIFNP( | | | 1/lXTE^Pf ( I )- 
IXTtMPUO)) 

IF ( PUlNT.EQ. 1 1 GO TO 61 
GO TO 57 
61 CONTINUE 
57 CONTINUE 

61 CALL POUT 2 INK, OTEMP.OPRES.ORHOl ,ORHPV,OS IGM.QXL AH.niMul ,OXMUV,GAN 
IMA, XMM,Nf)RI T,OTHCL ,POIN7 ,1 5UPP I 

If ( POINT ‘.EQ.U GO TO 10 
If 1 TERP.LE . TMAX I GO TO 99 
10 If IN6RIT.Kc.nl GO TO 6? 

NBR 1 I *1 

If (UNlTS.EQ.2t GO TO 12 
GO TO 62 
36 Al-XTEHPtll 
A7*XTENP(MMI 

IF (POINT.EQ.il GO TO 55 
TEMP*0.0 
GO TO 66 
55 ?H!H*O.Q 
THAXaO.O 

66 CUt POUT) 4 A1 , A2,NBR1T ,UH1TS# P0JNT,TEMP,TN1N, 1NAX,iOUT I 
6) If (UNITS.NE.2) GO TO 62 
If t NtlR J T.NE.O I GO TO 62 
N 8 RII*l 

IF (P0IKT.NE.1I GO TO )t 
GO TO 66 

62 RETURN 
END 


SUBROUTINE OXYGC I LSUPP.LOUT,* JRlT,POINf ,UNI TS.ORHOL, ORHOV#OXLAN#0 
1SIGM,0XNUL,0XHUV, TEMP, OPRES.XMK, GAMMA) 

DIMENSION XMUl(3Ot»XMUV(30l*HD<2#?>»XU*O(30!»OXMUU«}»GXMUV(«t 
DIMENSION XTEPPI 30 1 ,PRESS( )0 1 # RNUVI 30 ) , RH0L ( 301 • SIGMA ( )T I 
01 MENS ION OTHCll6>«0XLAM(0J.OTEMP<*)#OPRt*(9t#OHHOUftl#ORHOVl6l 
DIMENSION OSlGMfB) ,THECt t 301 
INTEGER POINT, UNITS 
16 FORMAT (MM 

15 FORMAT (//• NOTE INPUT TEMPERATURE RANGE HAS BEEN ADJUSTED* I 
MM* tl 
IND-O 
NBRI T *0 
GAMMA* 1.601 
XMM-31.9988 
1CRIT-156.7 
I MEET *56.8 
T«OIL*90.2 

DATA lH0(I,l>tl»l,2l,HD(l,2l/8HGXYGEN ,2H02/ 

PCRIT*. 507*10. **7 

DATMXTEMPI U# 1*1# 1 1 1/60,0,70.0,80. 0*90.0, 100. 0# 1 10.0# 120.0# DO. 0# 


ft 






1160.0,165.0*150.0/ 

DMA (PRESS! I)#l*l#It 1/6.601+2,6.601*3, 3.0GE +6, 9.60E+6, 2.5 )F+ 5, 5. 6 
lUE*5,9.90F*6,l.71F*6,?.71F*6,3.65f +6,6.60E*6/ 

DMA (kHUMI)»t*l,IW/ U2*»t M»l.?6t*),l.l<H«),Wl6fc + 3, 1.091*5, 1.0 
1 JtO,9.7ie*?,9.08E*2 f a.29E*2,?.6?E*2,6.85E*2/ 

DMA |RHOV( 1 1,1*1, ill/ l.3/L-l,1.19F-l,l.0lfc*0,6.l0f*U,l.n6E*l,2.3 
|UE*l,6.50fc»l,8.0DE+l, l.)3E*2»l.B0l+2#2.5iiM2/ 

0*TA I SIGMA 111, 1*1, 111/2. 06E-2, t. 82E-2 , l .58E -2 , U 16E -2 , 1 • 1 OE-2 # 8. 7 
m-),6.6Ot-1,6.07t-),?.l )E-J,l.28F-) t 5.m-6/ 

DATA IXLAMDll I, 1*1, 111/2.401*5, 2. 32E*5, 2. 26E*5,2.15f *5, 2.066*5, 1.9 
Ltf *5, l.75E*5,l.56E*5, 1.266*5, 1.016*5,7. QOE *6/ 

DA I A (XMUL m.IM.llJ/ 5.95F-6, 3.66E-4, 2.69E-6, 1 ■ 89E-6, 1.51E-6, 1.2 
lSE-4,1.1 JE-6, 1.056-4, 4. 9OE-5,9.6Ofc-5,9.40E-5/ 

• GAIA (XMUVUI, 1*1,111/ 5. OOE-b, 5.526-6, 6. 12f-6,6.78t-6,7.606-6,N.5 
10F-6,9.63F-6, l. 126-5, I. 316-5,1.666-*, 1.71E-5/ 

DATA I THECL lit, I * l, ill/ t.0)E-t , 1. 7?£- 1 , l .6U-1 , 1.69E-1 » 1.376-1,1.2 
15V -1, 1.161-1, 9. BOE-2. 8.106-2,7.106-2, 5. 801-2/ 

NIND-I 

IF (LSUPP.GT.il GO TO 12 
READ (5,301 POINT, UNI 1 S 
30 FORMAT (10(3) 

12 IF ( POINT- 1 ) 31,32,31 

32 IF IN8R1 T.EC.I I GO TO 66 
If (LSUPP.GT.l) GO TO 66 
READ (5,331 TEMP 

33 FORMAT (AMO. 6) 

64 IF mHP.tt.XTiHPmi GO TU 36 
If (UMP.r.f .XTEMPINHII GO 1U 16 
GU TU 11 

)1 READ (5,3)1 THIN, 1 MAX, OT 
IF tNMRIT.NE.il GU TO 5 
TM!N*ITM!N*65O.67l*5.0/9.0 
TMAX* ITMAX ♦659. 671*5.0/ 9.0 
0T*DT*5. 0/9.0 

5 IF (TMlN.GT.XTEMPTMMtl GO TO 36 
IF (TMAX.LT.XTENPMI) GU TO 36 
IF ITMIN.GE.RIEMPlin GO TO 63 
INO-l 

TMIN-XTEMPI U 

63 IF I TMAX.IE.XTF MPIMM) | GO TO 46 
1N0" 1 

TMAX*XTEHPIMM» 

66 TEMP.TMIN 

11 CALL P0UTl!HD,XNU,TMELT,TB01t,TCfUT,PCM1tN8RlT,GANNII,LSUPP) 

If (POINT.NE.il GO TO 20 
IF (TEMP.LT.XTEMPI1II GO TO 63 
IF UtMP.GT.XTFMPIMMI) GO 10 63 

20 (f (PCJINJ.EQ.il 60 TO 21 

IF UMIN.GT.XUHPIMNM GO TO 10 
IF ITMAX.ll.XttP.Pf II) GO TO 10 

21 If MND.EQ.OI GO TO 16 
*KJU (6,161 

MR I Tf 16,15} 

16 IF (PDINI-U 65,66,65 
66 OT'O.O 
TMAX.TEMP 

65 TEMP* IEMP-DT * 

99 OQ 51 Ml,# 




tFHP«ftHP*Or O 

IF (K.N6.1I OIJ Ui 83 CJ 

IF ITFMr.UT.TMAXI L;l 10 61 O 

88 IF (rtHP.GT.TMAXI CO TU 61 . ~3 fcJ 

NK-K > 

OQ 66 l-MNO.H* C 

IF (XfEMP(l)-TtMPI 68,59.60 ^ ^ 

OIENPIKI-TIMP ■ 03 ^ 

UPKESIK)»PR6SSt|l| (-V 

OK MOL I K I * RHQL (1X1 ^ 

GRHUVlKl-RHOVt III Q J 

GSIGHtXI-SIGMAMM O 

UXLAMIXI-XIAHDUII rO 

UXMUKKI-XMULI III fcj 

. uXMUVt K ) «XHUV( II) 

UTHCLIK |*f H6CL (111 
IF (P01NT.E0.II (.0 TU 61 
GO TO 57 

60 IF tt.EQ.U GO TO 59 
ll-l 

10*1-1 

OTEMPIKI-TEHP 

UPPESIK>*EXP((|ALOG|PRES5UIII-AlUGlPRESSllOII)*TEMPMXTEHPim<*AL 
ll’CtPRbSSf !OI>-XTEMPMa)*ALOG<PRE$S(IIIMI/UT6MP(in-XTEHP(IOlll 
JkH0UKI-RH0UIII»lRHaL(ni-RH0m0ll6|FEMP-XTtMP(lII|/CXTfcNPUI|. 
I XT l HP 1 10)1 

HRHQV |KI»E*PI 1 1 ALUGI RHOVl 1 1 I l-ALOCl RHOVl 10 1 ||*TE NPM XT6HPI III Ml 
LOGI KnOvl 101 l-XTENPI lOI^ALOGl RKOvI l III I I/I X TEMPI! 1 l-XUMPI IOIII 
QTHCL(K|«THECL(m*(THECLUn-THkClUO|>*UEHP-XIFHPU|ll/(XTFHPU 

III-XTEMPUOII 

OS I GH (K l-SI GMA( 1 1 ) # I $ IGHA ( || )-S I GNA { f 0 1 ) • ( TtHP-XTt MP ( I ( 1 1 / ( X T l MP ( I 
lll-XIEMMIOH 

OX L AH IK ) * XL AHOl 1 1 1 ♦ ( XL AHOI 1I)-XIAMU( 10 1 1 * ( Tt MP-XTE HP 1 11 1 1 / ( XHHPI I 
IM-XTfcHPIIOI! 

OXHUKKI-XHUMIMtlXHUU II l-XMUL (10 1 1 * ( TE V.P-XTE MP 111 II /< XTEMP I III- 
IXTEHPUQU 

OXMUVIK l-XMUV (1IIM XHUV1 1 1 l-XMUV( IQ I > ♦ I TEHP-XI EHPI II))/tXTEMP(ll)— 
lXTEHPIIOI) 

IF (POINT. EQ.l) GO TO 61 
CO in 57 
68 CONTINUE 
57 CONTINUE 

61 CALL POUT 2 (NK f OTEHP t OPKtS,OXHOL,ORHOV,OSlGH,OXlAH,OXHUUOXHUV,GAH 
IHA.XHNaNBRI T .OTHCL .POINT .LS'JPP I 

IF IPUINT.EQ.il GO TO 10 
IF HEHP.LE.TNAX) GO TO 99 
1C IF (NBRIT.NE.O) GO TO 62 
NSKIT-1 

IF (UNITS. F0.2I GO TO 12 
GO TO 62 
34 A1-XTEHP(1I 
A2-XTEM MHH1 

IF IIOlNT.EQ.ll GO TO 55 
TEMP :.p 
GO TL 

55 THfN-O.Q 
IMAX-0.0 



oo CALL HOUM (AUA2. NHR IT, UNITS, POINT, TEMP, THI „ THAX, LOUT | 

61 IF (UMTS. HE. /| GO 1(1 62 
IF (NHKir.Nfc.O) GU TU 62 
N8RIT-I 

IF IPUINT.NE.il GO TO 31 
GO TU 64 

62 ME TURN 
tNO 


SUBROUT INF WATER (ISUPP,L0UT.N8R II, POINT, UN I T S.ORHOL ,0RH0V,0KIAN # 0 
lSlGH,uXHUL,OXHUVtltHP,rPRES*XNWtGAHHA| 

DIMENSION XMUI *> I , XMU V I 30 1 , HOI 2.21, XL AMO f JO I , OXMUL ( 8 1 ,OXMUV ( B I . 
01 Hb NS ION XT EH. 4 30), PRESS! 301, RHOVl JOI,*HOLI 3UI,SIGMI 101 
DIMENSION OTHCLt8>,UXLAMm,UrEHP(8),UPRt$(8),ORHULm,ORHOV(|) 

Cl HENS ION OStGHIBl.THCCmOl 
INTEGER POINT , UNITS 
16 FOKHAT (MM 

IS FORMAT (//• NOTE INPUT TEMPERATURE RANGE HAS BEEN ‘ADJUSTED* I 

mm»W 

I ND»0 
N8RII-0 
GAHHA-I.n 
* *Hw-ie.o 
*fCKlT»647.0 
THELT-271.16 
TBQIL-373.U 

PCHII-218. 167*10 *2. 069.8 

DATA |H0II,1I,I«!,;nH0|1,2I/8HWATER , 3HH20/ 

DATA I XTEMPM I *{»l, 121/275.0, 300. 0, 325.0* 35U.O, 375. 0,400.0,425.0, 45 

10.0. 4 T5.0, 500.0, Sis. 0,550.0/ 

OAfAtr*fcSS( U, IM,U}/7.0tM,3.5EM,1.3E*4,4.2E*4,l.0e*5,2.4t45f4. 
WE«3,B.Sfc*5 c l.5F*6,2.6E»6,4.06«6«6.QE»6/ 

DATA (RHOL( Il,I*i,12l/T .QE* 3,9-956 *2 , 9. 85E*2,9.T6*Z*9. 66*2, 4.4E*2? 
l9.15t«2,6.9E*2 f 8.6E*2»B.3E*2«8.0E*2t7.9E*2/ 

DATA (RHOVl II, l-l* 121/5. 5E-3,2.5t-2,9.06-2,2.5E-l, 6.06*1,1. 3, 2.5,4 

1.6. 8.0. 1.36+ I, 2.26*1,3. 46*1/ 

OAT A(S IGHA ( 11,1*1,12) # 7. 56- 2 ,7. IE- 2, 6. 76-2, 6.206-2,5. BE-Z. 5. 356-02 
1,4. BSE-2, 4. 36-2,3.86-2, 3. ?5fc-2*2. 76-2, 2.26-2/ 

OATA ( XLAMOI I ) , 1*1 , 12 I/2.56+6, 2. 446*4,2. 376*6,2. 31E+6, 2.24E +6,2. 186 
1 *6,2. 16*6,2.016+6, I. 93E *6, 1.826*6, 1.7E *6, 1.556*6/ 

DATA ( XHUL 1 1 1 , !M« 121/ 1.66*3,8. 3E-4# 5. 2E-4, 3. 76-4, 2* IT •"■i .2. IE-4, 1. 
IBE-4, 1.5E-4, 1.3E-4, 1.2E-4, 1 .OE-4.9.0E-5/ 

OATA IXNUVC l I, 1-1,121/9.66-6, U02E-S,!. 122-5, 1.2F-5, 1.3*6-5,1.446- 
15*1.586-5, l • 746-5, 1 . 96-5,2. 1 26-5,2. 3 BE-5, 2.726-5/ 

OATAUHECK H. 1-1, 121/S. 86-1,6. 16-1, 6. 4E-1, 6. 7E-l,6.82E-W*.»5E-l, 
lb.0E-l,6. 76-1 ,6.566-1,6. 356-1,6. 12E-1, 5. 76E-1/ 

NINO-1 

IF (LSUPP.GT.tl GO TO 12 
KEAD 15,30) POINT, UNITS 

30 FORMAT 110131 * 

12 IF IPUiNT-ll 31,32,31 
32 IF (NQRIT.6U.il GO TO 64 
IF ILSUPP.GT.il UU TU 64 
RE AD (5.3M TtKP 
S3 FORMAT (8FI0.4I 



C.2-29 


C ' 

V 


64 IF (TEHP.ir.XTFHP(l)l CO TC 34 
IF t TERP.GT .XTEMPIMM) ) CO 10 34 
GO TU II 

31 KEAO ( 5 , 33 ) THIN, THU, DT 
IF (NtiRIT.NE.l) GO TO 5 
rH!N*(TKlNf 459 . 67 l* 5 . 0 / 9.0 
TMAX»(TMAX* 459 . 671 * 6 . 0 / 9.0 
0 T» 0 T* 5 . 0 / 9.0 

5 IF ( IHIN.Cr.XfFWPTHHM GO TO 34 
IF (THAX.LT .XTEHPtl)) GO TO 34 
IF (THIN. GE.XTEHP III) GO TO 43 
INO-I 

THIN*XTEMP(U 

43 IF { TMAX.LE.XTEHPCHH I ) GO TO 44 
INO* l 

THAX*XTERP(HM) 

44 TERP-THIN 

11 CALL PauTUH 0 ,XMW,THELT,T 8 aiL,TCRlT,PCRlT*NBR!T,GAHHA t LSUPPl 
IF (POINT. NE.l) GO TO 20 
IF nEHP.LT.XTEHP(l)) GO TO 63 
IF (TEMP.GT.XTEHPlHN) ) GO TO 63 

20 IF (POINT. EQ.X) GO TO 21 

[F ( THlN.GT .XTEHPlKHl I GO TO 10 
IF ‘(THAX.LT. XlEMPIll) GO TO 10 

21 IF UND.EO.O) GO TO 14 
WRITE ( 6 , 16 ) 

WRITE ( 6 , 15 ) 

14 IF (POINT -11 45 , 46,45 
46 DT* 0.0 
THAX-TEHP 

45 TEHP-TEHP-OT 
49 DO 57 K*> 1 ,6 

TEKP-TEHP+OT 
IF (K.NE.l ) GO TO 68 
IF (TERP.GT. TMAX) GO TO 63 
68 IF (TERP.GT. THAX) GO TO 61 
NK>K 

DO 60 1 -NIND,HH 

IF IXTENP( I )-TEHPl 68 , 59,60 

59 II-I 
OTEHP(KI«TEHP 
OPRESIX)-PRESS(II) 

ORHOUM-RHOUU) 

ORHOVIK )«RHOV( III 
OSIGH(K)-SIGHA( II) 

OXIAHCKI-KLAHOUII 
OXMUL ( K )*XHUL III) 

UXHUV | X )»XRUV (II) 

UTHCUKJ-THECLI II) 

IF (POINT. EQ. I) GO TO 61 
GO TO 57 

60 IF [I.EU.ll GU TO 59 

12-1 

10 - 4-1 

QTEHP(K)-TEMP 

0 PRE 5 (X)»EXPU ( ALOGI PRESS! II ))-Ak.QGl PRESS I lomMfcMPMXTfRPtl 1 )*AL 
10 GI PRESS! 10 M-XTEHPf IO)*ALOC< PRESS! II ))) )/( XTEMPU I l-XTEHPI I 0 >)> 
URHQL ( XI»RHOL (It ) ♦(RHQL( 1 1 l-AHQl ( 10 ) )*( TfcHP-XTfc HP| !I))/|XTEHP(II)- 



IXTENP(IOI) 

ORHOV(K)-tXP( (CAIOGI RHOVIII M-ALOGI RHOVI 10) ) I6TEHPM XTEMPI l 1 >*Al 
10G( R MOV (10) l-XTE HP ( I0)*A10G( RHUV( II IM )/(XTEMP(tI l-XTEHP( 10) 1 1 
GTHCl(K)-TH£Cll !!>*(THfcCL( II )-THECLl IQ H •( TEHP-XTEHP ( 1 1 ) ) /< XTCHP( I 
U)~XTEMPU0)1 

OSIGH(K)*SIGMA(III*(SIGNAm )-SIGHA( 10 > )* l TEHP-XTEMPI I II) / 1 XTfcPP ( I 
ll)-XTEMP(lO)) 

OXtAH(K)«XLAHO( II)f (XLARDI 1 ll-XLAMOI 10 1 1* ( TEHP-XTEHP (| I M/(XTEPP(l 
1D-XTENPI 10)) 

QXHUUtU*XHULm)*tXHUUn>-XMUUIO)mTEHP~XTEHpm )>/!XTERPHl)- 
IXTFKH(IOI) 

UXMUVIK)-XMUV( II l*IXHUV( II l-XMUVI 10 ) l» I TEMP-XTEHPC 1 1 ) I / ( XTEHPIM )- 
UIEHPUUI) 

IF (POINT. Eb.l) GO TO 61 

GO TO 57 
68 CONTINUE 
57 CONTINUE 

61 CALL POUT 2 (NR,OTEHP,OPRES,ORHOI.,ORHOV,0$lGH,OXLAH,OXHUL,OXNUV f GAN 
lHA,XNW,NaRIT,OTHCL, POINT, LSUPPi 

IF (POINT. EO.l) GO TO 10 
IF ITfcHP.LE.THAX) GO TO 99 
10 IF (NBRIT.NE.G) GO TO 62 

IF (UNI TS.EQ.2) GO TO 12 
GO TO 62 
34 Al-XTEHP(l) 

A2*XT tHP (HR) 

IF (POINT .EQ. 1 1 GO TO 55 
TEMP*0. C 
GO TO 66 
55 ruiN-o.o 
TMAX=O.C 

66 CALL P0UT3 I A1 »A2, NBA! T, UN i'TS, POINT ,TERP,TH|N,TRAX,LOUT ) 

63 IF (UNITS. NE.2I GO TO 62 
IF (NBRIT.NE.O) GO TO 62 
NBRIT-l 

IF (POINT. NE.l) GO TO 31 
GO TO 64 

62 RETURN 
END 



SUBROUTINE ARROW USUPP,LOUT,hBRIT,POINT,UNlTS,ORHCL,OAHOV,OXLAN,0 
lStGR,OXHUL,QXHUV v TERP,OPKE$,XHM,GAHHA) 

01 HENS ION XRUL (30 1 , XRUV ( 30 1 , HD( 2, 2 1 , XLAMDI 30 ) tUXHUL ( 8 ) ,OXHUV( 8 ) 
DIHENSIUN XTERP|30),PRESSI 301, RHUVI 30) *RHOL( 30 ) t SlGHAOO) 

01 HENS I UN OTHCL 18 ) ,0XLAN(6 ) , 0T6HP( 8),QPRES(8) ,QRHOLt 8 ) , 0RH0VI8 I 
01 RE NS I ON 0SIGH(8),TKECL( 30) 

INTEGER POINT, UNITS 
16 FORMAT (•!*) 

15 FORMAT (//« NOTE INPUT TEMPERATURE RANGE HAS BEEN ADJUSTED*) 
MH=24 

iun»o 

NBRIT'O 
GAMMA -1.31 
XMW»17.0 


\ 



I 


02-iTD 


TC-UT«4U*>. i 
l MEL T*l9*i.Q 
1ftOlL*239 .7 

PCK1T*U1. 1*1. 01325*10. 0**1 

DATA lHDtl,l), 1*1,2},HD( 1 , ? ) / 7HANMONI A, 3HNH3/ 

OAT A 4 XT EMP ( 1) , I - 1 , 24 » / 200 . 0.2 10.0, 220.0, 2 30.0. 2*0. 0, 710. 0. 260. 0, 
12IvJ.D*2 80.0,290.0, 300.0, 310.0, 320.0, 340. 0, 360.0,361.0,370.0, 175.0 
2, 3H0.0, 381.0,190.0, 395.0, 400. U , 406. 0/ 

DATA I PRESS < 1 >,1* 1,241/8.996*3, 1 . H2E *4 , 3.46E *4 ,6. 16E *4 , l .W *6, 
lU67t*6,2.144L*5,3.84t*5,5.56E*6, f.B8t*5, 1.0 If *6* 1.416*6, 1 - B4I *6, 
/3.0U*6.4.T1htC f 5-25E*-6,5.866*{,,l>.46E*6. / .OH »b,?.8UE*6,6.691r*b, 

19. J9tf6,1.00t W,l.U2fc*7/ 

• DATA (AHOLdl, I * 1 , 24) / 7 . 30E *2 , 7 . 166 ♦ 2 , 7 .056* J? . 6.<» )E *2 , 6. &U ♦ 2 , 6. 69 
it *2,6.566* 2 »6.43E*2,6.J0£*2»C.t5E*2,6.01E*2»1.b5fc*2,5.69E*2»5. 32«* 

. 22,4.B9E*2,4.77f *2, 4.65E *2 , 4.5QE* 2 ,4. 35E *2 ,4 . 1 7E*2 , J .986 *2 , 3. TIE *2 » 
33.406*2,2. 876*2/ 

DATA (KHOVdl , l m l,24l/8.90E-2, i.606- l* 3.206-1,5.406-1 »9.00C - 1 » 1 • 42 
16*0,2. 1 06*0, 3.016*0,4. 32£*0« 6.086 + 0, 8. 306*0, 1.116*1,1. 4 BE *l«2.43b* 
21,4. 106*1 ,4.60E*1 »5.20E *1 ,5.806*1 ,6.606*1, t .406* l , 8. 806* 1 , 1 .036* 2, 
31.306*2,2.306*2/ 

DATA I SIGH A ( 1 3 , l" 1 , 24 3 /4 .386-2 * 4 • 126-2, 3. 8*6- 2 , 3.62E-2, 3. 366-2 » 3.1 
l2E-2,2.8flE**2,2.63f-2,2.42F*2,2.20E-2,1.986-2,l.7SE-2,l.54t-2,l.l?t 
2-2, T.45E-3, 6.506-3,5. 686-3,4. 806-3, 3. 986-3, 3. 1 36- 3, 2. 286- 5, l.>Ofc- 3 
3, 0.70E-3, 0.006*0/ 

DATA (XLAMOd 1,1*1,241/ 1.486*6, 1.456* 6, 1.436*6, 1.406,6, 1.376*6,1. 3 
14E*6, 1.316*6,1.276*6, 1.24E + 6, 1.2 06 *6,1. 16E *6, 1.116*6,1. 076 •♦6, 9. 60E 
2*5* B.206* 5, T.806*5»7. 3T£*5,6.9Q£ *5,6.346*5, 5.806*5,5. 136*5,4. 306*5 
a,3.1fc*5,l.30t*5/ 

DATA (XMULd ) , 1*1,241/4. 3 2E-4, 3.606-4, 3.116-4, 2.886-4 ,2.686 -4, 2.56 
IE-4, 2. 506-4 , 2. 42E-4,2.346-4 ,2.246-4, 2.146-4 , 2. 026-4,1. 906-4,1.626“ 
24,1. 346—4, 1*256-4, l. 166-4, 1.086-4,9*886-5, 3* 756-1,7.506-1,6. 25L-5* 
34. 88E-5, 3.006-5/ 

OATA (XMUVUlt l*l,24|/6.78r-6, 7. 196-6,7.526* 6,7.946-6,8.356-6,8.76 
16-6, 9. 186-6, 9. 5 9E— 6, 1 .006-5, 1. 056-5, 1.1 OE-5, 1.156-5,1. 20E-5,1.3?£- 
21, 1.466-5, U50E-5. 1.556-5, 1.596-5, 1.65E-5, 1. UF-5, l . 78E-5. 1. 88E-5, 
32. 00E-5, 2. 306-5/ 

OATA UMECL (I I, 1*1, 24 1/5.316-1, 5.406-1 *5. 456-1 ,9.486-1,5.486-1,5.4 
176-1,5.456-1,5.416-1,5.356-1,5.256-1,5.126-1,4.976-1,4.816-1,4.466 
2-1, 4. 076-1, 3. 986-1,3. 896-1, 3.796-1, 3- 706-1, 3. 616-1 t 3. 516-1, 3. 326-1 
3,3.236-1,3. 146-1/ 

KIND*! 

IF I LSUrP.GT.il GO TO 12 
READ (5,301 POINT, UNITS 

30 FORMAT 1 10! 3 1 

12 IF (POINT-1) 31,32,31 

32 IF (N8RlT.EQ.il GO TO 64 
IF (LSUPP.GT.il GO TO 64 
READ 15,331 TEMP 

33 FORMAT ( 8F10.4 ) 

44 IF fTEMP.LT.XTEMPI 1)1 GO TO 34 
IF (TEMp.GT.XTfeMP(MM)l GO TO 34 
GO TO 11 

31 READ (5,331 TMIN, TMAX»DT 
IF (NBRIT.NE.il GO TO 5 
TMIN* I TMIN* 4 59. 67 1*5. 0/9.0 
r«AX«ITMAX*459. 671*5. 0/9.0 
DT*0T*5. 0/9.0 

5 IF (TMIN.GT.XTEMPCMMII GO TO 34 


if ithax.lt. A ll HP I 1 ) 1 U1 Hi 34 
if (tNiN.r.6.xn*p(m Go to 4* 

INDd 

TM|NcXT6HP( 1 1 

43 IF 1 1 MAX .1 6.XT6HPIMH)) GO TO 44 
I NO* 1 

lHAX*XT6MPtMM) 

44 T6HP*TM!N 

11 CALL PiHmtH0,XMW,TM£t?*TBQlt,TCMT,PtRlT*N8*tT,GAMMA,L$UPP> 

IF (PvINT.NE.l) GO TO 20 
IF UtMP.LT.XTtMPlUl GO TO 63 
IF (TLHP.G1.XTI HP ( MM |) GO TO 63 

20 IF l POINT .EQ. I I CO TO 21 

If (THIN.GT.XTEMPIMMM GO TO 10 
IF UHAX.LT.XTEMPUl) GO TO 10 

21 IF ( I NO. 6 0.01 CO TO 14 
HAUL (6,161 

HRITfc (6,11) 

14 IF (POINT-1) 45,46.45 
46 DT*Q.O 
1KAX-TEMP 

45 TEHP*TEMP-OT 
99 00 57 X*l,6 

TEMP-T6MP*0T 
IF (K.NE.U GO TO 88 
IF (TEMP.GT.TMAXI GO TO 63 
88 IF (TEMP.GT.TMAXI GO TO 61 
r,X*K 

00 68 I»N IND, MM 
IF (XTEMPUI-TEMP) 68,59,60 

59 U-l 
OT6KP!K)*T6HP 
(JPR6S(X1*PR6SM III 
URHOL ( X ) *RHOL (III 
UKHQV (Ml »R hOV (III 
OS I GH ( K ) * S I GMA I II ) 

OXLAM(K)*XLAMO( 111 
UXHUL(X)-XMUL( I I I 
UXHUV(K|»XMUVIII) 

UTHCL(K)*TH6CL( 1 1 1 

IF (P0INT.6Q.il GO TO 41 
GO TO 57 

60 IF ( I.EQ. 1 1 GO TO 59 
ll-l 

10-1-1 

UTEMP(Xl*T6MP 

OPRES IKI*EXP( ( ( ALOGIPRESSt 1 1 1 t-ALOGI PRESS! 10) I IPTEHPMXTEMPC I 1 )*AL 
IGGIPRESSdQII-XfEMPi 1Q|*AL0G(PRE$S( 1I)|))/(XTEMP(1I l-XTEMPl 10 1 1 ) 
ORHOL (KI-RHOU 1 1 »*(RHOL f 1 1 )*RHOL( 10) )• I TEMP-XTEMPI It ))/(XTEMPI 1 1 )- 
lXTEHPCIOU 

ORHOV ( K ) -EXP ( I ( AlCGI RHOVI 1 1 H-AlOGI RHOVI 101 1 l*TEMP* ( XTEMP (II )*AL 
LOG ( RHQV( (0) l-XTEMPl IO)*ALOG( RHOVI !f))H/(X TEMP (II)-X TEMP (10)1) 
UTHCL(K)*THECL( 1 1 )*CTH£CH II l-THECH 101)* I TEMP-XTEMPf I M )/( XTEMP( I 
in-XTtMPUO)) 

USIGM(K)*SI GHA( t T 1*1 S tGMA( 1 1 )-SIGMA( IU ) )*( TEMP-XTEMPf 1 1 1 1 /(XTEMP 1 1 
II l-XUMPUOM 

UXLAKIK I »XLAHD( 1 1 ) *IXLAH0( 1 1 )-XLAMO( 101 )♦( TEHP-XTCMPI 1 1 1 )/(XT£MP( I 
m-xnHPiiun 



' i 
:! 


j 


It -2D 


ux*uuki«xhuu t iimxmuli i i i-xmuli ioi »*(TFNP-xrfcMPt mmxfi wmi* 

IXTtMP'l ioi ) 

L'XMUVfKl-XHUVI 1 I )+(XMUV( II)-XHUV( 101 )•( UMP-XTFMPt I! l)/UU*P( lll- 

ixrfcMPfion 

IF IPUlNT.FQ.il GO T>) 61 
Gil TU 57 
68 CONTINUE 
57 CONTINUE 

61 CALL PUUT2 (NX,0TEHP,0PaeS,QRH0L,URH0V*0S!6M,0XLAH t UKMUL,UXM 

1UV, GAMMA f XMM, NBR I T ,QTHCL , PO INT, L SUPP ) 

IF IP01NT.EQ.II GO 10 10 
IF (TEMP.LE.THAX) GO TO 99 
10 IF (N8P1T.NE.0I GO TO 62 
N8RlT*l 

IF (UN! T£.EQ.2 I GO TO 12 
GO TO 62 
36 Al-XTEMPIU 
A2*XTEHP( MMJ 

IF lPCINT.EU.il GO TO 55 
TEMP-0.0 
GO TO 66 
55 TNIN-0.0 
TMAX-0.0 

66 CALL P0UT3 ( Al, A2.N8R I T, UNITS, POINT, TEMP, THIN, TNAX.LOUT I 
63 IF (UN ITS.NE.2) GO TO 62 
IF I NBRIT.NE.OI GO TO 62 
NBRITM 

IF (POINT.NE.il GO TO 31 
GO TO 64 

62 RETURN 
END 


SUBROUTINE HETtfA (LSUPP,L0UT,N8RI T, f OINT,UNm,ORHOl,QRHOV,OXLAH,0 
lSIGH,OXHUL,OXHUV,TEKP,OPRES*XHIi,GANKA) 

DIMENSION XMUL(30l*XMUVt30l t HD(2,2l,XLAM0(301,0XMULI8l,0XMUV(8> 

PI HENS ION XTJEMP(30I*PRESS( 30|,RH0V( 30)«RH0L( 30) » SIGMA (30 1 
DIMENSION OTHCL ( 0 I ,OXL AM { 8 1 ,QT£MP( 8 1 ,OPRES( 8} *ORH0L (8 1 tORHOVI 8 ) 
DIMENSION OS IGM( 8 I « THECL 1301 
INTEGER POINT ,UNITS 
16 FORMAt |*|«l 

15 FORMAT f //• NOTE INPUT TEMPERATURE RANGE HAS BEEN ADJUSTED*! 
NM-21 
IND-0 
N8RIT-0 
GAMMA-1.25 
XMW* 32.04 

TCR IT-513. 1 
TMFLT-175.3 
TBQI 1*337.8 

PC«1T-1156.9*6.B95*10.**3 

DATA |HD1I,IUI-1*2I#HD(1,2)/8MMETHAN01,5HCM30H/ 

DATA (XTEMPI 1 1 , 1*1 , 2 1 1/250.0,260.0, 2 70.0, 280.0, 290.0* 300.0 , 320. 0, 

1340.0. 360.0.380.0.400.0.420.0. 440. 0.450. 0.460. 0.470. 0.480. 0.485.0.' 

2490.0. 495.0.500.0/ 

OATA (PRESS (11,1*1,21 1/7.706+2, 1.576 *3, 2. 926 *3, 5. 156+ J,H. 7 IE+ 3,1.4 
l7P«4,4.25E + 4, l.0 3E+5 f ?.23E + 5,4.566+5,7.581+5,1.24f +6,2.02E+6,2.5If 





2*6, 3.09E *6,3. 7 ?E +6,4. 56E+6, 5.0IE+6, 5.49E + 6, 5.986 + 6,6.5 3F+6/ 

DAT A (kHOL ( 11,1*1*211/8. 336+2, tt. 2 56+2, 8. lttC+2, 8. 106+2,6.026 *2, 7.94 
Ik *2, 7. 75E *2. 7.546+2,7.306+2, 7. 06t +2* d. 796+2,6. 506+2, 6.1 7E+2 ,6. 00t+ 
2/,5.'i'}fc«2,5.5 76+2,5.326+2,6.m + 2,5.u0F+2,4.60fc+2,4.56E + 2/ 

OATAI <HOV( I I , l - 1,21 1/ 1 .24E-2, 2 • 10E-2, 3.90E-2, 7.10E-2, 1-26E-1.2.6G 
lfc-l,b.00F-Ifl.42E +0*2 .806+0, 6, 106+0, H. 806+0, 1.4 J6+ l #2. 7Q6+1 ,2 *98E+ 
21,3. 746+1,4.656+1, 5.826+1, 6.60E+ l , 7. 58E+ 1 ,<j.90E +1,1. 07E+2/ 

OATA ( SIGMA (I), 1-1,21 |/2. 70fc-2,2.60fc-2,2.50t-2, 2. 39E-2, 2. 29E-2.2-1 
lVE-2,l.9Rfc-2,i.78E-2, 1.57E-2, 1 .376-2* 1-166-2,9.606-3, 7. 50E-3.6.50E 
2-3, 5.50E- 3, 4. 50E-3, 3 .40E-3, 2 .90E-3, 2.40E-3, 1 .9QE-3, 1*406-3/ 

OATA (XLAMO ((),(-!, 211/1. 21E*6, 1.206+6, 1 .20E+6, 1. 19E+6, 1 . 18E+6, 1. 1 
17C+6, l.i«>F +6,1. 1QE +6, 1.066+6,1.006+6, 9. 35E + 5, 8.626 + 5, 7. 806 + 5, 7. 326 
2+5»6.80E+5,6. 2<:£*5,6.06E+5, 5.l7£+5,4.76E+5#4.?8E+5»3.68E+5/ 

OATA ( XMUL ( I), 1*1, 211/1. I9E-3, 9.90E-4, 8.40E-4, 7.15E-4, 6.10E-4,5. 30 
IE-4 ,6* 106-4 , 3. 4 OE- 4,2.526-4, 2.06E-4, 1.68E-4, 1.37E-4, 1. 11E-4,1.01E- 
*4, 9.00E-5, 7 .901 — 5* 6. 85E-5, 6.456-5,6-006-5, 5. 606-5,5.056-5/ 

OATA l XKUVUI ,1-1, 21 >/8. 096-6, 8. 406-6,8.716-6, 9, 036-6, 9. 37F-6,9. 69 
it-6, 1.03F-5, 1. 10E-5, X.16E-5, i.236-5, 1.29E-5,1.36fc-5*l-43E-5#1.47E- 
25,1. 51E-5, 1.556-5,1. 596-5, 1.62E-5, 1.646-5,1. 671-5, 1.69E-5/ 

DATA (THECL1 1 1 , 1-1,21 1/2.076- 1,2.066-1 ,2. 05E- 1, 2»05E-l,2.04E- 1 ,2-0 
13E-1 ,2.02E- 1 ,2.00E-1, 1.99E-1, 1 .98E- t , 1.966-1, i.956-1, 1.936-1, 1.93E 
2- 1,1. 926-1,1. 91E-l,1.906-i,1.90E-l,l.90E-l, I. 896-1, 1.89E-1/ 

NIND-1 

IF (LSUPP.Gl.il GO TO 12 
READ 15,?0» POINT, UNITS 

30 FORMAT l 10131 

12 IF ( POINT- 1 1 31,32,31 

32 IF (NBRIT.EQ.ll GO TO 64 
IF IISUPP.GT.U GO TU 64 
RE AO (5,331 IEMP 

33 FORMAT (BF10.4) 

64 IF (TEMP. IT.xTEHPt 1)1 GO TO 34 
IF ITEMP.GT.XTEMP(MMI) GO TO 34 
GO TO 11 

31 LEAD 15,33) TMIN,TMAX,OT 
IF (NBRiT.NE.il GO TO 5 
TN!N*(TM|N+ 4 59. 671+5.0/9.0 
TMAX=l T M aX +4 5 9. 67 1 *5. 0/9.0 
0T*DT*5. 0/9.0 

5 IF ( T HI N.GT.X Tl.* 1 ,* ( MM 1 1 GO TO 34 
IF (TMAX.LT.XTEMPUII GO TO 34 
IF (TMJN.Gfc.XTfcHPUn GO TO 43 
(ND-1 

TMIN*XTEMP(1I 

43 IF HHAX.LE.XTEMPIMMII GO 10 44 
INO-L 

TMAX-XTEMP(MM) 

44 TEMP* THIN 

II CALL POUT I (HD, XMW, TMELT ,TB0IL, TCR1T,PCRIT,N8R I T, GAMMA,LSUPP! 

IF (PUINT.N6.il GO TO 20 
IF ITfcMP.LT.XUMPllM GO TO 63 
IF ( TfcMP.GT »X TtMlPl MM ) ) GO TO 63 

20 IF IP01NT.E0.1) GO TO 21 

IF (TMIN.UT.XTEMPIMMI) GO TO 10 
IF UMAX.LT.XTEMPim GO TO 10 

21 IF UNO. 60.0) t.o TO 14 

Wit 1 T 1 16,16) 

mKITc (6,151 








Z£-?j 


14 IF (PlUNl-t) 45,46,45 
46 uno.o 
IMAX-TFMP 
45 I fcMP*T fcMP-OT 
99 UO 5? K- 1,8 
TEMP* l £ HP*DT 
IF (K.NE.l) CO TO 88 
IF f TEMP.G* .THAX) GO TO 63 
88 IF UEHP.GT.THAX) GO TO 61 
NK*K 

DO 68 I -Nl ND,HM 

IF I XTEMP 1 1 l-TEMP I 68,59,60 

59 1I-I 
OTEMPIK1-TEMP 
UPRE SIR) -PRESS! 111 
ORHOL(KI*RHQL( III 
ORHOV I K )*RHOV ( III 
USIGM(K)-SIGHA( II) 

UXLAMUI-XLAMOI 11) 

JXMULIK l-XMUL (til 
OXHUVUI-XMUV! III 
OTHtL(X)«THECl 111) 

IF (POINT. EQ. I) GO TO 61 
GO TO 57 

60 IF (I.EQ.l) GO TO 59 
11*1 

ib-i-i 

OTEMPIKI-TEMP 

OPRESCKl-EXPt UALOGI PRESS lim-HOGC PRESS! 10 )M*TEMP* I XTEMPUI )*AL 
1UGIPRESS! 10) )-XTEMP 1 101*** OG! PRESS! 1 1 ) ) ))/! XTEMP ! 1 1 )~XIEMP( IOM ) 
URHOUK)-RHOLUn*CRHOHI!>-^;:OU lOn^tTFMP-XIEHPMl n/(XTEMP( 11)' 
IXTEMPUO) I 

ORHOV IM-EXPU | AGOG! RHOVt 1 1 ) l-ALOGt RHUVlIO) ) »*TfeMP*( XTE*P(1 I >*At 
IOC I RHOVI 10 > »-X TEMPI IO)*ALOGI RHOVI I Ml )) /! XTEHPtl II-XTE MP 1 10) ) ) 
OTHCLIKI-THECLI 1 I )*( THECL! It)- TK6CL1 10 ) l«! TEMP-XTEMPI 1 1 ) >/( XTEMP (t 
m-xTEMPtion 

0SIGMIK1-S1GMAI II niSIGMAI Ill-S I GMA lion *| TEMP-XTEMPI ID l/IXTEMPI I 
lI)-XrEMP| 10) ) 

OXtAMIKI*XLAMDI in^TXLAMOl I I J-XLAMOI IOM* I TEMP-XTEMP 1 1 n 1/ I XTE MPI I 
1 1 l-XTEMP I 1U ) I 

OXMULIKI-XMUn IDAIXMOLI in-XMULIIOn*<TEMP-XTEMPIin)/IXTEMPIin- 
1XTEMP 1 10) I 

QXMUV|K)*XMUVnn*|XHUV( MI-XMUVI 10) I* ( TEMP-XTEHPI II ) I/IXTEMPIII 1- 
LXTEKPUOI) 

IF (POINT. EQ.l) GO TO 61 
GO TO 57 
68 CONTINUE 
57 CONTINUE 

61 CALL POUT? I NKfOTEMPfOPRES* ORHOL , ORHOV, OS 1GM,0XL AM,OXMUL,OXMUV, GAM 
1MA,XMU, NQR I TfOTHCL, POINT ,LSOPP) 

IF JPOINT.EQ.il GO TO 10 
IF ITEMP.LE.IMAXI GO TO 9*7 
10 IF INBRIT.NE.O) GO TO 62 
NBR1T-1 

IF IUNITS.EQ.2I GO TO 12 
GO TO 62 
34 Al-XTEMPIl) 

A2-XTEMPIMM) 


IF IPuINT.EO.II GO TO 55 
It- HP-0.0 
uO Tu 66 
55 THIN^O.O 
TMAX-0.0 

66 C/LL P0UT3 Ul, A2.N0RIT, UNITS, POINT, TeMP,TMIN,TMAX,LGUTI 
63 IF (UNITS. ME. 21 GO TO 62 
IF lNnRir.NE.OI GU TO 62 
<8*11=1 

IF (POINT. NE.l) GO TO 31 

GO TU 64 
62 RETURN 
END 



SUBROUTINE ACE TO I L SUPP , LOUT ,NBRIT, POINT, UNI TS,ORHOL, ORHOV, 0XIAM,0 
1 S ■ GM , OXMUL ,OXMUV, TEMP, OPR ES,XMH, GAMMA I . 

DIMENSION XMUL (301, XMUV ( 30) ,HD!2, 2) , XLAMDI30) ,QXKUL ( 8) •OXNUYt 8) 
DIMENSION XTFHPI 30), PRESS! 30 I, RHOVI 30) , RHOL I 30) , SIGMA 1 30) 

UIHENSION OTHCL (8 ) *OXLAM( 8 ) tOTEMPI 8) , OPRESI 8 1 ,ORHOL 1 8 ) , ORHOV! 6) 
DIMENSION OSIGMI 8 ), THE CL (30) 

INTEGER POINT, UNITS 

15 FORMAT I//' NOTE INPUT TEMPERATURE RANGE HAS BEEN ADJUSTED* ) 

16 FORMA T I ■ 1* | 

MM* L 3 
IND-O 
NBRIT-0 
GAMMA-1.1 
XMH = 5tt. C 
TCHIT-509.5 
TMEL T* 1 76.5 
lbUIL-329.7 

PCRi f=47. 2*1. 01325*10. 0**5 

DATA IlfOl 1,1), 1*1,2) ,HD( 1,2 )/7HACET0NE,5HC3HM)/ 

DATA (XTEMP | 11,1-1,131/325.0,335.0,345.0,355.0,365.0,375.0,385.0, 
1395.0,405.0,415.0,425.0,435.0,445.0/ 

DATA (PRESS(I), 1*1 , 13 1/8. 53E*4, 1.206+5, 1.61 £*5,2. 166*5,2.926*5,3.7 
i6E*5,4.84E*5,6.12E*5,7.466*5»9.33E*5,l. 13E*6, 1.366*6,1.656*6/ 

DATA IKHOL(l), 1-1,13)/ 7.556*2, 7.416*2, 7. 29E+2, 7.1 76+2, 7.026*2, 6. 9 
IDE *2, 6. 75E*2,6.&2E*2*6.4BE *2,6. 32E*2,6. 186*2,6.006*2, S.B2E*-?/ 

DATA I RHOVI I), 1-1,13)/ 1.896*0, 2. 56E*0. 3. 366*0, 4.376*0, 5.62E*0, 7.0 
18F*0, 8. 97E*0,1.09E*l,1.35E*l, 1.606*1, 1.926*1,2. 246* 1,2.60E*1/ 

DATA I SIGMM I), 1-1, 131/1. 896-2,1. 776-2, 1.646-2.1. 526-2, 1.40E-2, 1.2 
l8E-2,l.l6E-2,i.03E-2,9.iOE-3,7.80E-3,6.60E-3,5.40E-3,4.10E-3/ 

DATA (XLAMOt I), 1*1,131/5. I IE+5, 5.006*5, 4.86E*5, 4. 696*5, 4. 60E*5, 4.4 
16E+5,4.37E*5,4.2JE*5,4.14E*5,3.95E*5,3.83E*5,3.72E*5,3.60E*5/ 

DATA I XMUL (1), 1-1,13)/ 2. 50E-4, 2.236-4, 2.09E-4, 1.94E-4, 1.80E-4,i.T 
10E-4,I.62E-4,1.52E-4, l* 47E-4, 1 -40E-4, 1 • 34E-4, 1.29E-4, 1.24E-4/ 

OATA I XMUV I I), 1*1,131/ 6. 206-6,8. 50E-6.B.70E-6, 8. 906-6, 9.10E-6, 9.4 
10E-6, 9. 60E-6, 9. 70E-6, 9.90E-6, 1 .OOE-5, 1 . 026-5, 1.04E-5, 1.066-5/ 

DATA (THECLI II, I-l ,131/1.616-1,1. 576-1,1.546-1,?. 496-1, 1.456-1, 1.4 
126-1, l.37t-I,1.33E-l, 1.306-1,1. 256-1,1. 216-1, 1.16E- 1, l. 126-1/ 

Nl ND* l, 

IF (LSUPP.GT.il GO TO 12 
READ 1 5, )0 ) POINT, UNITS 
30 FORMAT (1013) 

12 IF (POINT-1) 31,32,31 



1 
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J • 


32 IF (NBRIT.EU.l) GO TO 64 
IF (LSUPP.GT.l I GU TO 64 
Kt AO 16,331 UHP 

33 FORMAT I BF 10.4) 

64 IF I r£Hp.LT.XTbNP(|)l GO TC 34 
IF ITEHP.GI.XUHP(MHl) GO TO 34 
GO TO 11 

31 READ (6,331 THIN, TMAX,DT 
IF (N8RIT.NE.il GO TO 5 
THIN* 1 1 MI N+4 59. 671*5.0/9.0 
THAX* (THAX* 459. 671*5.0/9.0 
DT-0T*5. 0/9.0 

5 IF ITMIN.GT.XTEMP(NM) I GO TO 34 
IF (TMAX.LT.XTfcMPI »)) GO TO 34 
IF ITHIN.GE.XTEMP(l)) GO TO 43 
I NO*! 

TMIN*XTEMPI 1 I 

43 IF (TMAX.LE.XTEMP(NH) > GO TO 44 
1 N 0*1 

THAX*XTEMP(HMI 

44 TEHP*THIN 

11 CALL POUT MHO, XHN, THELT.TBOll, TCRIT,PCRIT,NBRlf .GAMMA, LSUPP) 
IF (POINT. NE.ll GO 10 20 
IF CTEHP.LT. XTEHPUM GO TO 63 
IF (TEHP.GT.XTEHP(KH) I GO TO 63 

20 IF (P0lNT.EQ.il GO TO 21 

IF (fMIN.GT.XTEMP(HM)l GO TO 10 
IF UMAX. LT. XTEHPI ID GO TO 10 

21 15 (1ND.EC.0) GO TO 14 
WRITE (6,16) 

WRITE (6,161 

14 IF IPOINT-l I 45,46,46 
46 DT~0.0 
TNAX*TEHP 

45 TEHP*TEMP-OT 
99 DO 57 X-1,8 

TEHP-TEHP+DT 
IF (X. NE.ll GO TO 88 
IF (TEHP.GT.THAXI GO TO 63 
88 IF ( TEMP.M.TMAX 1 GO TO 61 
NK«K 

00 68 |-NINO»MM 
IF (XTEHPf I 1-TEHP I 66,59,60 

59 11*1 
OTEHPir.!*TEHP 
OPK£S(X)»PRESSI II) 

URHOL ( X I *RH0L (II) 

OKMOVlXl-RHOVdll 
0S1GNIK )*S(GHAll|) 

OXLAH(K)*XLAHDI 111 
UXHULIKI-XHULI II 1 
OXHUV ( X ) * XHUV 1 1 1 1 
UTHCl(K)*THECl(ll) 

IF I POINT. EQ. 1 ) GO T 0 61 
Co TO 97 

60 IF M.FQ.ll GO TO 59 
IIM 




U1 EMPtK) *TEHP 

IJPKF. j(K|*FXP(((ALUG(PKES6( 111 1-AL0G1 PRESS I 10)1 )• TEMP* ( XTEHPI I I l*AL 
IUG( PRESS! IO))-XrEMP(IO)*ALOG(PRESSU|l)))/IXrtHPMI|-XTEKPI 10 1 1) 
DKHOL I k )»r XPI ( IALUGI RHQL ( 1 1 1 l-ALOGI RHULIIO) I l*TfcMP+ ( XTEHP ( 1 1 |*AL 
lUCl KHOLIIO) >-XUMP( IU)*ALUG( RHUL (1 1 >)) )/( XTEHPt 1 1 l-XTEMPUOl ) I 
URHOVIKjst XP II (ALUM RHUVI IIJ)-AL0G( RHOVl 10) ) 1* TEMP* ( XTEMP I 1 1 )*AL 
1UC( KHOV ( 10) l-X TEMPI It) ) *ALOG( RHOVl 1 1 ) > ) ) / ( X TEMP ( 1 1 l-XTf HP( 101)1 
0TMCLIX1-THECLI ItlMTHECU Ill-THECLI 101 1*( TEMP-XTEMPI ID 1/ 1 XTEHPI I 
111-XTtHPUO) I 

0S!GM(K)*S1GHA( It )«(SIGHAI 1I)-S1 GMA I IU )>•< TEMP-XTEMPI 1 1 >)/ (XTEMP 1 1 
11 1-XTEMPI 101) 

CJXLAM{K)-XLAMOUII*(XLAHC;iII-XLAMDI I01 IMTEHP-XTENPf U )!/( XTEHPI I 
111-XTtHPt IUU 

OXHUL 1X)*XHUL (II ) *( XMUL (II l-XHULI 101 )*( TEMP-XTEMPI I1I)/(XTEHPIII>- 
1XTEHPI 101 1 

OX KUVt X ) *XHUV ( 1 1 1 ♦ * XMUV (II I— XHUVt 10) )•( TEMP-XTEMPI II )I/(XTEMP(1I)- 
IXTEHP(IO) 1 

IF tPUlNT.EQ.il GO TO 61 
GO TO 57 
68 CONTINUE 
57 CONTINUE 

61 CALL POUT 2 INK, OTEMP,OPRES,ORHOL»ORHOV,OSIGK,OXLAH, OXHUL, OXHUV, GAM 
IHA.XMW, NBR I T.GTHCL » POINT ,L5UPP) 

IF (POINT. EO.l) GO TO 10 
IF tlEMP.LE. THAX) GO TO 99 
10 IF (NBRIT.NE.OI GO TO 62 
NBRl F*l 

IF (UNITS. EO. 21 GO TO 12 
GO TO 62 
34 A1-XTEHP1U 
A2-XTENP(MM) 

IF IPOINT.EQ.il GO TO 55 
T6MP*0.0 
GO TU 66 
55 THIN«=0.0 
IMAX-0.0 

66 CALL POUT 3 I A1,A2,NBA IT, UNITS,POINT,TEHP,TMIN, THAX, LOUT) 

63 IF (UNITS. NE. 21 GU TO 62 
IF INBRIT.NE.O) GO TO 67 
NBR I T* 1 

IF (POINT.NE.il CO TO 31 

GO TO 64 

62 RETURN 
ENO 


SUBROUTINE FRE21 ( LSUPP, LOUT, NBR IT, POINT, UNI TS,0RH0L,0RH0V,0XLAH»0 
1SIGH, OXHUL, UXHUV,TEHP,OPRES,XHW,GAKMA) 

L)I HE NS I ON XMUt. ( 30 ) « XMUV( 301 ,HD( 2, 2 1 , XL AMD ( 30), OXHUL (8), OXHUV (6) 
DIMENSION XTFHP(i0),PRESSI30) , RHOVl 30) ,RHOLI 30) ,SIGHA| 30 » 

DIMENSION OTHCL ( 8 ) ,OXLAK( 8 ) ,OTEHP( 8) ,OPRES( 8 1 ,ORHOL (8) ,0RH0VI8 1 
DIMENSION USIGM(«) , THECLC30) 

INTEGtR POINT, UNITS 
16 FORMAT 1*1*1 

15 FURHAI (//• NOTE INPUT TEMPERATURE RANGE HAS SEEN ADJUSTED* I 

MM* 1 4 



fr£-£'J 


INU*U 

>*6K | r»o 

G4MNA*&.|75 

XMW-102.93 

ICRI1*451.7 

TMEL T* 1 38.2 

TBOtL-282.12 

PC*I T-750. *6. 895*10. **3 

CAT* lMD( I, 1I,1M,?),HD( l.2>/HHMtCS-2l ,6HCHU2F/ 

DATA IX IE MP| I ) , I* I « 141/235.0, ?4Q. 0,245.0, 250. D» 2 55. 0,260. 0,26 6. 0,2 
170.0.275.0,280.0,290.0, 300.0* JlU.U, 120.0/ 

DATA (PRESS! 11,1*1,141/ t. ME *4, 1.45t*4,I.m*4,?.4 6E*4 f 3.16F *4,4.0 
13E*4,5.12E*4,6.28E*4, 7.80E*4.9.62b*4,l.40E*5 t 1.96E*5,2.6afc* , » i i.65b 
2*5/ 

OAM (RMOK I>,lM,14>/1.50E*i,1.49E*), 1.4&E ♦ 3, 1 .4?E* 3, 1 .466 * 3 . 1 .45 
lfe*3,l.44E*3,1.43e*3,1.42E*3*1.4ie*3,l.»9fc*3,1.36fc*3,1.33E*3,l.31E* 
23 / 

DATA IRHOVI I ),!«!, 14J/5. 70E-1 , 7, 5QE-1 , 9.40E-1, 1.2 1 E*Q, 1.55E *0,1 .93 
IE*0,2.42E*0,2.9QE*0, J.fcOE *0, 4. 22E *0, 6. 88£ +0. 8. 1 OE *0, 1 . 1 4fc ♦ 1 , 1 . 50fc* 
21/ 

DATA (SIGMA III, l-l,l4l/2.75E-2,2.67E-2,2.58E-2,2.5lE-2*2.4?E-?,2.3 
15E-2,2.20E-2,2.2OE-2,2.l4t-2,;>.Q7E-2,l.92E-2,l.?8b-2,l.64E-2,1.48t 
12 / 

DATA { KLAHD ( 1 1, I* 1, 14)/2.65E*5,2.63E*5,2.6lE*5,2.58E*5,2.56E*5,2.5 
14E*5,2.51E*5, 2.49E*S,2.46E'-5,2.44E*5*2.3BE*5.2.33e*5,2.27E*S,2.20E 
2 * 8 / 

DATA I XMUl I 1 1, 1 - 1, 14 1 /6. 00E-4, 5.65E- 4, 5.34E-4, 5.06E-4,4.806-4,4.56 
16-4,4.33E-4,4.l5E-4,3.96E-4, 3, 806-4, 3.52E-4, 3.26E -4 , 3. 05E-4, 2. 85E- 
24/ 

DATA (XHUVI 11,1*1,14 1/8. 40E-6, 8. 65E- 6,8.95£-6,9. 206-6, 9. 50E -6. 9. 73 
lE-6,1.00E-5,1.02E-5,U05‘-5, I, 07E *5,1. i IE -5, 1. 156-5,1.196-5, 1.23E- 
25/ 

DATA 1THECI • « I, 1-1* 141/1. 40E-1, 1.386-1, 1.35E-U1.33E- 1,1. 30E-1,1. 2 
18E-1, 1. 25 E- l, 1-23E-1, 1.206-1,1. I 8E -l, 1.13E -l, 1.086-t,1.03E -1,9.806 
2 - 2 / 

NlNO-l 

IF USUPP.CT.il GO TO 12 
RE AO (5,301 POINT, UNITS 

30 format 1 101 3) 

12 IF ( POINT -l I 31,32,31 

32 16 (N8AIT.EQ.ll GO TO 64 
IF (LSUPP.GT.1I GO TO 64 
HEAD (5,33) TEMP 

33 FORMAT (BF10.4I 

64 IF ITtMP.lT.XTEMP(l)! GO TO 34 
IF ( TEMP.GT .XTEMP ( MM) ] GO TO 34 
GO TO 11 

31 READ (5,33) T M | N, TMAX ,DT 
IF lN8RIT.NE.il GO TO 5 
TMlN*ITM|N* 4 59. 671*5.0/9.0 
TMAX* | TMAX *459 .67 )*5.Q/9.Q 
01-01*5.0/9.0 

5 IF (IMIN.GT.XTIMPIMMI I GO TO 34 
IF (TNAX.lT.XT6 HP | 1 ) I GO TO 34 
I F' (IMIN. GE. XTE*P( II ) GO TO 4) 

I NO* 1 

THIN-XTEMPC 1 t 

43 IF ITMAX.LE.XTEMP(MM) I GO TO 44. 



INtfl 

f M A X-)frfcMP(MN) 

44 I f. MP- r H 1 N 

II CALL PQUTl I Hf), Xt*b, THEl T«TROIL , TCRI IfPCA 1T,NBK| T • GAMMA, L SUP PI 
16 (PUINT.K6.il GO TO 2Q 
If I TEMP. LT. /TEMPI 1)1 GO TO 63 
If (TEMP.GT.XUMPIMHM GO TO 63 

20 IF IPUINT.eO.il GO TO 21 

16 f TMJN.ifT.X TEMPI MM| ) GO TO 10 
IF | 7M4X.Lt .XTEMP J 111 GO TO 10 

21 IF t 1 ND.EC.O I GO TO 14 

H.UU I G, 16 ) 

WklU 14,15) 

14 IF (Pi 1NT-11 45,46,45 
46 OT *0.0 
THAX » TF HP 

45 TE HP»TE MP-OT 
99 UU 57 X*1 ,8 

T E MP* T EMP4 DT 
IF (K.NE.ll GO TO 88 
IF I ItMP.GT.TMAXl GO TO 63 
08 16 I TEMP.GT. TMAX) GU TO 61 
NK*K 

00 66 I*N 1N0, MM 
IF IXTEMPM J-TEMP) 68,59,60 

59 !!•] 

OT bHP ( K ) -T6MP 1 

UPRESIXI-PRESSIII I 

UR HOL I K ) * P.HOL (III 

OkHOVIX )«RH0V(1|I 

OSIGMUI-SIGMAMl) 

UXIAM|K) S XIAMG( II I 
(JX M JL IK) *7 HUL (II) 

OXHUVUI-XMUVU! ) ■ 

UrHCLU)*lHECLUll 

IF (POINT. EQ.l) GO TO 61 i 

GO TU 67 

60 IF (I. EQ.l) GO TO 59 
IT*! 

IO»I-l 

OT(MP(Kl*TEMP 

UPRESIKI*EXP(I( ALOG I PRESS 1 1 1 II- ALOG ( PR*SS( 101 1 l*TEMP* I XTEMP 1 1 ( )*AL 
1UG( PRESS ( 10) )-XTEMP( 10)*AL0G( PRESSI I1))))/(XTEMP (III -XTEMP (10)1) 

URHOL ( K ) *RHOL ( II » ♦ ( RHOL ( M l-RHOl. UO ))•( .EMP-XTEMP 1 1 II ) /( XTEMP !1 1)- 
IXTEMPl 10) I 

OKHOV ( K I *£ XP ( I I ALOG ( RHQV1 1 1 I l-ALOGI RHOV( 10)1 I *TFMP* ( XTEMPI 2 I )*AL 
1GGI RHOV( lO))-XTEMPnai*At.OG( RHUV( ( I ))))/( XTERPUI l-XTEHPI 1011) 

OTHCL I K )* THECL ( ! ( ) * I THECLT 1 1 l-THECL 1 10 ) )*( TEMP-XTOMP ( 1 1 1 1 / ( XTE**P( I 
ID-XTfcMPMOIl 

OSIGMIKI-SIGMAIM )*(SIGMA(m-SCCMAf 10) 1*1 TEMP-XTEMP( 1 1)1/1 XTEMPI 1 
iii-xTEMpnou : 

OXLAHIK)-XLAMOm)6lXLAMf)Un-XLAMOUOII*nEMP-XTtMP{nn/(XTeMPU 
II )-XTEMPt lai ) 

(JXMULIKJ-XMU*. ( m*(XMULIin-XMUL(101 )•( TEMP-XTEMPC 1 1 ) ) / (XTEMPI 1 1 1- 
IXTfcMPIIOI) 

OXMUV(K)*XMUV(ll)6(XMUV(m-XMUVnO) )*|TtMP-XTEMP(im/(XTEMPmi- 
IXTEMPl 1 0 1 t 

16 (PCINT.EQ.il GO r 0 61 


* 


■) 

A 
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CO ro 57 

6fl CONTINUE 

5 7 CONTINUE 

61 CAU POUT 2 (NK,OTEHP,OPRES,URHOL«URHOV,OSICH v nXLAH ( OXHUL ,MXMUV*l*AM 
lMA,XMW t NHRlT,OTHtL*PUINT,LSUPP) 

IF IPOlNI.FO.l) GO 10 10 
IF (TENP.LE.TMAXl GO TO 99 
10 IF (NBR 1 T .NE.O ) GO TO 62 

nbritm 

IF (UNITS. FQ. 2) GO TO 12 
GO Tu 62 
34 Al-XTEMPIU. 

A2-XTEMP I MM ) 

IF (POINJ.EQ.l) GO TO 55 
TENP-Q.O 
GO TO 66 
55 TH1 N-0. 0 
INAX-0.0 

06 CALI P0UT3 (A1,A2,N6RIT, UNITS. POINT, TEMP, THIN, TMAX.LOUT) 

63 IF (UNITS. NE. 2) GO TO 62 

IF (N8RIT.NE.0) GO TO 62 
NDRI I* 1 

IF (POINT. HE. 1> GO TO 51 
GO TO 64 

62 RETURN 
ENO 


SUBROUTINE SODIU ( LSUPP.LOUT.NBRIT.POI NT, UN I TS,ORHOL ,CRHOV,OXLAM,0 
l SIGH, OXHUL, OXNUV, TEMP, OPRES, XHW, GANNA I 
LI HENS I ON HO(2,2),XTEMP(2),OrHCL( 8 1 ,0X1 AM ( 8 ), 0 TENP( 8 ) .OPRES ( 8 ) 
DIMENSION ORHOU0I,ORHOV(0|,O$lGM(8),OXMUU8ltOXMUV(8) 

INTEGER POINT, UNITS 
16 FORMAT (MM 

15 FORMAT (//• NOTE INPUT TEMPERATURE RANGE HAS BEEN AOJUSTED* ) 
HM-2 
INU-0 
NBR 1 T»0 
GAMMA-5.0/3.0 
XMW-22.991 
TCRI T-2600.0 
THE l T-3M.0 
TBUIL* 1 156. 2 

PC R I T - 1 8 7.2* l .01325M0**5 

DATA (H0M,1),I*1,2),HD( 1,2)/0HSODIUM ,2HNA/ 

DATA (XTEHPI I) 9 I- 1,21/371.0,2200.0/ 

MNO-l 

IF (ISUPP.GT.l) GO TO 12 
READ 15,301 POINT, UNITS 
30 FOIHAT ( 1013) 
t2 IF (POINT-11 31,32,31 

32 IT (NBR IT.EO.lt GO TO 64 
lr 'ISUPP.GT.l) 00 TO 64 
kC AD 15,33) TFMP 

33 FORMAT f BF 10.4) 

64 IF (TEMP.LT.XTEMPim GO TO 34 


IF (TEMP.GT.XTEMP(MM)) GO TO 34 

te» To 11 

U MEAD ( 5 , 33 ) THIN, TMAX ,DT 
IF (NBMT.NE.il GO TO 5 
(MIN- (TM 1 N» 459 . 67 165 . 0 / 9.0 
T MAX- | TMAX ♦ 459 . 67 )* 5 . 0 / 9.0 
Ur-DT 65 . 0 / 9.0 

5 IF ( TMIN.GT .XTENP'MM I ) GO TO 34 
IF ITHAX.LT. XTEMPIIM GO TO 34 
IF I TMIN.GF.XTE HP ( 1 )) GO 10 43 
I NO- l 

TMIN-XTEMPIll 

43 IF WMAX.LE.XTF HP (MM)) CO TO 44 
I NO- 1 

TMAX-XTEMP(MM) 

44 TEMP-THIN 

11 CALL P 0 UTUH 0 ,XMK,TMELT,T 80 a,TC*lT,PCRIT,NBRIT,GAMNA,LSUPPI 
IF (POINT. NE. I) GO 1 C 20 
IF CEMP.LT.XUMPnn GU TO 63 
IF ( TEMP.GT.XTEMP(MM) ) GO TO 63 

20 IF (POINT. EO.ll GO TO 21 

IF ( TMIN.GT.XTEMP(MH) ) GO TO 10 
IF (TMAX. IT. XTEMPIIM GO TO 10 

21 IF (INO.bO.O) GO TO 14 
WRIU ( 6 , 16 ) 

WRITE ( 6 , 1)1 

14 IF (POINT- 1 ) 45 , 46,45 
46 OT-O.O 
T MAX- T EMP 

45 TEHP-/EMP-DT 
99 UO 57 K- 1,8 

TEMP-TEMP 40 T 
IF (K.NF.l) GO TO 08 
IF ( TEMP. GT. TMAX) GO TO 63 
00 IF ( TEMP. GT. TMAX) GO TO 61 
NK-K 

CTEMP(K|-TEMP 

UPRES (X ) ■ ( EXP I 9.9831 75-109 10 .06/TEMP-68623 1.9/ TEMP**? ) ) 6 10. 066$ 
URHOL (K )-( 1.01 36 30-0.0002 3 5 044 5* TEMP -0.0000000009661 0406 TEMP662 )*l 
10.06*3 

ORHOV (K I- (EXP 1 1 . 000785 - 10129 . 14 /TIMP- 575469 . 0 /TEMP 662 > 161000 . 
TF-l.h 6 IEMP- 459 . 6 T 

OTHCL(K)-' t. 4296 - 0 . 0006 1 Z 926 TF 60 . 000000 11127 * 1 F** 2 - 0. 0000000000306 
133 MF** 3 »* 7 O .07 

US IGMIK )• ( 2 2 0 . 0 * 0 . 091 6 TFMP I * 0 . 00 | 

UXLAM(R). ( 4 170 . 049 * 0 . 202984 I*TtMP- 0 . 0004765964 »TENP** 2 )* 1000.0 
UXMUL(K)-O. 093 #EXP( 15 17 . 0/1 TEMP *!. 98591 )/ 10 . 0**4 
UXNUV ( K ) ■ 1 0.0000 705571 * 0 » 00000015 619 6 & 6 TCMP- 0 . 000000000019 ) 81 3 J*TE 
1 MP** 2 )* 0.1 

IF CPQINT.EQ.il GO TO 61 
57 CONTINUE 

61 CALL POUT? INK, OTEMP, OPRES, ORHOL,ORHOV,OSIGM,OXIAM,OXMUL,OXMUV,GAM 
lHA,Xh 4 ,N 8 RIT« 0 THl.L,P 0 INT,LSUPPt 
IF (PUINT.EQ.l) tt! TO 10 
IF ( IEMP.LE.TMAX) GO TO 99 
10 IF (NhMT.NE.Of GO TU 62 
WMU l-l 
NINO- l 


I 

l 

i 




t 


. 





IF (UNITS. EQ. 21 GU TU 12 
GO TU 62 
S * i Al*XftMP| 1 1 
A2*Xl feHPlMH) 

JF (PUiNT.EO.il GO TO 55 
TtHP*0.0 
uO TU 66 
55 IMIN*0.0 
fMAX*O.U 

66 CALL t»GUM I Al,A2,NtH I T ,UNI TS, POINT, TEMP, TWIN. INAX.LnuT) 
63 IF IUNITS.NE.2I GO TO 62 
IF (NRRIT.NE.OI GU TO 62 
.‘■BRI T-l 

IF (POINT. HE. 1) GO TO 31 
GO TU 66 
62 KETURN 
eNO 


SUBROUTINE POT AS IL SUPP . LOUT.NHR I T , POINT , UNI T S . ORHOl , ORHOV,UXL AM, 0 
1 SI GH f OXH'JL*OXNUV* TEMP ,OPRES,XHW, GAMMA ) 

DIMENSION HOI 2,2 1 »XTEMP( 2 1 ,OTHCl ( 0)»OXLAMt6) ,OTEHP| 81 »OPRE SIB I 
DIMENSION 0RH0LISI,0RM0Vm.0SlGM(8),0XMUL(Bl,0XHUYm 
INTEGER POINT, UNITS 

15 FORMAT !//• NOTE INPUT TEMPERATURE RANGE HAS BEEN ADJUSTED* ) 

16 FORMAT 1*1*1 
MM-2 

INO-O 

NBRIT-O 

GAMMA-5.0/3.0 

XMW-J9.1 

ICRI 1*2300.0 

TMELT-336.9 

TBOIL-1033.2 

PCRIT-0.0 

DATA (hDI 1,11,1*1,2) ,HDI 1,2 I /9MP0T ASSIUM, 1HK/ 

DATA (XTEMPI | I , I" 1,2)2336.9, 190U.C/ 

N I ND a 1 

IF fLSUPP.GT.il GO TO 12 
READ 15,30) POINT, UNITS 

30 FORMAT ( 101 3 > 

12 IF (POINT- 1 1 31,32,31 
32 IF INBRIT.EQ.il GO TO 66 
IF lL5UPP.GT.il GU TO 66 
READ 15,331 TEMP 
3) FORMAT (8F10.6) 

66 IF ITEMP.LT.XTEMPmi GU TO 36 
IF (TEMP. GT. XTEMPI MM 1 ) ,0 TO 36 
GO TO 11 

31 READ (5,33) THIN, TMAX, DT 
IF (NBR1T.NE.1I GO TO 5 
TMIN*ITMIN*659. 671*5.0/9.0 
TMAX* (TMAX*659. 671*5.0/9.0 
0T*DT*5. 0/9.0 

3 IF (TNIN.GT.XTEMP(HM)) GO TO 36 
IF ( TMAX.L T .XTEMP (1 ) ) GU TO 36 
IF ITHIN.GE.XTtHPm) GO TO 63 



ir.Mlal 

TH l N* XT EMP ( 1 ) 

63 IF (TMAX.LF.XTEMPCMMII GU TO 66 
1 NO* l 

IMAX-XTEMP(MH) 

66 TEMP-THIN 

11 CALL PCUIl(HO,XMN,TMELT,rrtOIL,TC«IT,PCRIT,NBRlT,GAHMA,LSUPP» 

IF (POlNT.NE.il GO TO 20 
IF ITtMP.LT.XTEMPim GO TO 63 
IF ( TEMP. GT. XTEMP (MM I ) GO TO 63 

20 IF (PLlNT.EQ.il GO TO 21 

IF ( f M Ili.GT. XTEMP (HH) ) GC TO 10 
IF UMAX.LT.XTEMPiUI GO TO 10 

21 IF ( I NO. E 0*0 1 GO TO 16 
WRITE (6,16) 

WRIU (6,15) 

16 IF (POINT-!) 65,66,63 
66 UT*0.0 
THAX-TEMP 
65 T EMP* TfcNP-DT 
99 DO 5/ K* 1,0 
TEMP-TEHP+DT 
IF (K.NE.l) GO TO 88 
IF (TEMP. GT. TMAX) GO TO 63 
88 IF (TEMP. GT. TMAX) GO TO 61 
NK*K 

OTEHP(K)»TFMP 

UPRFSI K)*l EXPC 9. 19 186 3-9030. 992/ TEMP-633038. Q/TEMP**2 1)*10.0**5 
ORHOL ( K I* ( 0. 900 35 78-0.00022665 36*TEMP-0. 00000001 2 7661 7*TEMP**2 1*10 
100.0 

0RH0V(K|*(EXP(0. 8135762-8261. 15/TEMP-626986.1/TEMP**2) 1*1000.0 
IF»1.8*TEMP-659.67 

OTHCL ( X )* ( 0.96689-0 .00067906* TF *0.000000 1 3778* TF**2-0. 000000000026 
1886*TF** 3 1*70.87 

OS IGMIM* (16590.0*10. 0086029-0. 3000028169*TFM*0.001 
UXLAHIK )* (2269.079-0. 131Q665*TEMP-0.0002003039*TEMP**2 1*1000.0 
UXMUL (K I » (— 0.0006 39 0506*2. 02 H65 2/ TEMP- 56 1 .0966/1 EMP** 2* 166680. 6/TE 
LMp**3)*0. I 

OXMUVI K )* l 0 .00 00 3 6 70096*0. 000000 1982 508* TEMP- 0.000000 00006528 3 3*TE 
IMP**2 1*0.1 

IF ( POINT.EQ. 1 1 GO TO 61 
57 CONTINUE 

61 CALL POUT 2 ( NK, OTEMP,OPRES,ORHOL* ORHOV, OSIGM,OXLAN,OXMUl,OXMUV,GAM 
1HA, XMW,N8R IT, OTHCL, POINT ,LSUPP I 
IF (POINT. EQ.U GO TO 10 
IF (TEMP.IE.TNAXI GO TO 99 
1J IF (N8RIT.NE.0) GO TO 62 
NSRIT-l 
NINO-l 

IF (UNITS. EG. 2) GO TO 12 
GO TU 62 
36 A I "XTEMP ( 1 1 
A?*XTfcMP( MM) 

IF (POINT. EQ.l) GO TO 55 
IEHP-0.0 
go rn 66 
55 TM1N-0.0 
TMAX-0.0 
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66 CALL POUT 3 I 61 1 A?, NBK I T .UNITS, POINT , f EMP, THl* , TMAX, LOUT ) 
63 IF | UN l I5.NE.21 GO TO 62 
IP INBRlT.NE.Ot GO TO 62 
N*Ml f«\ 

IF (PdINT.NE.il GO TO 31 
GO TO 66 
62 RETURN 
END 


SUBROUTINE L 1 T H I ( LSUPP .LOUT ,N BRIT, POINT, UNI TS« ORHOl .ORHOV.OXL AM, Q 
lSiGN,OXHUt,OXhUV,TENP,OPRES,XHW,GAMMA) 

01 HENS I ON HDf2,21,XT£HP|2),OTHCUet.OXlAMl8l»OTEHP(8»,OPft£$lBl 
DIMENSION ORHOl ( 8 1 .ORHOVl B 1 »0S |G*( 8 1 »OXMUL C 8 1 , QXMUVt 8 1 
INTEGER POINT. UNITS 
16 FORMAT (M* I 

15 FORMAT (//• NOTE INPUT TEMPERATURE RANGE HAS BEEN AOJUSTEDM 
MM* 2 
1NQ-U 
NER1T-0 
GAHMA-5.0/3.0 
XMW-6.939 
fC«IT-3500.G 
TKELT-453.0 
TBOI L* It ' 3*0 
PCRIT-0.0 

OATA ( HD (1, 11,1-1,2), HDt l,2)/7hl I THIUM.2HLI / 

DATA (XTEHPI 1 1,1- 1,21/453.0,2100.0/ 

NlND-l 

IF (LSL’PP.GT • 1 ) GO TO 12 
READ 15,301 POINT, UNITS 

30 FORMAT (10131 

12 IF ( POINT- i 1 31,32,31 

32 IF (NUR1T.EQ.11 GO TO 66 
IF (LSUPP.GT.il GO TO 66 
READ 15,331 TEMP 

33 FORMAT I8F10.6) 

66 IF f TEMP. lT.XTEMPI 111 CO TO 36 
IF { TEMP.GT* XTfcMPl H?i I ) GO TO 3** 

GO TO 11 

31 READ (5,331 TMIN.TMAX.OT 
IF (NORIT.NE.il GO TO 5 
ININ- (TM1N+659. 67 1*6. 0/9.0 
TMAX- (TMAX+459. 67 1*5. 0/9.0 
DT«UT*5. 0/9.0 

5 IF ( TMlN.GT.Xf E MP(MMI 1 GO TO 36 
IF (TMAX.LT.XTt HP (11) Gn TO 36 
IF (TNIN.Ge.XftMPmi GO TO 63 
IND-i 

TMiN*xrfcMPm 

63 IF ( TMAX. LE.X TEMPI MM 1 1 GO TO 66 
INO-1 

f MAX* X I trip t KM1 
66 TEMP-1MIN 

11 CALL PCt;rilHC,XMN,rMFtr,fB0IL,rCRir,PCR|T,NBKrr # GAMN6,LSUPPI 
IF (POINT. NE.l) GO TO 20 


IF irfcMp.LT.XTtMPI HI GO TU 63 
IF (ItMP.GT.XTEMPfMMl I GO TO 63 

20 IF (PulNT.EO.il GU TO 21 

IF (TKlN.GT.XTEKP(MM)l go TO 10 
IF lTH*x.lT.XTfcMPim GQ TO 10 

21 IF I I NO. EC .01 GO T*J 16 
WRITE (6,16) 

WRITE (6,151 

16 IF I POINT-. 1 1 65.66,45 
66 UT-0.0 
TMAX*UHP 

65 TcMP-TEMP-OT 
99 DO 57 K-1,8 

f EMP-TEMP+OT 
IF (K.NE.ll GO TO 93 
IF (TEMP.GT.TMAX1 GO TO 61 
88 IT ( Tt Kp.GT • TMAX 1 GO TO 61 
NK*K 

OTEMPUl-TEMP 

OPRESI K) x 1.01 356*10.0** I 6.8»31* , 7877.9/TEMP1* 10.0**5 
ORHOL (K I* (0.546398-0.0000938 l 799*TEMP+0.00000000931$74l*TEHP**2 1*1 
10**3 

ORHOV Ul- l EXPfO. 4324234-15605. 72/TEMP- I 124864. 0/TEMP**2 11*10. 0**3 
TF-TEKP*9. 0/5. 0-459. 67 

OTHCL ( K 1 * < 0 . 69998+0 .0002 7992*7F*O.OOOOOOO22565*TF**2-O.OOeOOOOO002 
i4606*TF** 31 *70.87 

dS IGMf K 1- ( 454. 494B-0. 1356226*TEMP*0«Q000016|5487*TEMP**21*0.001 
UXLAHlK ) = t 26390.f>~$. 325*TCMP+0.000625*TEMP**2 1*1000.0 
UXMUL ( M* (0. 002926367-2 .668556/TEKP* 2995.26 l/TEMP**2-530257.6/T€MP 
l-*3)*0.1 

0XMUV1 Kl-IO. 000036 73815*O.OOOOQOU67182*TENP-O.OOOOOOOOOOH35025*T 
ICMP**? 1 *0. 1 

IF (PU1NT.EQ.I1 GO TO 61 
57 CONTINUE 

61 CALL POUT 2 (NK,OTEMP,OPRES,ORHQL«ORHOV v O$IGM,OXLAM v OXMUL 9 OXMUV > GAP( 
lHA,XMW,NBRIT,0THCL,P0INT 9 LSUPPI 

IF (POINT.EQ.il GO 10 10 
IF ITEMP.LE.TMAX! GO TO 99 
10 If (NBRIT.NE.O) GO TO 62 
N8RIT-1 
NlND-l 

IF (UNITS.E0.2t GO TO 12 
GO TO 62 
36 Al-XTEMPtll 
A2-XTtMP(HMl 

IF (POINT.EQ.il GO TO 55 
RMP-0.0 
GO TO 66 
55 f MIN-0.0 
r MAX=0.0 

66 CALL POUT 3 ( AI.A2 , N8R I T, UNI T$, POINT, TLMP.TMIN, TMAX, LOUT) 

63 IF (UMTS. NE. 21 GU TO 62 

IF (NPRIT.NE.Ol GU TU 62 
.mHRIT*! 

IF (PUlNT.NE.il GO TO 31 

• . • r ■> nc 

62 Rt IIM.'V • 
t 
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SU0RUUT INF MC.RCU U SUPP.LOUI ,NR« I T,PUINT,UNI TS ,<IRHOL,t;RHr,Y,UXLAM,0 
1SIGM,UXMUL,0XMUV, TEMP.OPRES.XMW, GAMMA) 

U l HE NS I ON HOlZ,21,XTEHPC?),lHHCLU).OXLAMC01,OUMPCM,CPRE$Cr.) 

DI HENSI ON ORMULIfl) f ORHOVC8) f (jSlGHl0» ? OXhUUHl.OXMUVl8) 

INTEGER POINT «UNl T S 
16 FORMAT IM'I 

lb FORMAT I//* NOTE INPUT TEMPEkATUKE RANGE HAS BEEN ADJUST) D* I 
MM«2 
INO-O 
NftRlT-0 
GAMHA-5. 0/3-0 
XKH*Z00.61 
TCRIT-1735.0 
TMElT-234.33 
16011*630.2 
PCRI T*i .Q5* 10.0**9 

OATA IHDI I* 1). 1-1,2), H0( 1,2 1/7HMERCURY, 2HHG/ 

DATA CXTEMPCI),l-l,2)/234.33,900.0/ 

NINO* L 

IF ILSUPP.GT.il GO TO 12 ^ 

READ 15,301 POINT , UNITS O JSQ 

30 FORMAT 1 1013) & hr] 

12 IF I POINT-1 ) 31,32.31 ►=? h? 

32 IF CN0RIT.EO.il GO TO 64 2 td 

IF CLSUPP.GT.il GO TO 64 Q 

READ <5,331 TEHP g 

33 FORMAT I8F10.4) rZ g 

64 IF CTEHP.LT .XTEMPC 1)1 GO TO 34 g 

IF ( T ENP.GT. XTEHPC HN) I GO TO 34 g 

GO TO 11 . ZZ 

31 Read C5,:3» tmin.tmax.dt JT? 

IF (N0RI i-NE. 1 ) GO TO 5 g 

THIN-CTMIN+459. 671*5. 0/9.0 g 

TMAX* ( THAX+459. 67 1*5.0/ 9.0 

0T-0T*5. 0/9.0 _ 

5 IF ITMIN.GT.XTEMPINHI) GO TO 34 D 

IF CTHAX.LT .XTEHPI 111 GO TO 34 ^ 

IF (THIN.GE.XTEHPC1M GO TO 43 . 

1ND-1 g 

TKIN«XTEHP(n ^ 

43 IF (THAX.LE.XTEMPIKH ) I GO TO 44 
1ND-1 

THAX-XTEHPCHMI 

44 TEKP-TH1N 

11 CALL POUTl ( HO, XMW, THELT , T0OIL , TCRI T ,PCRI T.N6RI T ,GAMHA ,LSUPP) 

IF C POINT .NE.ll GO TO 20 
IF ITEHP.LT. XTENPIUI GO TO 63 
IF (TEHP.GT.XTEMPCHH)I CO TO 63 

20 IF (POINT. EQ.l) GO TO 21 

IF ITMIN.GT.XTEHPtHH)) GO TO 10 
(F CTHAX.LT. XTEHPCLII GO tO 10 

21 IF (INO.EO.OI GO TO 14 
WRITE 16,16) 

WRITE (6,15) 
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14 IF I Pli INT-1 ) v3, 46,45 
46 UT-0.0 
TMAX* TEMP 
4b TtMPxTbHP-DT 
99 DU 57 K*1 » B 
temp* n MP + OT 
IF (K.NE.l) GO TO 08 
IF HFMP.GT.TMAX) GO TO 63 
08 IF I 7 EMP.GT - TMAX ) GU TO 61 
:4K*K 

UTfcHP(K)*TEMP 

UPRESU)-(lO.O**ClO.3735-3308.O/TEHP-O.8*ALOGl0CTEHPI | 1*133.3 
LIQUID-METALS HANDBOOK (ATOHIC ENERGY COMMISSION) 

U*H0LU1»C 14. 3B1 76-0. 002Bbl766*T£MP+C.0000003763475*TEMP**2 >*1000. 
URHOV(K)* (EXP (3. 24 3496-4559. 02/TEHP-607443.0/TEMP**2> 1*1000. 
UTHClCX)* 10.) 4648003+50. 0368/TEMP-82000.5/TEHP**2+3Z629500.0/TEMP* 
1*3-44 3661 0000.0/TrMP**4>* IOC. 0 

OS IGHU)* (407. 6255+0. 0013279*TENP~O.C002456797*TEMP**2>*0.001 
UXLAH(K)*( -0. 024* TEHP *8 3. 1566 ) *4184.0 

UXMUL CK1*C 0*00003658 7- 3. 1 9 88 39/ TEMP *2971. 399/ T6HP**2* 35440 8. 7/TEMP 
1**3) *0.1 

OXMUV IK )*( 0.000071 4 3205+0 .000000630029* TEHP+Q.000000000337 34 75*T£M 
1 P **2 ) * 0 . 1 

IF (POINT. EQ.l) GO TO 61 
57 CONTINUE 

61 CALL P0UT2 (NK c OTEMP,OPRES,ORHOL*QRHOV,OSIGH t OXLAM«OXMU|. v OXMUV,GAM 
IMAfXMWfNBK! T,OTHCL, POINT, LSUPP) 

IF CP01NT.EQ.il GO TO 10 
IF (TEMP.LE.TMAX) GO TO 99 
10 IF CNBRIT.NE.O) GO TO 62 
NOR ! T = 1 
NlND-i 

IF (UNITS. EQ. 2) GO TO 12 
GO TO 62 
34 Al-XIEMP(l) 

A2-X TfcMP( HMI 

IF CPOlNT.EQ.il GO TO 55 
TEMP-O.O 
GO TO 66 
55 TMIN-0.0 
TMAX-0.0 

66 CALL POUT 3 t A1 , A2.NBRIT, UNITS.POINT, TEMP ,TMIN, TMAX, LOUT | 

63 IF (UNITS. NE.Z) GO TO 62 
IF CNBRIT.NE.O! GO TO 62 
N0RIT-1 

IF CP01NT.NE.il GO TO 31 
GU TO 6 V 

62 RETURN 
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C. 3 VARIABLE CONDUCTANCE HEAT PIPE ANALYSIS CODE USER'S MANUAL 


C.3.1 Introduction 

This section describes the utilization of a digital computer code for Vari- 
able Conductance Heat Pipe Analysis (VCHPA). This computer program considers 
the steady-state performance of cold-wicked reservoir gas-controlled heat pipes. The 
reservoir and inactive condenser section are assumed to be in thermal equilibrium 
with the sink temperature. Flat-front analysis is used with conduction and diffusion 
effects being assumed negligible. This code consists of the following analyses: 

• Design Analysis 
e Performance Analysis 

In the Design Analysis, storage volume and gas charge requirements are 
calculated parametrically as a function of sink temperature range and allowable vapor 
temperature range (control sensitivity). The Performance Analysis presents a para- 
metric study of performance for a system within the range of specified maximum con- 
ditions. The analysis and formulation of the equations used in the program are pre- 
sented in Section C. 3. 2. A general description of the program is presented in Section 
C, 3. 3 along with a description of the program's input and output. Nomenclature is 
listed at the end of this section. 

The flow diagrams, program listing, and sample problems are presented 
in the Appendices. A listing of FORTRAN names with enginee_ mg quantities is also 
included as an Appendix. 

C. 3. 2 Analysis 

Figure C.3.1 shows a schematic of a cold-wicked gas-controlled heat pipe 
and its assumed temperature distribution. A steady-state analysis has been performed 
to determine storage volume requirements and performance within the design range. 
The following assumptions were made in performing the analysis: 

• Flat front analysis is applicable; i.e. , thermal conduction 
and mass diffusion effects are negligible. 
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Figure C.3-1. Schematic cf a Cold 7vic«.ed Gas -Controlled Heat Pipe 

and Its Assumed Temperature Distribution 
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• . The entire condenser length is active at the maximum condition, 

• The entire condenser is blocked at the minimum condition. 

• The inactive part of the condenser and the reservoir are at the 
sink temperature T^. 

• The noncondensable gas obeys the Ideal Gas Equation of State. 


Using these assumptions, the following equations apply: 


• Conservation of Mass 


m - m + m 
g g, r g, c 


C. 3-1 


Law of Additive Partial Pres ju res (applicable to inactive part of the 
condenser and also the storage volume) 


P ~ P + p 
v *v,o 


(Inactive Condenser) 


p = p + p (Reservoir) 

v v, r g, r ' ' 


e Ideal Gas Equation of State 


C. 3-2 


(p V) = (m R T) 

o b 


C. 3-3 


The above equations yield the following relationship: 


s\ (V .v' , 

V \ T o ) T v - c 


C. 3-4 


where V is the volume of the vapor space in the inactive condenser at the minimum 
v* c 

condition. Thus at the maximum condition: 


P - P 
v v, r 


m g l R T / r 

V g r /max 


C. 3-5 


At the minimum condition: 


C.3-3 


x 


. V ' , fv3.r \ 

* v v« *•'. v vv 


The ratio of reservoir volume to condenser void (inactive volume of condenser) is ob- 
tained by equating Equation (\ 3-5 and C. 


f/ ^v, max ^v t o, max \ / r o t min \ 1 

L\ *\\ min P v, o, min / V 1 o t m;ix / J 

The ratio of the mass of the noncondensing gas to the condenser volume is obtained 
from Equation C.a-t?: 

m 1 / V \ p . - p 

iL _ / , . L_ 1 VillHn. v. min 

v* n 1 1 v' / t . 

v, e g \ v, c / o t min 

The heat transfer from the condenser to the sink can bo defined as: 


Q - C. ( T - T ) 

C V o 


c.a-9 


The minimum required conductance of the condenser Is therefore: 


*c, mux (T - T ) 

' v t max o # max 


Combining Equations C. 3-1) and C. 3-10, the heat transfer at any vapor and sink tem- 
peratures within the control range is: 


Q Cl tj (T - T ) 

c, mux ' v o 


where *? is given by: 


G V 

c . ' r 

G = 1 ' V' 

e, max v f c 


U m T 
g o 


V' (p - p ) 
V, c v r v # o 


c. a-t 


The Design Analysis consists of solving Equations C. 3-7 and C. 3-8 for 

the storage volume and gas charge requirements as a function of the sink temperature 

range and control sensitivity A T y . The calculations are performed for a specified 

nominal operating temperature AT and minimum sink condition T . with 

v,n o,mm 

maximum sink temperature as a parameter. For a given AT y the maximum and 
minimum vapor temperatures are: 


AT y 

T v, max T v,n + 2 


C. 3-13 


AT v 

T v, min T v, n 2 


C. 3-14 


The control range is decreased in accordance with a specified AT v incre- 
ment until the storage requirement exceeds a specified maximum or becomes negative 
implying that a cold-wicked reservoir cannot provide the desired control. The analy- 
sis is repeated for successively increasing sink temperature ranges until the maximum 
specified range is reached. 

In the Performance Analysis, the heat transport is calculated as a function 
of vapor and sink temperature . Performance calculations arc performed for a spe- 
cific design (based on Equation C. 3-7), for a specified acceptable control range T^ ^ ^ 

and T . , and for specified extreme sink conditions T and T . . In the 

v, nun o, max o, mm 

event that the specified control range cannot be accommodated with a cold reservoir 

system (i. e. , V /V 1 < 0), the range of vapor temperatures is increased by a i°C 

r v, c 

increment and this process is repeated until control can be obtained. Once the stor- 
age requirement 1ms been determined, the maximum required condenser conductance 
G max is calculated from Equation C.3-10 and the parametric analysis is initiated for 
the minimum sink condition. At a given sink temperature, the heat transport Q is 
calculated for increasing values of vapor temperature above the minimum specified 
value until the condenser is fully active (i. e. , V 1). The analysis is then repeated 
for successively increasing values of sink temperature until the specified maximum 
sink temperature is reached. 
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c. :t. a 


Program Oescript ion 



C\ a. a. 1 General 

A listing of the program Is presented in Appendix C. The program was 
written in KOUTUAN V and was designed to operate on the UN1VAC 11CS system. The 
Property Data Acquisition Code is required as a subroutine to determine working fluid 
vapor pressures. The KOUTUAN names and the physical quantities they represent 
are listed in Appendix H. Storage requirements are on the enter of SO, 000 words 
(octal). 

The flow diagram is included Jn Appendix A us an aid In the overall pro- 
gram logic. Input for this program is in a NAMELIST format. This allows the user 
to run multiple eases by changing only the variables which are different from those 
of the previous case. Ail other inputs are reinitialized as for the preceding case. 

The program contains options to perform cither the design analysis and/or perfor- 
mance analysis. Basically the program reads the input data, performs the required 
analysis, and outputs the data. 

The dock setup as shown in Figure O. 3-3 consists of job control cards, 
the program source dock (which includes fluid property data acquisition code), addi- 
tional control cards followed by the input data and program termination cards. 

< \3.3.3 Input Description 

The entries to bo made on the various input cards are described in Table 
O. :i-l. The NAMELIST group name, KOUTUAN name, format, and units to be used 
are indicated for each entry. A listing of sample input data of an ammonia heat pipe, 
using nitrogen as the noneondonsable gas, is presented in Table C.3-2 for two cases. 

Table C.3-2 Sample Input Data 

AvmONIA 

NIT^Oi^h 

i* CON TRO M0*t=2tLlQUID=5tLCI=2*LC2*lS 

iiJAT A IN TOMAX = 3O8.OfTOMlN*268*O*DTO=5.OtTVMAX = 328.O»TVMlN = 3O0«Of 
l»TV = l ♦ D#DOTV=UOt«G = 296.b9*QMAX=6.0*VRVC=50.0$ 

AvjmonIA 
N I 1 ORBt N 

► COMTPO MOP£ = :itLCl = l f LC2=2S 
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Table C.3-1. Input Data Description 


Input 

Card 

No. 

Format 

Fortran 

Name 

Description 

Units 

1 

2A6 

HDl 

HD2 

Headings (working fluid) 

- 

2 

2A6 

HD3 

HD4 

Headings (non -condensible gas) 

- 

3* 

- 

CONTRO 

NAMELIST group name 

- 


integers 

MORE 

Control Point, integer = 0 for last case; 
otherwise integer = 2 

- 



LC 1 

Control Paint, integer = 2 for both 
design and performance analyses; 
otherwise integer = 1 

— 



LC 2 

Control Point, integer = 2 for perfor- 
mance analysis; otherwise integer = 1 

- 

4* 

- 

DATAIN 

NAMELIST group name 

- 


Floating 

Point 

Constants 

TOMAX 

TOMIN 

Maximum sink temperature 
Minimum sink temperature 

°K 

°K 



DTO 

Increment of the sink temperature 

°K 



TVMAX 

Maximum vapor temperature 

°K 



TVMIN 

Minimum vapor temperature 

°K 



DTV 

Increment of the vapor temperature 

°K 



DDTV 

Operating temperature control range 

°K 



RG 

Gas constant of non-condensible gas 

N m/kg 



QMAX 

Maximum heat transport 

W 



VRVC 

Maximum allowable ratio of * sservoir 
volume to condenser void 

- 


* Start with a $ in column 2 followed immediately by NAMELIST group name. New 
data follow and are ended with $. 
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The first case requires both the design and performance analysis. In the second case, 
only the performance analysis Is requested. 

C. 3.3. 3 Output Desc ript ion 

The program outputs essentially all input data. In the Design Analysis, 

the reservoir-to-eondonsor volume ratio am! mass-to-eondenser volume ratio are 

output in tabular form versus operating temperature control range AT for T . < 
r « « v o, mm 

T < T .If the required temperature control cannot be satisfied by a particular* 

o “ o, max 

sink temperature, the following statement will appear in the table. "RESERVOIR ('AN 
NO LONGER PRO VIDE DESIRED TEMPER ATI 1 RE CONTROL". When the Performance 
Analysis is requested, the required reservolr-to-eondenser volume ratio correspond- 
ing to the specified input is output along with the noneondensable mass-to-condenser 
volume ratio. This is followed by tables of the heat dissipated and the ratio of active 
condenser length to total condenser length versus operating vapor temperature for 
successive sink temperatures. If the specified vajx>r temperature range had to be 
adjusted in order t o have a working cold-reservoir system, the statement "VAPOR 
TEMPERATURE RANGE UAS REEN ADJUSTED" will appear before the tabular data. 
The storage requirements corresponding to these adjusted temperatures are calcu- 
lated and used in the analysis. A listing of typical output is given in Appendix D with 
the sample problem. 


(\ 3. 3. -1 Nomenclature 

The following is a listing of the nomenclature and associated symbols used 
in the variable conduct anee heat pipe analysis. 


\ 
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NOMENCLATURE 


Symbol 

Description 

Units 

A 

Heat rejection area per unit length of condenser 

m*/m 

G 

Thermal conductance 

W/°K 

h 

Heat transfer coefficient 

W/°K m 1 

La 

Active length of condenser 

m 

m 

Mass of non- condensible gas 

kg 

P 

Pressure 

N/m* 

Q 

Hcnt transport 

w 

R 

Gas constant of non-condensible gas 

N m/kg * 

T 

Temperatui'e 

°K 

V 

Volume 

m s 

Subscript: 



c 

Condenser 


off 

Effective 


£ 

Non-condensible gas 


n 

Nominal 


max. 

Maximum condition 


min. 

Minimum condition 


o 

Sink 


i 

Reservoir 


V 

Vapor 




Appendix A. Flow Diagram 


for the Variable Conductance Heat Pipe Analysis Code 
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„ 
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COLO RESERVOIR CAJ 
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Nl-I 

OTMAX * OTMAX + DTO I 


CALL OAC POP VAPOR 
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(v,/V c > fl(m/v c ) 


otvmax + ddtv 
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Appendix B. FORTRAN Names 
and Associated Physical or Engineering Quantities 
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Fortran Name 


Description 


Units 


OUT 1 

Operating temperature control range (AT V ) 

K 

OUT 2 

Reservoir - condenser volume ratio (V r /V c ) 

- - 

OUT 3 

Ratio of the mass of the gas to the condenser 
volume (m/V c ) 

kg/m 3 

CONTROL 

A NAMELIST group name 

- 

MORE 

Control Point, integer = 0 for last set of data, 
otherwise integer = 2 

- 

LIQUID 

Control Point, type of working fluid 

- 


1 

Hydrogen 

2 

Nitrogen 

3 

Oxygen 

4 

Water 

5 

Ammonia 

6 

Methanol 

7 

Acetone 

8 

Freon-21 

9 

Sodium 

10 

Potassium 

11 

Lithium 

12 

Mercury 


LC1 

Control Point, integer = 2 for both Design Analysis 
and Performance Analysis, otherwise integer 2 


LC 2 

Control Point, integer = 2 for Performance 
Analysis only, otherwise integer ¥ 2 

- 

DATAIN 

A NAMELIST group name 

- 

TOMAX 

Maximum sink temperature 

o 

K 

TOMIN 

Minimum sink temperature 

°K 

DTO 

Increment of the sink temperature 

°K 

TVMAX 

Maximum vapor temperature 

°K 

TVMIN 

Minimum vapor temperature 

o 

K 



Fortran Name 

Description 

Units 

DTV 

Increment of the vapor temperature 

°K 

DDTV 

Operating temperature control range 

°K 

RC 

Gas constant of the non-condensible gas 

N m/kg 

QMAX 

Maximum heat transport 

W 

VRVC 

Maximum allowable ratio of reservoir volume 
to condenser vapor volume 

- 

HD1, HD2 

Headings (working fluid) 

-■ 

HD3, HD4 

Headings (non -condensible gas) 

* - • 

KOUNT 

Control Point for calculations 

- 

LSUPP 

Control Point for the Property Data Acquisition 
Code 

- 

LOUT 

Control Point for the Property Data Acquisition 
Code 

- 

TVN 

Nominal operating temperature 

°K 

IND 

Control Point for different analysis 

°K 

OTMAX 

Sink temperature 

°K 

NN 

Control-number of calculations 


MM 

Control-number of calculations 

- 

DAC 

Property Data Acquisition Code 

- 

ARHOL* 

ORHOL 

Liquid density 

kg/m 3 

ARHOV* 

ORHOV 

Vapor density 

kg/m 5 

A XL AM* 
OXLAM 

Latent heat of evaporization 

W s/kg 


* Data not used in this program but included as a general format for calling up 
the Property Data Acquisition Code, 



Fortran Name 


ASIGM* 

OSIGM 

AXMUL* 

OXMUL 

AXMUV* 

OXMUV 

P' r OMA 

XMW 

GAMMA 

PVOMI 

VTD 

VTMAX 

VTMIN 

PVMAX 


Description 

Surface tension 


Units 

N/m 


Dynamic liquid viscosity 


kg/m s 


Dynamic vapor viscosity 


kg/m s 


Vapor pressure at the maximum sink 
temperature 

Molecular weight of the working fluid 

Ratio of specific her.ts of the working fluid 

Vapor pressure at the minimum sink temperature 

Increment of vapor temperature 

Maximum vapor temperature 

Minimum vapor temperature 

Vapor pressure at the maximum vapor tempera- 
ture 


N/m* 

kg/mole 

N/m* 

°K 

°K 

°K 

N/m* 


. PVMIN 

VCDVR 

VRDVC 

NNN 

NNP 

CTVMA 

CTVMI 

GMAX 


Vapor pressure at the minimum vapor tempera- 
ture 

Condenser-reservoir volume ratio 
Reservoir-condenser volume ratio 
Control-number of calculations 
Control-number for temporary storage 
Maximum vapor temperature 
Minimum vapor temperature 
Maximum thermal conductance 


N/m* 


K 


K 


JO 

W/ K 


* Data not used in this program but included as a general format for calling up 
the Property Data Acquisition Code. 
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Fortran Name 

Description 

Un 

NA 

Control-number of calculations 

- 

NB 

Control-number of calculations 

- 

TO 

Sink temperature 

°K 

KT 

Counter for calculations 

- 

PVO 

Vapor pressure at sink temperature 

N/m a 

PV 

Vapor pressure 

N/m* 

F 

Effective length 

m 

Q 

Heat transport 

W 
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Appendix C. Program Listing 
of the Variable Conductance Heat Pipe Analysis 
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u(Nf nsiun uut nsu),uwf/{ *»oi,<jur m su 
named si/cuntmi /mokf , i Iuimo.iC i,tC -> 

SAHU I ST /OAT A IN/ I UMAX, T..N|!i,(: 1 ii t IVMX, I VMI N»(J I V ,f*OTV,KG,-.}MAX . VftvC 

n fim*ai {* 1*1 

U FJRMAJ I////* V A K { A M L ( tLVUUCfA^Cfc H I A t 
l PIPE 4 N * L t S I St VLHP* } • I 
U f-UHMAj (//* 1C I) 10 • » I C H J B U ^ K H‘ 

l I #* t * * 

W EOaMAI (////* 0 \ S 1 ft * 4 {. A t t S I SM 

;o fdkmat t //• input o a r t *t 

21 l-UKMAl l//« HAX1-UM AUUWA0U SINK TtHPU) MINIMUM SINK ftHPt 
Ul*) 

2 ) FURMM (//• NOMJNAL OPtRAIlNG Tt*PtM MAX, OPf RATING TCMPU) 

1 H1N. UPE RATING TEMPI*)*) 

25 FORMAT I//* MAXIMUM AUOwAbU RATIO UP RFSCkVOIK VOLUME 10 COHDEN 
15M VOID* t 

26 FORMAT I//* OUTPUT D A I 4 *) 

2 ? FORMAT (//• MAX. SINK TEMPI*} { OT ) VAPOR < 1C ) (VR/V 

1C) <H/VC ) (KG/H-M-h) * ) 

5J FORMAT i //• RESERVOIR CAR NO LONGER PROVIDE OESIREO TEMPERATURE C 
10NTRU* J 

16 FORMAT 1////* PERFORMANCE ANALYSIS*) 

121 FORMAT I//* MAXIMUM SINK TEMPI*) MINIMUM SINK TfMP(K)*) 

57 FORMAT I//* MAXIMUM Hf AT TRANSPORT ( W ) * J 

123 FORMAT J //■ MAX. VAPOR TEMPI*) MIN. VAPOR TEMPOS*) 

23* FORMAT (//• VAPOR TEMPERATURE RANGE HAS SEEN A0JUS7E0 *1 
30 FORMAT f//« RATIO OF RESERVOIR VOLUME TO CONOENSER VOJDM 
39 FQKMAt (//• RATIO OF THE MASS Of NON-C ONDE*S I SLf TO CONOENSER VCJ 
IDI KG/M-m-m) * i 

AO FORMAT (//* SINK TEMP IK ) VAPuR TEMP(K) HEAT TRANSPO 

UUVJ ACTIVE CONDENSER lENGTH/TClTAL CONDENSER LENGTH* | 

10 FORMA! I2A6I 

1* FORMAT ( / / // I 0H MORKING FLUIO *246) 

15 FORMAT t//2*H NQN-COUDENS I BL E GAS ,2A6) 

16 FUrtMAT (//20M CAS CONSTANT (N-H/X&-K) ,F0.3) 

22 FORMAT I/FIS.2*IRK#F1S«2) 

2* FORMAT {/F19.2,7K,f V9.2,7K,H9.2) 

2« FORMAT ( /f 20. 5) 

34 FORMAT </2OX,F20.3,2t?0.6) 

122 FORMAT \r, 1*,3,?X,F1*.3I 

50 FORMAT T/F16.*) 

42 FORMAT l/?OX,F20.3,€2O.6,2OX,FlQ,5t 

51 READ (5,101 HDJrHD? 

RE AO 13,10) HD3,H04 
RC AO (S,CONTRO ) 

CE-AO (l.OATAIM 
MOUNT «CJ 

LSUPP-2 

LOUT-0 

TVN-UVMAX-7YM1N)/2.0*TVMIN 
WRITE (6,11) 

WRITE (6,12) 

*R I Tt (6,13) 

WR 1 7 c (6,2*1 HDI.H02 
MR I Tt 16,13) HO J, HO* 

WRITE (6,161 MG 
1N0*0 

IF (LCI.NF.2) CO TO 17 


IN*'* l 

)? IF ( tC2.EC.2t Ui TO 11 
wRIIt 16, III 

*«Ut (6,14) 

<*«m (6,201 

N.'UU (6*21 ) 

•Ml ft (6,221 TOMAX,TOMIN 
MR 1 T F (6,2 ) > 

MRIIv (6,'*) T VN# T VMA X , T VH IN 

MHlTfc (6,231 

wKITi (6, SO) VR VC 

OTHAX^TOMlN-OTa 

6N-( TOMAX*JOMlN)/OIO*l 

MM* ( TVMAX-TVMIN) /DDT V* I 

00 24 1*1, NN 

UTMAK«QTMAX»010 

IF lOTMAX.Gf .IOMAX) GO TO 30 

CALL OAtlt I QUID, l SUPP ,L OUT,ORHOL tORKOV,OKLAM,OSIGM,OXMUL ,U)( 

1 UV, (JIM AX, PVOMA.XMw, GAMMA) 

If l LOUT .GT, 1 ) GO TO IS 

CAU OAC(LIUUlO,LSUPP,LOOT,URHOt,0»HOV,OKtAM,OSlGM,OKMUL,UA 

1UV , I ON JN* PVOMI ,X MW, GAMMA) 

IF tLOUT.GT.n GO ra 15 

VTO*tTVMAX-TVMlN)*DOTV 

NNN*0 

DO 31 1 , MM 

VTD*VTC*OOTV 

wtmax«t VN*y TD/2.0 

VTM(N*TVN-V10/?.0 

IT 1 VTMIN.EO* VTMAX I GO TO 12 

CALL OAC( LIQUID, L*“PP, LOOT, ORHOl,ORHOV,0*L4N,05lGM,OXMUlfOX 

1UV» VTMAX PVMAX, XMM, GAMMA ’ 

IF (LUUT.GT.1I CO TO 35 

CAU DAClLlQU!D,lSUPP,t{H!T,OAHOL,ORMQV t OjaA*,OStG*,OXNUl»DX 

1UV,VTMJN,PVM IN, KMM, GAMMA) 

IF (LOUr.GT.U GO TO 35 

VCDVR-( (PVMAX-PVOMA)*TOM|N/| IRVHIN-PVUMi IRQTrAXJ-l.O) 

vrovc* l.o/vc dvr 

XMOVC« ( l. ♦VRDVU*! PVMIN*PVOMUF( rOMNRRG) 

OUT \\ J>«VTO 

OUT 2 ( J ) ■ VRDVC 

OUT >1 J)*XMDVC 

IF ( vkpVC.l T.o. 0) GQ TO 52 

IF ( VRDVC. Gt.VRVC) GO TO 32 

NNN-J 

31 CCNf (NUt 

32 IF (NNN.EQ.Ot GO TO 99 
NNP-*NNK/Z 

00 3 3 K**UNNP 

STOR)*OUniX) 

SIQRZ*0UT2(K> 

STDR3-0UMIK) 

NP*NNN*1-K 
UuTl *KI«OUTl(NPI 
UUT2IK )-UUT2(«P) 

UUf3(* l-OUTKNP) 

UUTKNPI-STORl 
OU T 2 ( NP l * S T OR ? 



le-ro 


{ 

k 


UUT3lhP)-STOR3 ^ 

Jl CONTINUE 

99 WRITE 16,111 Ct [ 

WRtlfc (6,26) t 

■RITE 16.2?) 

WRITE (6,26) OTNAX : ' 

IF (NNh.NE.QI CO TO 52 

WRITE (6,53 1 hr; 

CO TO 29 

52 WRITE (6,361 I OUT III 1 1 ,0UT2U 1 1 .OUT 31 1 1 1 , 1 1 -l ,NMN1 ^ 

29 CONTINUE . 

30 IF UNO.NE. 1) CO TO 35 - 

18 WRITE' (6,11 I 

NR I It 16,36) 

WRITE (6,20) 

WRITE (6,1211 
WRITE (6,122) TOMAX,TOMIN 
WRITE 16,371 
WRITE (6,26) QNAX 

ctvna-tvnax 

CTVH1-TVH1N 

96 GHAX-QMAX/ | CTVNA-TOHAX ) 

CALL DACILlQU(0,lSUPP,LOUT,ORHOl,ORHQV,aXLAH,OSIGH,OXHUL,OX 

1UV,CTVWA,PVNAX,XMW, GANNA) 

IF (LOUT.GT.t) GO TO 35 

CALL OACIL1QUIU,15UPP,LOUT,ORHOL,ORHOV,OXLAM,OS1GN,OXNUL,OX 

1UV,CTVM1,PVH1N,XNW,CANNAI 
IF ILOUT.GT.I) GO TO 35 

CALL nACUUUID,LSUPP,lOUT,ORHQL.OkHOV,OXLAN,OSIGM t OXNUL,OX 

1UV,T0NAX,PV0NA,XNW, GANNA) 

IF (LOUT.GT.il GO TO 35 

CALL DAC(LlQUID,LSUPP,LOUT,ORHOL,URHOV,OXLAN r OSIGH t OXNUL,UX 

1UV.T0NIN,PV0HI,XNW,GANHA) 
if ILOUT.GT.I) GO TO 35 

VCDVR»JPVMAX-PVONAI*TON1NM(PVNIN-PVOM116TOHAX)-UO 

VROVC-l.J/VCOVR 

IF (VROVL.GT.O.O) CO TO 56 

CrVPA-CTVNA*1.0 

CTVNI«CTVNl-l.O 

XOUNT»F?UNT ♦ 1 

IF IKUUNf.GT.20) CO 10 135 
GO TO 56 

J35 WRITE (6,53' 

GO TC 35 

54 WRITE 16,1231 

WKlTt (6,122) CTVHA.CTVH! 

XMOVC* ( l .0* VRDVC I •( PVN1 N-PVUH1 1/ ( TON !N*RG ) 

IF (KOUNT.EO.O) GO TO 253 
■KITE (6,254) 

253 WRITE (6,11) 

WRITE (6,26) 

WRITE (6,38) 

■KITE (6,50) VKOVC 
WRITE (4,391 
wPITi 1 6, 5*' I XNOVC 
70*7 uN] N-UTU 
:iA«(H.NAX-inN|S)/DTU*l 
'jr. M |-1, NA 



Nh-(CTVNA-CTVHII/0TV41 

Tu>flt*OTO 

KT -0 

IF (IO.GT.TONAX) GO TO 35 

CALL OACaiUUlO,LSUPP,LOU?,OflHOl t ORHOV.OXLAN,OSIGN,OXNUL,OX 

LUV, TO, PVO.XMW, GANNA) 

IF ILLUT.CT.1I CO TO 35 
WRITE (6,11) 

WRITE (6,40) 

WRITE (6,28) TO 
1V*CI VNl-DTV 
DO 141 J-I.N8 
TV»T V*DTV 

IF ( TV.CT.CTVNAI CO TO 141 
IF ( TV.EQ.TOt GO 10 141 

CALL DAC(LlflUIO,LSUPP,LOUT,ORHOL,ORHOV,OXLAN,OSICN,OXNUL,OX 

1UV, TV, PV,XMW, GANNA) 

IF (LOUT.GT.I) tO TO 35 

F-1.0*VR0VC-RG*XH0VC*T0/(PV-PV0I 

U»GHAX*F*ITV-Tui 

IF I U.LT.O. 0) CO TO 141 

TV1-IV 

K T ■ K T ♦ 1 

IF ( KT.GT .6) GO TO 201 
IF (F.GT.1.0) GO TO 149 
141 CONTINUE 

201 TV-CTVNI-DTV 
DO 202 L-l.NB 
TV-TV*0TV 

IF (TV.GT.CTVNAI GO TO 41 

IF (TV.EQ.TOl GO TO 202 

CALL dac(liuuio,lsupp,lout,qrool,orhov,oxlan,osicn,oxnul,ox 

1UV, TV. PV.XNW, GANNA) 

IF (LOUT.GT.U GO TO 35 
F«WD*VRDVC-RG*XN0VC*T0/(PV-PV0> 

Q»GNAX*F*(TV-TU) 

IF lG.LT. 0.01 CO TO 202 

If (F.GT.1.0) CO T r J 41 

WRIU ( 6 , 42 ) TV, 0 ,F 

202 CONTINUE 
GO TO 41 

149 NfW*.0*0TV 

?V«m-6.0*DTV 
DO 152 K* 1 ,N6 
IV»TV*DTV/“.0 
If ITV.EO-TOI CO TO 152 

CALL DAC IL 1 WUI 0 ,L SUPP, LOUT ,nRHOL,ORHOV, OXL AN, 0 SIGN, OXNUL,UX 

IUV» W fPV.XNW, GANNA) 

IF (lOUT.GT.il GO TO 35 

F- I . J*VROVC-RG*XROVC*TO/(PV-PVO) 

U-l,HAX*F*(TV-TU) 

IF (y.Lf.C.O) GO TO 152 
If (F.GT.1.0) GO TO 41 
WRITE (6,421 T V ,Q« F 
152 CONTINUE 
41 LUNIlNUt 

14 IE (MORE. GT. || GU TO 51 
>T;.P 


Appendix D. Sample Problem 



The sample problem consists of determining storage volume requirements 
(Design Analysis) for a given set of sink conditions and of predicting the performance 
(Performance Analysis) within the range of these sink conditions and within a maximum 
acceptable range of vapor temperatures. The heat pipe working fluid is ammonia at a 
nominal operating temperature of 45°C, and the noncondensable gas Is nitrogen. The 
sink conditions and maximum control range are specified as: 


T 

o, max 

308°K 

T . 

o, mm 

268°K 

T 

v, max 

328°K 

T 

v, mm 

308°K 


The heat that must be dissipated at maximum conditions (i. e. , T and T ) 

v, max o, max 

is specified as: 


Q = 6 w 


and the maximum allowable reservoir to condenser vapor volume ratio is: 


V 



v,c 


The calculation increments are specified as: 

AT = 1°C 

v 

A(AT v ) = 1°C 

AT = 10°C 

o 

The associated data cards are listed in Table C. 3.D-1. The resulting computer output 
data follow. The results are also plotted in Figures C. 3. D-i and C. 3. D-2, In Figure 
C. 3. D-l, the storage requirements are shown as a function of control sensitivity for 
various sink temperature ranges. As indicated in the computer printout, the cold- 
reservoir cannot provide control with the specified volume ratio at maximum sink 


servoir-fo -Condenser Ratio V r /V 











Operating Vapor Temperature T, (°K) 


1 

S^Kjempera 

lure 3g8 0 K!^jL° ^ 


298°Kt77^Fl 




268*K(23 0 _F> 


Heat Dissipated Q (watts) 


Figure C. 3. D-2. Performance Analysis: Operating Temperature 

Versus Heat Dissipated 


Nominal Vapor Temperature 
ttangc of Sink Temperatures 
Working Fluid 


318°K (113°F) 

2G8°K - 308°K (23°F - 95° F) 
Ammonia 
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Table C* 3. D-l List of Associated Data Cards 


AMMONIA 

NITROGEN 

SCONTRO MORE=2fLI0UID*:5tLCl = 2fLC2*lS 

SDATAIN TOMAXmlOS. Ot TOM IN*268.0fDTO* 10-0 » T VM AX«328 . 0 1 T VM IN*308 • 0 • 
OTV=U0tDOTV*1.0tRG»296 # 59fQMAX*6 # 0tVKVC*50.0S 
AMMONIA 
NITO^GEN 

5C0NTR0 MORE=0*LCl=lfLC2a2$ 


temperatures above 298°K. The results from the Performance Analysis are given in 
Figure C. D-2 which shows the heat pipe operating temperature T y versus the heat 
dissipated Q at different sink temperatures. The allowable vapor temperature con- 
trol range had to be adjusted in order to obtain control with a cold- reservoir at the 
specified extreme sink temperatures. This adjustment is indicated on the first page 
of output data for the Performance Analysis. Also, the calculation increment for the 
vapor temperature is adjusted internally so that at least six data points are calculated 
for a given sin!; temperature. This is to guarantee that sufficient data is available for 
curve plotting. 
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C. 3-28 


VARI ABLE CONDUCT ANCE HEAT PIPE A N a «. Y s | SIVCHPaI 


ICOL p-A lCKE 0.. R L-?-£ A-Y_u J-JLL-- 


R0RKIN6 FLUID AMHONJA 


NON-CONDENSIBLE CaS NlTMOfiLN. 


CAS CONSTANT IN-M/KQ-k) ?VA.SVQ 


-a. 




L>> 


DCSIbN ANALYSIS 


INPUT D AT A 


,3? 

.9 y 




MAXIMUM ALLqRAMLE SINK TEMPIK) MlNl«UH 5t H K TEHp J K) 


300*00 


240*00 


-Haiti MAL-QP£K A T U tfr—lLHE I K ) HAX. OPERATING TEnPIkI M | N. OPERATING TCnPCkI 

>10*00 >21*00 : Iflfl" SQ 


MAXtHUM ALLORABIE RaT|0 Of NcSeRVOIR VOLUHe TO CONDENSER VOID 


SO. 0000 


fp 


U 1 PUT DATA 

HU. UIUL-Uilf IXL 

**••0 00 

1*PML .... 

*•000 

Ufl/VC) 

— 

• l*Ull*Oi 

lillllMOl 

•2*1004*03 

i 



— 1-iDttB. 

i»apa_ _ 

...... *Jsfi44*JlL 

•4rf«jfi*gi 

; (I7lll»*0) 

• 1 JQf 4 | *P3 

f 1 (11124*9] 

„ _ 

Ax ODD- 

jX4*h?*qi 

.•■07jf*u7 


— 

iiua 

*12*441*01 ... 

*7Jp*2iq*07 

i 

— 

ItOOfl 

u»ui*ai. . 


j 



. .. »«ooo.. 

_ *mm*oi _ 

•4ftf4***02 

i 


iJJfOOO, _ 



*121M1*02 

> 


IltOOO 

__ __ (111171*01 

• IHfll 1*1*02 

; 

, 

ll»O00_ 


•40Ol|0«O2 

\ 


l*.000_ 

. _ • |M407*Q| 

•40*42**0? 









c ly; ; . 

1 4.000 

'-Ul*l r >tSll 

0717*1*07 v 



. 1 % .000 


•11*072*0? 

i 

. 

. ._l*»ooo 

• 1.1* | 44*0J 

.3271 3*»*02 


- 1 

17.000 

iU1*H*UI .. . 

• JJ4*I**02 

* 

— - 

If *000 


•21*43**02 


- 

If.nuo 
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